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ABSTRACT 


The  objective  of  this  thesis  was  to  contribute  to  the  development  of  a  second-order 
approximation  to  the  steady  and  oscillatory  lifting  flow  past  low  aspect  ratio  wings  by  establishing 
the  validity  of  a  NASA-Ames  developed  incompressible  three-dimensional  flow  panel  code  named 
PMARC,  which  would  subsequently  be  used  to  test  the  range  of  validity  of  this  second-order  theory. 
The  steady  state  configurations  modelled  in  the  validation  process  include  swept  back,  delta  and 
rectangular  wings,  an  F5  wing  and  three  equivalent  bodies  of  revolution.  Oscillatory  cases  were  also 
run  with  one  delta  and  the  F5  wings  and  their  respective  equivalent  bodies  of  revolution,  as  well  as 
with  one  spindle.  In  steady  state,  comparison  with  experimental  and  theoretical  data  proved 
PMARC  to  be  very  accurate  for  lift  and  pressure  calculations,  but  revealed  a  discrepancy  in  the 
velocity  distribution  calculation  around  delta  wings.  This  finding  was  corroborated  by  applying  the 
slender  body/slender  wing  theory  (Oswatitsch-Keune  theory)  to  the  delta  wing.  The  unsteady  state 
results  are  presented,  but  their  validation  is  left  for  future  work.  As  part  of  its  primary  objective,  this 
thesis  also  presents  a  computer  code  that  generates  the  F5  wing  equivalent  body  of  revolution  from 
its  chordwise  section  definition.  In  order  to  be  used  as  base  software  to  the  second-order  theory, 
PMARC  will  have  to  be  corrected,  or  a  new  software  will  have  to  be  validated. 
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I.  INTRODUCTION 


Oyer  the  past  quarter  century,  great  advances  have  been  made  in  the  computation  of  the 
aerodynamic  forces  on  aircraft  and  missiles.  These  advances  were  made  possible  by  the 
development  of  high-speed  computers  which ,  in  turn,  made  possible  the  numerical  integration  of 
the  compressible  flow  equations.  However,  the  numerical  computation  of  the  flow  field  around  a 
typical  military  aircraft  or  missile  configuration  is  still  a  very  time  consuming  and  expensive  task 
because  it  requires  the  development  of  a  complex  three-dimensional  computational  grid  and  the 
solution  of  a  system  of  nonlinear  partial  differential  equations  by  means  of  a  suitable  numerical 
integration  scheme.  Furthermore,  the  design  of  aircraft  or  missiles  requires  not  only  the 
determination  of  the  aerodynamic  forces  in  steady  flight  but  also  of  the  transient  and  oscillatory 
forces  in  order  to  assure  the  flight  vehicle's  dynamic  and  aeroelastic  stability  and  to  predict  the 
dynamic  loads  due  to  wind  gusts.  For  these  reasons  there  exists  a  continuing  need  to  develop 
approaches  which  are  sufficiently  accurate,  yet  are  computationally  more  efficient  than  full-fledged 
three-dimensional  solutions  of  the  inviscid  compressible  flow  equations  (Euler  equations)  or  of  the 
viscous  compressible  flow  equations  (Navier-StokeS  equations). 

It  is  the  objective  of  this  thesis  to  contribute  to  the  development  of  a  new  approach  to 
determine  the  steady  and  unsteady  airloads  on  missiles  and  on  low-aspect  ratio  aircraft 
configurations  typically  used  for  high-performance  fighter  aircraft.  This  approach  is  based  on  an 
extension  of  the  well-known  slender-body/slender-wing  theory  developed  by  Oswatitsch,  Keune, 
R.T.  Jones  and  others  in  the  1940's  and  1950's  (see  e.g.  the  book  by  Ashley  and  Landahl,  Ref.  [8]). 
Oswatitsch  and  Keune  showed  that,  in  a  first  approximation,  the  non-lifting  flow  over  bodies  of 
revolution  and  low  aspect  ratio  configurations  can  be  obtained  by  the  superposition  of  two  flows, 
namely  the  flow  in  the  planes  normal  to  the  fteestream  direction  (the  cross-flow)  and  the  "spatial 
influence"  which  accounts  for  the  contribution  upstream  and,  in  subsonic  flow,  also  downstream  of 
the  cross-flow  plane  under  consideration.  The  cross-flow  can  be  treated  as  a  purely  incompressible 
flow  (even  though  the  missile  or  aircraft  flies  at  transonic  or  supersonic  speeds),  whereas  the  spatial 
influence  is  Mach  number  dependent  but  has  the  remarkable  feature  that  the  spatial  influence  is  the 
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same  for  all  configurations  which  have  the  same  cross-sectional  area  distribution  in  the  streamwise 
direction.  The  lifting  flow,  on  the  other  hand,  can  be  obtained  by  merely  considering  the  cross-flow, 
hence  reducing  the  analysis  to  a  purely  incompressible  flow  analysis  for  high-speed  flight  vehicles 
of  low  aspect  ratio. 

M.F.  Platzer  has  published  a  NASA  technical  memorandum  (NASA  TM  X-53295,  April, 
1965)  on  his  development  of  a  second-order  approximation  to  the  steady  and  oscillatory  lifting  flow 
past  low  aspect  ratio  wings  which  exhibits  the  same  fundamental  division  into  an  incompressible 
cross-flow  and  a  Mach-number  dependent  spatial  Influence.  To  test  the  range  of  validity  of  this 
second-order  theory  it  is  necessary  to  compare  it  with  an  exact  three-dimensional  flow  solution. 
Fortunately,  such  a  solution  has  become  available  for  low  speed  (incompressible)  flow  in  recent 
years  with  the  NASA-Ames  developed  PMARC  code.  Therefore,  in  this  thesis,  the  following 
specific  objectives  were  attempted,  namely: 

1 .  Demonstrate  and  document  the  PMARC  code, 

2.  Compare  the  PMARC  solutions  with  other  available  solutions  in  order  to  establish 
its  validity, 

3 .  Compute  the  incompressible  flow  over  a  delta  wing  at  zero  angle  of  attack  using  the 
Oswatitsch-Keune  theory, 

4.  Obtain  the  equivalent  area  body  of  revolution  of  an  F-5  wing, 

5.  Compare  the  Oswatitsch-Keune  theory  with  the  PMARC  predictions,  and 

6.  Generalize  the  non-lifting  theory  to  incompressible  flow  past  oscillating  bodies  of 
revolution  and  low-aspect  ratio  wings. 

These  tasks  will  be  described  in  the  subsequent  chapters. 
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II.  SOFTWARE  OVERVIEW 


A.  THEORY  BEHIND  PMARC 

PMARC  stands  for  Panel  Method  Ames  Research  Center.  It  is  a  potential  flow  panel  code 
developed  at  NASA  Ames  Research  Center  to  numerically  predict  flow  fields  around  complex 
three-dimensional  geometries.  The  general  idea  behind  this  panel  code  is  to  cover  the  body  under 
study  with  panels  over  which  there  is  an  unknown  distribution  of  doublets  and  sources.  (Recall  that 
the  doublet  and  source  flows  are  two  of  several  elementary  flows  used  to  model  the  flow  field 
around  arbitrary  shapes.)  These  unknovras  are  solved  through  a  set  of  simultaneous  linear  algebraic 
equations  generated  by  calculating  the  induced  velocity  at  each  panel's  control  point  and  applying 
the  flow  tangency  condition. 

Two  kinds  of  simplifications  were  made  to  allow  this  software  to  remain  computationally 
cost  effective.  The  first  one  has  to  do  with  the  strength  of  the  elementary  flows  (also  called 
singularities).  In  PMARC,  the  strength  of  the  singularities  is  constant  over  each  panel  allowing 
factorization  of  the  doublet  and  source  strengths  in  the  thus  greatly  simplified  analytical  solution  of 
the  flow.  The  second  kind  of  simplification  is  in  the  assumption  of  an  inviscid,  irrotational,  and 
incompressible  flow  field,  thus  satisfying  Laplace's  equation: 

V2$=0  (2.1) 

With  these  simplifications,  PMARC  proceeds  to  solve  for  the  potential  at  any  point  P  by 
applying  Green's  Theorem,  which  results  in  the  following  integral  equation: 

[(®-®>-V(4)]  dS  -  [(4)«(V®-V®^)  (2.2) 

4nJJs^w*s„  r  4tcJjs*w*s„  f 
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where  r  is  the  distance  from  the  point  P  to  the  element  dS  on  the  surface  and  n  is  the  unit  vector 
normal  to  the  surface  pointing  into  the  flow  field  of  interest.  In  this  equation  the  first  integral 
represents  the  disturbance  potential  from  a  surface  distribution  of  doublets  with  strength  'phi-phi,' 
per  unit  area  and  the  second  integral  represents  the  contribution  from  a  surface  distribution  of 
sources.  As  shown  on  Figure  2.1,  phi  corresponds  to  the  external  region  (or  flow  field  of  interest), 
and  phi,  is  the  internal  region  (or  fictitious  flow).  This  arrangement  is  required  to  model  external 
flow  problems  such  as  a  wing  in  a  uniform  stream. 


(From  Ref.[l]) 

PMARC  actually  solves  a  simpler  form  of  Green's  Theorem  than  the  one  presented  above. 
The  simplifications  are  made  possible  by  noting  that  at  the  surface  at  infinity,  the  perturbation 
potential  due  to  the  configuration  is  essentially  zero,  leaving  only  the  potential  due  to  the  uniform 
onset  flow.  Also,  it  is  assumed  that  the  wake  is  thin  and  there  is  no  entrainment  so  the  source  term 
for  the  wake  disappears  and  the  jump  in  normal  velocity  across  the  wake  is  zero. 
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The  theory  described  above  makes  for  three  major  steps  in  the  PMARC  solution  process. 
First,  it  solves  for  the  doublet  strengths  leading  to  known  values  of  the  singularities  on  each  of  the 
panels.  Second,  it  computes  the  velocity  by  differentiating  the  velocity  potential  expression  in  the 
appropriate  direction  for  each  tangential  component  of  velocity.  Using  this  information  it  calculates 
the  resultant  velocity  at  each  panel  control  point.  With  this  result  it  finally  derives  the  pressure 
distribution  using  Equation  (2.3).  The  forces  and  moments  on  the  body  are  then  obtained  by 
performing  a  surface  pressure  integration. 


A? 


(2.3) 


B.  PMARC  OPERATING  MODE 

Instead  of  giving  all  the  details  pertaining  to  the  operation  of  the  PMARC  software  within 
this  section,  only  the  generalities  will  be  presented.  For  practicality  and  ease  of  understanding,  the 
information  and  instructions  for  construction  of  a  specific  wing  configuration  will  be  presented  in 
the  software  validation  sections,  each  time  a  new  configuration  is  introduced.  For  additional  details, 
the  reader  should  refer  to  Refs.  [1],  [2],  and  [3].  Installation  instructions  are  also  found  in  Ref.  [1]. 

1.  Generalities 

There  exists  more  than  one  way  to  construct  a  structure  in  PMARC.  One  is  to  make  use  of 
three  levels  of  coordinates,  as  shown  in  Figure  2.2 ,  in  order  to  join  components  built  separately. 
Another  way  would  be  to  build  the  whole  structure  at  once  referring  all  points  to  one  coordinate 
system. 


5 


Figure  2.2.  PMARC  Hierarchy  of  Coordinate  Systems. 
(From  Ref.  [2]) 


But  no  matter  which  method  is  preferred,  the  physical  representation  of  each  component 
is  made  with  sheets  of  panels  called  patches  used  to  solve  the  required  simultaneous  flow  equations 
for  that  particular  group  of  panels .  The  panel  density  of  these  patches  and  the  manner  in  which  all 
the  patches  are  joined  together  are  defined  within  the  input  file  as  will  be  shown  in  later  sections. 
The  shape  of  theses  patches  is  given  by  providing  a  point  definition,  or  path,  around  or  along  which 
the  patches  will  be  laid  down  or  wrapped.  Figure  2.3  shows  how  a  patch  is  rolled  up  on  top  of  itself 
to  make  a  surface  in  which  two  of  the  four  sides  are  joined  in  a  seam.  There  are  a  few  rules  about 
the  numbering  of  each  panel,  and  although  most  of  them  are  transparent  to  the  users  upon 
construction  of  a  model,  others  shall  be  known  for  output  interpretation  purposes.  These  rules  will 
be  explained  later. 


6 


Figure  2.3.  PMARC  Patch  of  Panels. 
(From  Ref.  [2]) 


The  previous  figure  is  a  representation  of  hoAv  wings  are  modelled.  Again,  PMARC  offers 
some  flexibility  as  there  are  two  ways  to  define  a  wing  section.  The  first  one  is  done  by  entering  a 
series  of  points  that  define  the  wing  cross-sectional  shape.  PMARC  will  then  wrap  the  patch  on  top 
of  that  definition.  The  second  method  will  automatically  generate  the  point  coordinates  of  a  NACA 
Four  Digit  Series  airfoil.  All  that  is  required  are  the  airfoil's  maximum  thickness  ratio,  and  the 
maximum  camber  and  its  location.  A  later  section  will  show  how  to  use  these  two  options. 

Once  the  shape  of  a  section  is  defined,  the  third  dimension  of  the  component  has  to  be 
modelled.  Again,  several  options  are  available  to  this  end.  One  is  to  position  the  leading  edge  origin 
of  each  end  of  the  component,  for  instance  the  root  and  the  tip  cross-section  of  a  wing,  in  some 
chosen  coordinate  system,  and  then  to  scale  each  of  these  ends  to  the  desired  size.  Another  method 
is  to  define  the  cross-section  of  each  end  of  the  component  in  relation  to  a  coordinate  system.  This 
way,  each  cross-section  is  built  originally  to  its  desired  size,  instead  of  being  scaled  from  the  non- 
dimensionalized  definition  as  with  the  previous  method.  Depending  on  the  application,  each  of  these 
methods,  or  a  combination  thereof,  show  significant  advantages. 


PMARC  incorporates  routines  to  automatically  generate  wing  tips,  whether  they  need  to  be 
flat  or  rounded.  Also,  by  toggling  the  value  of  certain  variables,  complete  components  can  be  copied. 
A  multitude  of  other  features  are  available  among  which  is  found  the  wake  generation  and  the 
streamline  generation  routine.  The  interested  reader  should  refer  to  Refs.  [1],  [2],  and  [3]. 

2.  PMARC  Input  File 

Several  decisions  have  to  be  made  when  building  a  geometry  with  PMARC.  The  number 
of  panels  spanwise  and  chordwise  and  their  distribution  are  two  important  ones.  More  panels  means 
better  accuracy,  but  also  means  proportionately  longer  computational  time.  The  optimal  number  of 
panels  should  be  sought  to  make  the  results  sufficiently  accurate.  The  panels  can  be  spaced 
following  four  types  of  distribution:  full  cosine,  two  types  of  half  cosine,  and  equal  spacing.  These 
allow  the  user  to  concentrate  more  panels  in  areas  of  expected  greater  gradients.  These  decisions  will 
be  justified  in  the  validation  section  each  time  a  new  configuration  is  tested. 

Other  information,  geometry  specifics,  and  test  conditions  are  entered  into  PMARC  via  its 
input  file.  It  is  called  DATA5  and  is  FORTRAN  based.  Over  200  fields  or  variables  are  required 
to  be  filled  for  PMARC  to  successfully  read  the  file.  Although  not  all  the  fields  are  pertinent  to  all 
cases,  an  uncompleted  field  will  cause  a  software  crash.  The  following  paragraphs  will  laconically 
describe  the  DATA5  input  file,  identifying  which  fields  should  be  filled  with  correct  values  to  obtain 
correct  output  in  the  case  where  we  are  interested  in  the  Cp,  Vx,  Vy,  Vz,  V,  CL,  and  Cl  but  without 
wake  or  streamline  generation.  The  reader  should  refer  to  Figure  2.4  before  continuing.  Also,  each 
field  of  the  DATA5  input  file  is  represented  by  an  acronym.  The  definition  of  these  acronyms  along 
with  the  valid  entries  are  listed  in  Ref.  [1]  starting  on  page  33. 
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Figure  2.4  PM  ARC  Input  File  Sample. 


&VS3  Xl=  -0,1000,  YU  1.5000,  ZU  0.1000,  MPTU  0, 

iiVS-l  X2=  -0.1000,  Y2=  1.5000,  Z2»  0.1000,  MPT2=  20, 


a.  Basic  Input  Section 

This  section  allow  the  user  to  entitle  the  file,  prescribe  specific  job  requirements, 
control  some  solver  and  time-step  parameters,  define  some  of  the  testing  conditions  as  far  as  the 
movement  of  the  model  and  the  ffeestream  are  concerned,  and  also  to  input  reference  dimensions. 
All  the  above  make  up  for  sixty  fields.  Not  all  of  them  require  actual  testing  values  in  order  to 
obtain  valid  output  for  the  parameters  listed  in  subsection  2  above.  These  fields  can  be  filled  with 
any  number  without  altering  the  results.  In  the  steady  case,  the  critical  fields  are  the  type  of  run 
(LENRUN)  that  determines  if  a  full  flow  calculation  will  be  performed  or  just  a  geometry 
generation,  the  angle  of  attack  (ALDEG),  and  the  reference  surface  area  (SREF),  the  latter  being 
used  in  the  local  lift  coefficient  calculations.  Other  fields  become  critical  when  dealing  with 
unsteady  state  cases.  They  will  be  discussed  later  in  the  oscillation  chapter.  Two  of  these  fields 
(NTSTPS,  DTSTEP)  also  require  some  attention  in  the  steady  state  cases  as  will  be  seen. 

One  should  always  study  the  output  carefully  as  some  of  the  results  might  be  the 
consequence  of  fields  that  were  filled  with  arbitrary  values  instead  of  values  corresponding  to  the 
test  case.  A  prudent  habit  for  the  analysis  of  any  kind  of  wing  or  wing-body  combination  fixed  in 
space  is  to  fill  the  following  list  of  fields  with  the  information  corresponding  to  the  object  under 
study:  LENRUN,  VINF,  VSOUND,  ALDEG,  CBAR,  SREF,  and  SSPAN.  The  other  fields  can 
remain  the  same  as  the  values  used  in  example  files  presented  in  the  PMARC  manual  (Ref.  [1]). 

b.  Input  Section  for  Surface  Geometry 

This  portion  of  the  input  file  includes  the  headings  &ASEM  1 ,  &ASEM2,  &COMP  1 , 
&COMP2,  &PATCH1,  &SECT1,  and  &BPNODE.  For  the  example  analysis  described,  the  first  four 
sets  of  fields  don't  change  from  what  is  shown  on  Figure  2.4.  They  will  change  if  two  or  more 
independent  components  are  studied.  The  last  three  sets  (&PATCH1,  &SECT1  and  &BPNODE)  are 
critical  in  the  construction  of  components  and  the  definition  of  its  panel  layout.  Each  individual 
field  within  these  sets  has  to  be  entered  with  values  specifically  corresponding  to  the  model  being 
built  and  the  conditions  of  the  tests.  As  will  be  seen  later,  the  sets  of  fields  listed  above  are 
duplicated  as  often  as  required  to  define  more  complex  components  or  assemblies.  Also,  &PATCH2, 
&PATCH3,  &SECT2,  and  &SECT3  are  other  sets  of  field  required  in  certain  cases. 
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It  becomes  a  complex  matter  to  identify  which  cases  require  which  ones  of  the  sets 
of  fields  discussed  above.  The  uninitiated  should  study  the  various  test  cases  in  the  validation 
sections  to  acquire  a  broader  view  of  the  many  possible  applications. 

c.  Time-stepping  Wakes  and  On-body  Streamlines 

Wake  analysis  with  PMARC  is  a  challenging  undertaking.  If  wakes  are  to  be 
modelled,  each  of  the  fields  in  the  wake  portion  of  the  input  file  is  important.  These  fields  lose  all 
significance  when  the  wake  routine  is  turned  off  (IDWAK  =  0).  The  same  comments  are  true  for  the 
on-body  streamlines. 

3.  PMARC  Output  Files 

Once  the  PMARC  has  successfully  read  the  input  file  DATA5,  it  performs  the  flow 
calculations  and  generates  two  output  files  named  DATA6  and  DATA22.  The  latter  is  an  unreadable 
file  that  contains  all  the  information  required  for  the  visualization  of  the  various  parameter  gradients 
along  the  body  under  study.  This  visualization  is  made  possible  via  a  satellite  software  named  GVS. 
General  information  on  GVS  can  be  found  in  Ref  [3]  and  a  detailed  user  guide  is  also  available  in 
Ref  [18]. 

DATA6  is  the  output  file  formatted  for  interpretation  by  the  user.  Its  length  could  range 
fi’om  a  few  kilobytes  for  the  simple  geometry  computation  to  several  megabytes  for  full  flow 
computation  using  several  time-steps.  DATA6  is  a  well  formatted  and  structured  output  file.  The 
following  provides  a  succinct  description  of  its  structure. 

a.  Input  list  for  Basic  Data 

The  top  of  that  file  reiterate  some  of  the  important  geometry  and  test  conditions 
prescribed  by  the  user  within  the  input  file,  DATA5.  It  also  lays  out  the  position  of  the  assembly  and 
component  coordinate  systems.  This  portion  of  the  file  ends  with  the  cross-sections  definitions  and 
a  list  of  the  basic  patch  data. 
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b.  Time-step  zero  solution 

For  the  first  moment  of  the  test,  the  following  output  will  be  produced: 

(1)  Basic  wake  information: 

(2)  Solution  iteration  history: 

(3)  Aerodynamic  data  for  patch  one  which  consists  of  a  series  of  tables 
that  lists  for  each  of  the  panels  in  one  chordwise  row  the  x,  y  and  z  coordinate  of  the  panel  control 
point,  as  well  as  Cp,  Vx,  Vy,  Vz,  V,  DUB,  DPHI/DT,  and  MACH.  This  process  is  repeated  for  each 
row  of  panels  of  a  patch. 

(4)  Section  parameters  for  the  series  of  rows  of  panels  which 
aerodynamic  data  were  just  listed.  The  section  parameters  are  CL,  CD,  CY,  Cm,  Cn,  and  Cl.  These 
values  can  be  regarded  as  the  two-dimensional  values  for  that  strip  of  panels.  They  are  given  in  two 
reference  sets:  wind,  and  body  axes. 

(5)  Aerodynamic  data  and  section  parameters  for  all  other  patches.  The 
same  format  as  presented  for  the  first  patch  is  followed  for  all  remaining  patches. 

(6)  Wind  axis  force  and  moment  coefficients.  In  this  section,  total  CL, 
CD,  CY,  and  Cm  are  given  for  each  of  the  patches,  each  component,  and  finally  for  the  whole 
assembly. 

(7)  Body  axis  force  and  moment  coefficients.  In  this  section,  total  CN, 
CA,  CY,  and  Cm  are  given  for  each  of  the  patches,  each  component,  and  finally  for  the  whole 
assembly. 


c.  Time-step  1  and  above 

For  the  remaining  time  step  except  the  last  one,  DATA6  will  produce  a  condensed 
output  in  that  it  will  not  generate  all  the  tables  of  aerodynamic  data  as  for  time  step  zero.  Thus,  all 
the  steps  listed  above  are  executed  except  (3)  and  (5). 


d.  The  last  time-step 

For  the  last  time-step,  all  of  the  seven  steps  listed  in  time-step  zero  are  executed 
again.  The  results  will  be  different  than  when  the  flow  just  started  at  time  step  zero. 
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c. 


CONCLUSION 


The  theoretical  basis  for  PMARC,  a  low  order  potential-flow  panel  code  for  modelling 
complex  three-dimensional  geometries  has  been  outlined.  The  input  file  was  also  detailed  with  the 
emphasis  on  the  fact  that  although  not  all  the  fields  apply  to  all  cases,  they  all  have  to  be  filled  with 
values  for  the  software  to  run  successfully.  Finally,  the  output  file  was  reviewed.  Its  format  was 
studied  along  with  the  type  of  results  that  it  holds. 
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III.  SOFTWARE  VALIDATION  -  CONVENTIONAL  WINGS 


A.  INTRODUCTION 

The  objective  of  this  chapter  is  to  show  that  PMARC  is  a  panel  code  that  closely  recreates 
most  of  the  results  from  experimental  testing,  higher  order  panel  codes,  and  theoretical  analytical 
solutions.  As  it  will  be  seen,  there  are  restrictions  to  its  capabilities  in  the  area  of  velocity 
distribution  over  delta  wings  at  zero  angle  of  attack.  This  PMARC  review  will  be  done  by  modelling 
various  types  of  wings  for  which  analyses  have  been  conducted  and  results  are  available. 

Specifically,  two  types  of  swept  back  wings,  four  delta  wings,  various  rectangular  wings, 
an  F5  wing  and  a  parabolic  arc  delta  wing  will  be  studied.  The  lift-curve  slopes,  the  pressure 
coefficient  distribution  in  both  chordwise  and  spanwise  direction,  the  velocity  distribution  and  the 
equivalence  rule  results  will  all  be  comparison  factors  for  our  validation  process.  The  following 
chapter  will  conduct  a  similar  study  on  three  bodies  of  revolutions. 

As  mentioned  in  the  previous  chapter,  PMARC  operating  instructions  will  also  be  presented 
here.  Therefore  the  structure  of  the  following  sections,  each  corresponding  to  a  wing  configuration, 
will  be  as  follows.  A  description  of  the  wing  studied  will  first  be  given.  Then,  an  overview  of  the 
theoiy  and  results  serving  as  the  reference  data  will  be  presented,  followed  by  the  results  obtained 
with  PMARC,  shown  compared  to  the  reference  data,  and  complemented  by  explanatory  comments 
and  discussion.  Also,  the  specific  PMARC  input  steps  for  modelling  the  subject  wing  will  be 
explained.  Finally,  the  output  data  retrieval  technique  is  presented. 

B.  SWEPT  BACK  WING  I 

A  wing  with  relatively  simple  geometry  was  considered  for  the  first  comparison.  The  wing 
has  an  aspect  ratio  (AR)  of  5,  a  taper  ratio  of  unity,  and  an  uncambered  flat  plate  section.  Since  the 
taper  ratio  is  unity,  the  leading  edge,  the  quarter-chord  line,  the  three-quarter-chord  line,  and  the 
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trailing  edge  all  have  the  same  sweep,  45  degrees.  All  the  dimensions  and  coordinates  being  in  terms 
of  the  span  b,  an  arbitrary  value  of  b  will  be  assumed  at  10.  The  finite  geometry  of  the  wing  is  as 
shown  in  Figure  3.1. 


Figure  3.1  Four-Panel  Representation  of  a  Swept  Planar  Wing,  Taper  Ratio  of 
Unity,  AR=5,  Sweep=45  Degrees.  (From  Ref.  [4]) 
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1.  Theory,  Bacl^round  and  Reference  Data 

Several  lifting  flow  field  computation  methods  have  been  developed.  Probably  one  of  the 
simplest  is  the  Vortex  Lattice  Method  (VLM).  The  VLM  represents  the  wing  as  a  planar  surface  on 
which  a  grid  of  vortex  panels  is  superimposed.  The  velocities  induced  by  each  horseshoe  vortex  at 
a  specified  control  point  are  calculated  using  the  Bio-Savart  law.  Equation  (3.1). 

—  r  jc 
dv  =  — 

47tr^ 

A  summation  is  performed  for  all  control  points  on  the  wing  to  produce  a  set  of  linear  algebraic 
equations  for  the  horseshoe  vortex  strengths  that  satisfy  the  boundary  condition  of  no  flow  through 
the  wing.  The  vortex  strengths  are  related  to  the  wing  circulation  and  the  pressure  differential 
between  the  upper  and  lower  wing  surfaces.  Berlin  &  Smith  (Ref.  [4])  present  the  derivation  of  the 
equations  from  the  Bio-Savart  law,  and  apply  them  to  the  case  of  the  swept  back  flat  plate  shown 
above.  The  experimental  results  against  which  the  PMARC  and  VLM  results  will  be  compared  are 
reported  in  Ref.  [5]. 

2.  Results  and  Discussion 

As  shown  on  Figure  3.2,  the  coefficient  results  obtained  by  PMARC  correspond  very  closely 
to  the  ones  calculated  with  the  Vortex  Lattice  Method.  However,  although  good  for  the  low  angles 
of  attack,  both  sets  seem  to  gradually  increase  their  spread  with  the  experimental  results.  That  does 
not  discredit  PMARC  or  VLM,  as  this  behavior  was  to  be  expected.  PMARC  and  VLM  are  both 
linear  methods  that  do  not  account  for  boundary  layer  separation  and  stall.  That  is  why  the  slope  of 
the  experimental  results  shown  on  Figure  3.2  start  to  decrease  as  the  stall  angle  approaches,  while 
the  other  two  methods  keep  showing  a  steady  increase  of  lift  coefficient  with  increase  in  angle  of 
attack. 


r) 


(3.1) 
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Figure  3.2  Comparison  of  The  Lift  Coefficient  as  a  Fvinction  of  Angle  of  Attack. 
Swept  Planar  Wing,  Taper  Ratio  of  Unity,  AR=5,  Sweep=45  Degrees. 


3.  PMARC  Input  Instructions 

The  following  details  complement  Chapter  n  which  presented  the  theory  behind  PMARC 
and  its  major  operating  instructions.  This  subsection  will  reiterate  certain  important  details  and 
introduce  new  factors  that  were  intentionally  left  for  this  section  for  simplicity  reasons.  The 
following  headings  refer  to  the  PMARC  input  file,  DATA5,  shown  in  Figure  3.3. 
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a.  Basic  Input  Section 

When  referring  to  the  PMARC  manual  (Ref.  [1])  for  the  definition  of  each  of  the 
variables,  it  becomes  rather  straightforward  which  values  are  to  be  entered  in  each  field  of  the  Basic 
Input  Section.  As  mentioned  in  Chapter  n,  for  most  of  the  cases,  only  the  following  fields  require 
attention:  LENRUN,  NTSTPS,  DTSTEP,  ALDEG,  and  SREF,  although  we  will  also  continually 
fill  CBAR,  and  SSPAN.  Some  others  like  VINF  and  VSOUND  can  be  entered  with  actual  values 
to  represent  the  fi'eestream  velocity  and  the  altitude  but  these  won't  have  any  impact  on  the  lift 
coefficients  and  pressure  distribution  results.  It  is  to  be  noted  that  the  lift,  normal  force,  moment 
and  drag  coefficients  are  sensitive  to  the  fields  CBAR,  SREF,  and  SSPAN,  as  the  velocity 
components,  the  pressure  coefficient,  DUB,  and  dPHI/dT  are  not.  The  remaining  fields  should  be 
filled  with  zeros  or  numbers  taken  fi'om  existing  DATA5  files.  Although  they  too  don't  have  any 
impact  on  the  above  mentioned  coefficients,  numbers  have  to  be  entered  for  the  software  to  operate. 

One  important  point  is  the  steady  state  status.  Because  of  the  time-stepping 
capability  of  PMARC,  results  are  generated  right  from  the  first  moment  of  the  'flight',  while  the 
circulation  is  initiated  and  the  starting  vortex  generated.  If  the  user  is  interested  in  the  steady  state 
results,  which  is  when  the  starting  vortex  is  far  enough  from  the  trailing  edge  (five  or  six  chord 
length),  a  proper  combination  of  number  of  time-steps  (NTSTPS)  and  distance  between  each  time- 
step  (DTSTEP)  will  be  required.  Using  the  freestream  velocity,  this  calculation  is  easy  to  perform. 
For  instance,  a  wing  with  a  chord  of  1 1  units  would  require  the  starting  vortex  to  be  at  a  minimum 
of  55  units  behind  the  trailing  edge  in  order  to  assume  steady  state.  If  the  freestream  velocity  is  1 1 
units/sec,  then  5  seconds  of  flight  are  needed  (55  units  divided  by  11  units/sec).  Then  any 
combination  of  NTSTPS  and  DTSTEP  whose  product  is  at  least  5  units  will  ensure  steady  state. 

b.  Assembly  and  Component  Coordinate  System  Information 

All  the  fields  but  two  can  be  copied  from  those  shown  in  Figure  3.3  for  most  of  the 
cases  we  will  be  concerned  with.  The  two  fields  that  may  require  special  attention  are  NODEA  and 
NODEC.  Their  values  indicate  respectively  if  another  assembly  coordinate  or  component  coordinate 
will  be  defined.  Such  additional  coordinate  systems  could  be  required  when  modelling  a  wing-body 
combination,  or  two  wings  alone. 
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c.  Patch  and  Section  Coordinate  System  Information 
All  the  fields  in  PATCH  1  are  straight  forward  when  referring  to  the  PMARC  manual 
(Ref.  [1]).  The  point  to  note  here  is  the  existence  of  more  than  one  PATCH  1  section  within  the  input 
file.  One  introduces  the  patch  that  will  be  used  to  create  the  wing  according  to  the  geometric 
specifications  given  in  the  Section  Coordinate  section  (SECTl  and  SECT2),  and  the  second  one  is 
for  the  wing  tip  which  is  further  defined  in  PATCH2. 

In  SECTl,  two  fields  are  not  relevant  in  steady  state:  ALF  and  THETA.  All  the 
others  are  crucial  for  a  correct  model  definition.  In  this  case,  the  origin  of  the  root  airfoil  section 
definition  is  at  coordinates  (0, 0, 0).  The  root  chord  is  two  units.  INMODE  was  set  at  5  because  in 
SECT2  we  will  generate  a  NACA  4  digit  airfoil  section  for  which  PMARC  has  automatic  section 
coordinate  generation  routines.  All  that  is  required  are  the  thickness  to  chord  ratio  (RTC),  the 
maximum  camber  to  chord  ratio  (RMC),  and  the  chordwise  position  of  the  airfoil's  maximum 
camber  (RPC).  In  this  case,  we  are  dealing  with  a  flat  plate.  These  three  combinations  all  produce 
the  same  results: 

RTC=0.01  RTC=0.01  RTC=0.01 

RMC=0  RMC=0  RMC=0 

RPC=0  RPC=1  RPC=0.5 

In  our  case,  the  second  approach  was  taken  which  corresponds  to  the  NACA  designation  of  0001. 
It  becomes  obvious  why  the  three  combinations  shown  above  work  by  noting  that  for  a  flat  plate, 
there  is  no  camber  (RMC),  therefore  the  position  of  maximum  camber  as  described  by  the  input  field 
RPC  has  no  impact,  no  matter  what  value  is  entered. 

The  second  SECTl  section  defines  the  location  and  section  of  the  wing  tip,  or  more 
specifically  that  flat  surface  where  the  wing  tip  would  be  attached.  We  see  that  in  our  case,  the 
origin  of  that  last  patch's  section  is  at  the  coordinates  (5, 5, 0),  and  that  its  chord  is  also  2.  As  defined 
in  Ref.  [1],  INMODE  and  TNODS  seen  in  SECTl  and  PATCH2  respectively  are  important  fields. 
INMODE  set  to  zero  in  the  second  SECTl  tells  PMARC  to  copy  the  section  definition  of  the 
previous  SECTl.  The  TNODS  fields  indicate  the  status  of  the  patch  being  defined:  first, 
intermediate,  last,  etc. 
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Ten  panels  per  side  (TNPS=10)  are  requested  spanwise  for  the  specified  patch.  That 
means  20  panels  spanwise  for  the  complete  wing  on  each  surface.  Fifteen  panels  are  prescribed 
chordwise  (TNPC=15)  for  the  specified  patch.  That  transpose  into  15  panels  on  each  surface  of  the 
wing.  Note  that  TNPS  can  be  specified  only  in  the  second  SECTl  section.  If  both  TNPS  were  set 
to  zero,  PMARC  would  use  the  points  of  the  section  definition  it  created  to  determine  the  number 
of  panels. 


As  mentioned  earlier,  the  second  PATCH  1  section  along  with  PATCH2  are  for  the 
wing  tip,  and  a  circular  arc  tip  patch  was  requested  (ITYP=2). 

Two  fields  were  not  discussed  yet  but  become  very  important  for  the  accuracy  of 
the  results.  They  are  TINTS  and  TINTC.  They  are  found  in  SECTl,  SECT2,  and  PATCH2.  They 
define  the  type  of  panel  spacing  on  each  patch  (or  wing  and  body  surfaces).  As  it  will  be  seen  later 
for  other  configurations,  the  user  might  want  to  distribute  more  panels  in  areas  where  a  sharper 
gradient  of  a  parameter  is  known  to  exist.  Figure  3.4  shows  the  available  options  of  panel  spacing 
for  an  arbitrary  wing. 


d.  Wake  and  Streamlines 

The  streamlines  are  not  required  for  our  studies.  However,  the  wake  should  be 
turned  on  (IDWAK=1)  in  all  wing  cases  as  it  affects  the  flow  field  all  around  the  wing.  For  the 
simple  cases  studied  in  this  thesis,  the  type  of  wake  selected  (rigid  or  flexible)  in  the  IFLXW  field 
is  not  a  crucial  decision  as  this  will  not  change  the  results. 

4.  PMARC  Output  -  Data  Retrieval 

The  results  we  are  interested  in  are  the  lift  coefficients  (Cj,)  for  various  angles  of  attack.  In 
order  to  do  so,  our  input  file  has  to  be  run  repetitively,  each  time  modifying  the  angle  of  attack.  This 
will  create  a  series  of  output  files,  DATA6,  from  which  the  will  be  read.  Several  lift  coefficients 
are  presented  in  this  output  file:  assembly,  component,  and  total,  in  term  of  wind  or  body  axes.  For 
the  present  case,  we  need  the  total  wing  lift  coefficient  in  terms  of  the  wind  axis  which  refers  to  the 
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freestream.  This  value  is  found  towards  the  veiy  end  of  the  file  under  the  heading  Wind  Axes  -  Total 
coefficients.  The  lift  coefficient  corresponding  to  our  case  is  labelled  CL.  Note  that  when 
streamlines  are  requested,  the  related  results  are  placed  at  the  end  of  the  file,  following  the 
coefficient  results. 


Figure  3.4  PMARC  Options  of  Chordwise  and  Spanwise  Panel 
Distributions. 


c. 


SWEPT  BACK  WING  II 


Our  second  project  was  to  model  a  cambered  swept  back  wing.  The  wing  considered  had 
a  NACA  64A010  section,  an  AR  of  3,  a  tip  to  root  chord  ratio  of  0.5  and  a  45  degrees  swept  back 
of  the  quarter-chord  line.  A  value  of  the  root  chord  Cr  was  arbitrarily  set  at  10,  allowing  the  finite 
coordinates  as  shown  on  Figure  3.5. 


Figure  3.5  Swept  Back  Wing,  NACA  64A01 0  Section,  AR=3,  Taper 
Ratio=0.5,  Sweep=45  Degrees. 


1.  Theory,  Background,  and  Reference  Data 

Margason  et  al.  (Ref.  [6])  compared  computed  aerodynamic  coefficients  using  the  Vortex 
Lattice  Method,  source  panel  method,  low-order  surface  potential  distributions,  and  high-order 
surface  potential  distributions  to  experimental  data.  This  provided  an  ideal  opportunity  to  see  how 
PMARC  would  measure  up. 
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2. 


Results  and  Discussion 


The  results  are  shown  on  Figure  3.6.  It  can  be  observed  that  all  surface  potential  and  source 
panel  methods  overpredict  the  experimental  data.  However,  each  of  them  show  little  difference 
between  their  lift  coefficients,  PMARC  included.  The  results  for  the  Vortex  Lattice  Method  are  not 
presented,  but  Margason  reports  that  they  predict  the  experimental  data  very  well,  due  to  the  fact  that 
VLM  neglects  both  thickness  and  viscosity  effects. 


Figure  3.6  Comparison  of  The  Lift  Coefficient  as  a  Function  of  Angle  of  Attack. 
Aft  Swept  Wing,  NACA  64 AO  10  Section,  AR=3,  Taper  Ratio=0.5,  Sweep=45 
Degrees. 


Although  PMARC  did  not  recreate  exactly  the  experimental  data  for  the  NACA  64 AO  10 
swept  back  wing,  its  deviation  from  the  latter  is  still  acceptable.  Furthermore,  its  results  compare 
very  well  with  other  computer  methods,  even  with  panel  codes  of  higher  order. 
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3. 


PMARC  Input  Instructions 


The  following  sub-subsections  will  not  repeat  information  presented  for  the  swept  back  flat 
plate  that  still  hold  for  the  cambered  swept-back  wing.  Instead,  we  will  concentrate  on  the 
differences,  the  particularities,  and  some  lessons  learned.  The  input  file  is  shown  in  Figure  3.7. 

a.  Basic  Input  Section 

All  the  fields  in  this  section  can  remain  the  same  as  for  the  flat  plate  case,  except 
for  CBAR,  SREF,  and  SSPAN.  Of  course,  the  angle  of  attack  (ALDEG)  has  to  be  adjusted  as 
required.  Also,  VINE  and  VSOUND  could  be  changed,  although  there  is  no  need  to  do  so. 

b.  Assembly  and  Component  Coordinate  System  Information 

No  changes  are  required  in  these  two  sections.  As  it  can  be  seen  in  the  definition 
section  of  Ref.  [1],  these  terms  become  important  when  multiple  components  are  involved,  or  when 
a  component  needs  to  be  rotated  or  translated. 

c.  Patch  and  Section  Coordinate  System  Information 

The  PATCHl  sections  for  both  the  flat  plate  and  cambered  cases  are  identical.  The 
PATCH2  lines  have  for  only  difference  the  panel  field  (TNPS).  It  could  have  been  the  same  as  in 
each  case  it  was  an  arbitrary  decision.  This  field  defines  the  number  of  chordwise  (longitudinal)  cuts 
on  the  wing  tip,  which  along  with  the  vertical  cuts  will  define  the  number  of  panels  on  the  tip  (Figure 
3.8).  TNPS  of  PATCH2  can  be  kept  low  without  altering  the  results.  Fewer  panels  speeds  up  the 
computation  and  in  this  case  alleviates  a  message  in  the  output  file  warning  the  user  that  some  panels 
have  virtually  no  area.  Although  this  warning  did  not  seem  to  affect  any  of  the  results,  caution 
should  be  exercised  and  more  investigation  should  be  pursued  in  this  area.  It  is  to  be  noted  that 
Figure  3.8  shows  a  flat  tip  only  for  panels  illustration  purposes.  All  our  wings  have  rounded  tips 
(except  in  the  case  of  the  delta  wings). 
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Figure  3.7  PMARC  Input  File  for  an  Aft  Swept  Wing,  NACA  64 AO  10  Section,  AR=3, 
Taper  Ratio=0.5,  Sweep=45  Degrees 
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Figure  3.7  Cont'd  PMARC  Input  File  for  an  Aft  Swept  Wing,  NACA  64A010  Section, 
AR=3,  Taper  Ratio=0.5,  Sweep=45  Degrees 
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(From  Ref.  [1]) 

The  SECTl  and  SECT2  combination  is  different  for  this  case  of  cambered  wing. 
Instead  of  a  four  digit  NACA  airfoil  section,  we  now  have  a  section  definition  for  which  PMARC 
has  no  automatic  construction  routine.  We  then  have  to  enter  the  coordinates  of  the  airfoil  section 
definition.  For  this  reason,  SECT2  is  not  required  anymore,  as  it  is  only  used  in  the  case  where  a  four 
digit  NACA  airfoil  is  wanted. 

We  will  still  observe  two  SECTl  sections:  one  for  each  end  of  the  patch  (or  wing). 
The  differences  between  the  two  reside  in  the  coordinates  of  the  origins  of  each  section  starting  and 
ending  the  patch.  Also,  the  scaling  is  different,  meaning  that  the  outside  edge  of  the  patch  (wing  tip) 
has  a  different  length  than  the  inside  edge  (wing  root).  An  important  field  here  is  INMODE.  It  is  set 
at  four  in  the  first  SECTl  of  the  cambered  wing  input  file  to  indicate  that  the  patch  coordinates 
follow.  In  the  flat  plate  input  file,  the  first  INMODE  was  set  at  five  requesting  that  PMARC  generate 
an  NACA  four  digit  airfoil  section.  The  second  INMODE  is  set  at  zero,  instructing  PMARC  to  copy 
the  section  definition  of  the  previous  section.  The  same  thing  could  have  been  done  with  the 
cambered  wing  instead  of  redefining  the  section  definition  for  the  tip. 
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The  manual  entry  of  the  section  definition  should  be  done  carefully.  Some  entry  will 
not  cause  a  system  error,  but  will  produce  bad  results.  For  instance,  if  the  patch  coordinates  are 
entered  starting  at  the  origin  (0  0  0),  PMARC  will  try  to  satisfy  the  Kutta  condition  at  the  leading 
edge  instead  of  at  the  trailing  edge.  Therefore,  an  important  lesson  learned  while  exploring  PMARC 
was  to  understand  how  the  patches  are  defined  and  folded  (Ref.  [1]).  For  the  case  of  any  wing,  the 
section  definition  always  starts  at  the  trailing  edge,  going  under  the  wing,  around  the  leading  edge 
and  up  towards  the  trailing  edge  (Figure  3.9). 


Side  3  of  folded  wing  patch  with  15  panels  on 
both  upper  and  lower  surfaces. 


Comer  points  of  panels  1  •  15  define  first  section 
of  dp  patch. 

Comer  points  of  pa'*els  30-16  (note  reverse  order) 
define  last  secdon  of  dp  patcL 


Figure  3.9  Wing  Cross-Section  Definition  Sequence. 
(From  Ref.  [1]) 
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There  is  more  than  one  way  to  use  the  coordinate  system  to  define  a  section.  That 
is  controlled  by  the  field  INMODE  in  SECTl .  Depending  on  its  value,  the  three  columns  of  data  will 
have  different  meanings  (refer  to  page  53  of  Ref  [1]).  In  the  present  case,  the  three  columns 
respectively  represent  the  x,  y,  and  z  coordinates  of  each  point  defining  the  section.  This  section 
definition  of  the  NACA  64A010  was  taken  in  Ref.  [7]. 

As  well  illustrated  in  Ref.  [2],  and  in  the  test  case  of  Ref.  [1],  more  than  one  means 
are  available  to  shape  the  required  dimensions  into  our  wing  sections.  In  the  present  case,  we  elected 
to  non-dimensionlize  the  section  definition  and  then  scale  it  to  the  chord  and  tip  dimensions 
respectively.  Another  way  would  have  been  to  enter  the  global  coordinates  of  the  root  and  tip 
sections. 


d.  Wake  and  Streamlines 

As  mentioned  in  the  previous  section,  we  need  to  leave  the  wake  active  (IDWAK=1) 
as  we  do  not  need  the  streamlines  option. 

4.  PMARC  Output  -  Data  Retrieval 

Once  again,  as  for  the  swept  back  wing  of  the  previous  section,  we  are  interested  in  are  the 
lift  coefficients  (Q)  for  various  angles  of  attack.  In  order  to  do  so,  our  input  file  has  to  be  run 
repetitively,  each  time  modifying  the  angle  of  attack.  This  will  create  a  series  of  output  files, 
DATA6,  from  which  the  Q  will  be  read.  Several  lift  coefficients  are  presented  in  this  output  file: 
assembly,  component,  and  total,  in  term  of  wind  or  body  axes.  For  the  present  case,  we  need  the 
total  wing  lift  coefficient  in  terms  of  the  wind  axis  which  refers  to  the  freestream.  This  value  is 
found  towards  the  very  end  of  the  file  under  the  heading  Wind  Axes  -  Total  coefficients.  The  lift 
coefficient  corresponding  to  our  case  is  labelled  CL.  Note  that  when  streamlines  are  requested,  the 
related  results  are  placed  at  the  end  of  the  file,  following  the  coefficient  results. 
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D.  DELTA  WINGS 


A  good  amount  of  experimental  data  is  available  for  delta  wings,  allowing  a  wider  spectrum 
of  comparison.  Four  different  delta  wings  were  modeled,  each  having  different  aspect  ratio  (AR), 
but  all  with  aNACA  0012  airfoil  section  designation.  In  order  to  obtain  finite  geometric  parameters, 
the  span  was  arbitrarily  set  at  10  in  all  cases.  Figure  3.10  shows  those  configurations. 


1.  Theory,  Background  and  Reference  Data 

Prediction  and  control  of  flow  separation  is  a  major  aerodynamic  consideration  in  wing 
design.  Effort  to  avoid  undesirable  compressibility  effects  leads  to  greater  sweep  angle  and  reduced 
section  thickness.  However,  these  efforts  are  offset  by  the  increasing  difficulty  to  prevent  boundaiy- 
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layer  separation.  Schlichting  and  Truckenbrodt  (Ref.  [8]),  and  also  Ashley  and  Landahl  (Ref  [9]) 
studied  the  influence  of  the  delta  wing  aspect  ratio  on  the  lift  coefficients  and  the  boundary  layer 
separation.  The  purpose  of  our  study  is  not  to  draw  any  conclusion  with  regard  to  flow  separation 
and  AR  influence  using  our  PMARC  computation  but  rather  to  try  to  recreate  some  of  the 
experimental  results  obtained  by  the  above  mentioned  researchers. 

Our  comparative  study  will  present  five  sets  of  data.  We  will  first  compare  the  lift 
coefficients  of  the  delta  wings  of  AR  0.83  and  3.97  obtained  by  PMARC  with  the  available 
experimental  data.  Then  a  plot  of  the  lift  curve  slope  as  a  function  of  the  AR  will  be  built  with  the 
PMARC  results  of  four  different  AR  wings.  Finally,  the  chordwise  and  spanwise  pressure 
distributions  for  the  delta  wing  of  AR=3  will  be  established  using  PMARC  and  compared  to  the 
literature. 


2.  Results  and  Discussion 

The  first  delta  wing  modeled  had  an  aspect  ratio  of  0.83.  The  lift  coefficient  for  various 
angles  of  attack  as  obtained  with  PMARC  is  shown  compared  to  the  experimental  data  on  Figure 
3.11.  Very  good  agreement  exists  at  low  angles  of  attack  but  above  fO  degrees,  PMARC 
underpredicted  the  experimental  data.  This  situation  was  encountered  before  in  this  paper  and  as  it 
was  mentioned  then,  this  disparity  at  higher  angles  of  attack  does  not  invalidate  PMARC.  These 
results  were  expected  as  PMARC  only  models  linear  angle  of  attack  effects  and  therefore  does  not 
account  for  nonlinear  phenomena  encountered  at  high  angles  of  attack.  Therefore,  for  the  lower 
angles  of  attack,  PMARC  recreates  the  experimental  data  very  well. 

The  same  phenomena  occur  with  the  delta  wing  of  aspect  ratio  3.97  where  the  CL  obtained 
by  PMARC  steadily  and  linearly  increases  with  angle  of  attack  (Figure  3.12).  The  difference  here 
is  that  PMARC  overpredicted  the  lift  coefficient  values.  Again,  for  the  linear  portion  of  the 
experimental  data,  PMARC  is  in  close  agreement. 
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Angle  of  attack  (Degrees) 

Figure  3.11  Comparison  of  The  Lift  Coefficient  as  a  Function  of  Angle  of 
Attack.  Delta  Wing,  NACA  0012  Section,  AR=0.83. 


Figure  3.12  Comparison  of  The  Lift  Coefficient  as  a  Function  of  Angle  of 
Attack.  Delta  Wing,  NACA  0012  Section,  AR=3.97. 
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Lift-Curve  Slope,  dCI/d(atpha)  -  /radian 


The  next  figure  (Figure  3.13)  shows  the  lift  curve  slope  as  a  ftinction  of  aspect  ratio.  Again, 
the  PMARC  results  display  good  agreement  with  past  work  (Ref.  [9]).  It  is  shown  there  that  the 
lifting-surface  theory  predicted  very  well  two  sets  of  experimental  data  produced  at  different  time 
by  different  people.  We  superposed  the  PMARC  output  onto  these  combined  results  to  show  the 
very  good  agreement. 

We  have  so  far  observed  that  PMARC  seems  to  compute  fairly  accurately  the  lift 
coefficient.  Another  parameter  should  now  be  introduced  in  order  to  better  assess  this  panel  code. 
As  such,  the  last  two  sets  of  results  presented  in  this  section  will  deal  with  the  pressure  distribution 
for  a  delta  wing  of  aspect  ratio  equal  to  three  and  an  angle  of  attack  of  10  degrees. 


Figure  3.13  Comparison  of  Lift-Curve  Slope  as  a  Function  of  Aspect  Ratio. 
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The  chordwise  and  spanwise  distributions  are  shown  on  Figures  3.14  and  3.15.  Figures 
3.16  and  3.17  taken  from  Ref.  [17]  show  experimental  spanwise  pressure  distributions.  Some 
irregularities  appear  to  be  present  in  the  PMARC  output  of  the  spanwise  distribution  towards  the 
edges  of  the  spanwise  sections  considered.  The  pressure  do  not  seem  to  come  back  down  as  shown 
on  Figure  3.17.  Either  we  do  not  have  enough  points  or  the  panel  layout  can  not  pick  up  the  correct 
behavior  of  the  flow  in  these  regions  of  greater  gradient.  Note  that  the  PMARC  spanwise 
distribution  plotted  are  in  reality  approximations  to  the  real  spanwise  distribution.  What  is  plotted 
is  actually  the  pressure  distribution  along  a  slanted  spanwise  cross-section  that  corresponds  to  panel 
numbers  easy  to  retrieve  out  of  the  output  file.  More  detail  on  this  topic  will  be  given  in  subsection 
4.  Also,  the  reader  should  keep  in  mind  that  Figures  3.16  and  3.17  represent  visous  flows. 


Figure  3.14  Chordwise  Pressure  Distribution  at  10  Degrees  Angle  of  Attack.  Delta 
Wing,  NACA  0012  Section,  AR=3. 
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Figure  3.15  Approximate  Top  Spanwise  Pressure  Distribution  at  10  Degrees  Angle 
of  Attack  for  1/3  and  2/3  Chord  Locations.  Delta  Wing,  NACA  0012  Section,  AR=3. 


Figure  3.16  Spanwise  Pressure  Distribution  Over  a  Flat  and  Cambered  Delta 
Wing,  AR=1.15,  AOA=25  Degrees.  (From  Ref  [17]) 
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3.  PMARC  Input  Instructions 

Five  input  files  were  required  to  obtain  all  the  results  presented  in  this  section.  They  were 
all  identical  except  for  the  section  sizing  and  coordinates  (SECTl)  and  SREF.  These  sets  of  fields 
correspond  to  the  specific  wing  geometiy  as  well  as  to  the  root  and  tip  section  origins.  One  sample 
for  the  delta  wing  of  aspect  ratio  0.83  is  shown  in  Figure  3.18.  The  fields  mentioned  above  are  also 
the  only  differences  between  the  delta  wing  input  files  and  the  swept  back  flat  plate  input  file  of  our 
first  test  case.  Therefore,  no  time  will  be  spent  here  to  go  over  each  section  of  the  input  file.  Instead, 
we  will  discuss  a  particularity  of  the  delta  wing  modelling  in  PMARC.  Also  a  sensitivity  test 
performed  on  some  fields  will  be  discussed. 
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a.  Shaping  a  delta  wing  tip 

In  order  to  bring  our  NACA  0012  airfoil  section  from  its  full  chord  length  at  the  root 
to  a  point  at  the  tip,  an  initial  attempt  was  made  by  scaling  it  to  zero  (SCALE=0);  but  the  system 
crashed.  To  alleviate  this  situation,  it  was  found  that  a  number  other  than  zero  had  to  be  entered.  So 
to  keep  the  integrity  of  oiu'  geometry  a  very  small  number  was  selected,  SCALE=0.05.  In  fact,  that 
was  the  smallest  value  that  would  allow  the  system  to  run  successfully.  This  number  could  vary  as 
it  will  be  seen  in  a  later  case  where  it  was  successfully  set  at  0.01 .  This  variation  is  a  function  of  the 
overall  geometry  of  the  wing,  and  may  be  related  to  the  number  of  panels  on  the  surface. 

Although  it  has  not  been  confirmed  yet,  the  above  situation  may  be  related  to  the 
size  of  the  panel  in  the  region  of  the  tip.  Panels  that  are  too  small  may  cause  problems.  This  area 
needs  more  investigation  as  we  did  run  cases  where  the  output  warned  the  user  about  panels  having 
virtually  no  area,  but  still  produced  valid  lift  coefficients  and  pressure  distributions. 

b.  Testing  Certain  Fields 

The  software  sensitivity  to  the  number  of  panels  on  the  tip  was  brought  up  before, 
but  we  want  to  warn  potential  PMARC  users  that  it  is  even  more  true  in  the  case  of  very  small  tips 
like  for  the  delta  wing.  The  wing  tip  section  of  the  input  file,  the  TNPS  field  of  the  PATCH2  line 
should  be  set  at  a  very  low  value.  The  smaller  the  tip,  the  smaller  the  value.  In  the  case  of  the  delta 
wing,  TNPS=1  seemed  to  be  the  most  efficient  value  for  a  successful  run. 

4.  PMARC  Output  -  Data  Retrieval 

We  already  discussed  where  to  find  the  lift  coefficient  in  the  output  file  DATA6.  We  need 
now  to  indicate  how  to  obtain  the  data  for  the  pressure  distribution  shown  in  the  above  plots.  It  will 
be  different  depending  along  which  direction  (spanwise  or  chordwise)  we  need  the  distribution. 
Either  way,  a  clear  picture  of  how  the  panels  are  labelled  is  crucial.  Figure  3.19  details  this 
important  information  for  an  arbitrary  case  where  TNPS  is  8  and  TNPC  is  12. 
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Figure  3.19  PMARC  Panel  Labelling  Around  a  Delta  Wing. 
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a.  Chordwise  Pressure  Distribution 

In  order  to  plot  the  steady  state  chordwise  pressure  distribution ,  the  user  needs  to 
scroll  down  the  output  file  until  the  time-step  corresponding  to  steady  state  is  reached.  If  the  time- 
step  fields  NTSTPS  and  DTSTEP  are  filled  in  accordance  with  the  instructions  provided  in  Chapter 
n,  the  data  corresponding  to  the  last  time-step  (towards  the  bottom  of  the  file)  are  within  the  steady 
state  regime.  In  the  present  case  we  wanted  to  plot  the  distribution  at  the  root,  and  the  output  labels 
the  root  panels  fi-om  one  to  thirty  under  column  one. 

As  will  be  seen  later,  chordwise  pressure  distribution  at  spanwise  locations  other 
than  the  root  can  also  be  plotted  using  the  PMARC  output.  The  only  additional  step  is  to  locate  the 
column  of  data  with  a  constant  y  value  corresponding  to  the  spanwise  location  desired. 

b.  Spanwise  Pressure  Distribution 

In  order  to  precisely  retrieve  the  Cp  value  for  each  spanwise  station  along  a 
spanwise  cross-section,  one  has  to  select  one  panel  per  column  whose  coordinates  are  the  closest  to 
the  cross-section  and  station  under  consideration.  This  could  be  a  lengthy  process,  but  some 
computer  program  can  be  written  to  do  it  automatically.  What  was  done  with  the  delta  wing 
spanwise  pressure  distribution  was  an  approximation  believed  to  be  good  enough  to  analyze  the 
trend. 


The  fifth  to  last  panel  from  each  of  the  ten  columns  was  selected  and  the  Cp  value 
recorded  for  each.  The  process  was  repeated  with  the  tenth  from  last  panel  of  each  column.  As  the 
panel  selected  in  each  column  are  closer  to  the  trailing  edge,  the  panels  form  a  line  more  and  more 
parallel  with  a  normal  spanwise  distribution.  As  the  panels  get  away  from  the  trailing  edge,  this  line 
moves  more  and  more  towards  the  leading  edge,  ultimately  merging  with  it.  Consequently,  the 
results  shown  for  the  2/3  chord  location  are  more  accurate  than  those  of  the  1/3  chord  location,  but 
the  trend  is  reasonably  the  same. 
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E.  RECTANGULAR  WINGS 


The  rectangular  wing  is  the  simplest  from  the  point  of  view  of  structural  design  and  analysis. 
In  this  section  we  will  put  PMARC  to  the  test  with  rectangular  wings  of  aspect  ratio  1,  4.5,  and  2, 
respectively.  All  wings  will  have  a  NACA  0012  cross-section,  except  for  a  study  shown  at  the  end 
where  the  wing  of  AR=7  is  gradually  thinned  out  to  observe  the  impact  on  the  pressure  distribution. 
Figure  3.20  shows  the  relative  size  of  the  three  wings.  For  the  PMARC  input  file,  we  arbitrarily 
selected  all  chord  lengths  to  be  1,  leading  to  a  span  of  7, 4.5,  and  2  respectively. 
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Figure  3.20  Rectangular  Wings  of  Various  Aspect  Ratio. 


1.  Theory,  Background  and  Reference  Data 

Two  main  references  were  consulted  in  order  to  compare  the  results  obtained  with  PMARC. 
The  first  one  is  a  Vortex-Lattice  software  output  produced  in  1971  by  Margason  and  Lamar  (Ref 
[10]).  The  program  represents  the  lifting  planform  with  a  vortex  lattice  and  solves  the  system  of 
equations  in  a  manner  summarized  in  section  ni-B-l  and  detailed  in  the  above  mentioned  reference. 
The  Vortex  Lattice  Method  is  not  the  most  accurate  method  for  estimating  subsonic  aerod5mamic 
characteristics  of  lifting  planforms.  Consequently,  the  validation  of  the  PMARC  results  should 
account  for  this  fact. 
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The  second  reference  is  a  two-dimensional  computer  panel  code  named  PANEL  developed 
by  Dr.  M.  Platzer  (Ref.  [11]).  For  any  airfoil  coordinates  and  a  selected  number  of  panels,  it 
computes  the  major  aerodynamic  coefficients.  We  will  use  PANEL  pressure  distribution  results  for 
a  NACA  0012  to  compare  with  the  PMARC  output  of  a  large  aspect  ratio  (20)  NACA  0012 
rectangular  wing.  Although  PMARC  is  a  three-dimensional  code,  the  mid-span  results  for  a  high- 
aspect-ratio  wing  can  be  regarded  as  the  same  as  the  two-dimensional  values  as  the  wing  tips  have 
negligible  effect  at  the  root. 

2.  Results  and  Discussion 

a.  Comparing  with  the  Vortex-Lattice  Results 

In  addition  to  comparing  the  lift-curve  slope  of  PMARC  and  Ref  [10],  the  following 
results  will  show  the  influence  of  the  number  of  panels  on  the  final  output.  Figure  3.21  presents  the 
results  obtained  by  Margason  &  Lamar  superposed  by  those  computed  by  PMARC.  For  all  three 


Figure  3.21  Effect  of  Panel  Arrangement  on  Lift-Curve  Slope  for  Rectangular  Wings 
of  Various  Apsect  Ratio. 


44 


wings,  the  lift-curve  slopes  as  obtained  by  both  methods  asymptotically  decrease  as  the  panel  density 
increases.  In  all  cases  PMARC  overpredicted  the  Vortex-Lattice  results  by  five  to  10  percent,  the 
largest  disparities  belonging  to  the  smallest  aspect  ratio.  Also  shown  on  the  figure  is  the  result  when 
the  chordwise  panel  density  is  significantly  higher  than  what  is  presented.  It  is  suspected  that  the 
latter  results  are  the  most  accurate  simply  because  more  panels  allow  for  better  wing  definition.  But 
again,  these  results  do  not  differ  from  Margason  and  Lamar's  highest  panel  density  results  by  more 
than  7.7%. 


b.  Comparing  Chordwise  Cp  Distribution  with  the  2-D  Panel  Code 

In  order  to  efficiently  use  PMARC  to  get  exact  two-dimensional  results,  a  wing  of 
infinite  span  should  be  constructed.  As  discussed  earlier,  a  good  approximation  is  a  wing  of 
sufficiently  large  aspect  ratio.  For  the  case  of  the  rectangular  wing,  we  elected  to  test  a  20  unit  wing 
span  with  an  aspect  ratio  of  20.  The  object  of  the  tests  was  to  vary  the  thickness  of  the  airfoil  section 
from  that  of  a  NACA  0012  to  a  NACA  0001,  and  to  compare  the  pressure  distribution  of  each  wing 
at  an  angle  of  attack  of  ten  degrees,  as  computed  by  PMARC  and  PANEL. 

As  shown  on  Figures  3.22  to  3.28,  very  good  agreement  between  the  two  codes 
exists  for  all  airfoils.  It  is  to  be  noted  that  the  lack  of  smoothness  in  the  PMARC  Cp  curves  is  due 
to  the  small  number  of  panels  used  in  the  computation.  That  also  explains  why  the  peaks  are 
sometimes  cut-off.  Additional  panels  would  have  meant  more  data  points  resulting  in  a  better 
resolution  of  the  curve. 

Figure  3.29  is  a  composite  Cp  plot  for  all  seven  wings  as  computed  by  PMARC. 
This  plot  was  included  for  general  information  to  show  on  a  relative  scale  how  the  pressure 
distribution  changes  as  the  wing  gets  thinner. 
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Location.  Rectangular  Wing,  AR=20,  AOA=10  Degrees,  NACA  0012  Section. 


Figure  3.23  Comparison  of  The  Pressure  Coefficient  as  a  Function  of  Chord 
Location.  Rectangular  Wing,  AR=20,  AOA=10  Degrees,  NACA  0010  Section. 
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Figure  3.24  Comparison  of  The  Pressure  Coefficient  as  a  Function  of  Chord 
Location.  Rectangular  Wing,  AR=20,  AOA=10  Degrees,  NACA  0008  Section. 


Figure  3.25  Comparison  of  The  Pressure  Coefficient  as  a  Function  of  Chord 
Location.  Rectangular  Wing,  AR=20,  AOA=10  Degrees,  NACA  0006  Section. 
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Figure  3.26  Comparison  of  The  Pressure  Coefficient  as  a  Function  of  Chord  Location. 
Rectangular  Wing,  AR=20,  AOA=10  Degrees,  NACA  0004  Section. 


Figure  3.27  Comparison  of  The  Pressure  Coefficient  as  a  Function  of  Chord 
Location.  Rectangular  Wing,  AR=20,  AOA=10  Degrees,  NACA  0002  Section. 
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Figure  3.28  Comparison  of  The  Pressure  Coefficient  as  a  Function  of  Chord  Location. 
Rectangular  Wing,  AR=20,  AOA=10  Degrees,  NACA  0001  Section. 
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Figure  3.29  PMARC  Chordwise  Pressure  Distribution  for  Various  NACA 
Rectangular  Wing  Sections,  AR=20,  AOA=10  Degrees. 
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c.  Spanwise  Cp  Distribution 

Two  wing  thicknesses  were  tested  in  this  section  for  the  rectangular  planform  of 
aspect  ratio  20.  Figure  3.30  shows  the  comparative  pressure  distributions  for  the  NACA  0001  and 
0012,  taken  at  one  third  of  the  chord  length  of  each  wing.  Figures  3.3 1  and  3.32  show  each  wing 
spanwise  Cp  individually  but  taken  at  two  different  chordwise  locations. 


Spanwise  location  -  Unit  length 

Figure  3.30  Spanwise  Top  Pressure  Distribution  at  1/3  Chord  for  Two  NACA  Sections. 
Rectangular  Wing,  AR=20,  AOA=10  Degrees. 


Although  we  do  not  present  experimental  data  to  corroborate  the  PMARC  spanwise 
pressure  distribution  results,  an  intuitive  analysis  leads  to  a  positively  conclusive  validation  of  the 
PMARC  output.  First,  in  all  cases  PMARC  produced  the  expected  elliptical  distribution.  Then,  in 
Figure  3.30,  the  plots  show  higher  Cp's  for  thicker  wings,  which  was  also  to  be  expected.  Figures 
3.31  and  3.32  again  present  the  expected  behavior.  However,  one  area  of  concern  remains  towards 
the  tips  where  a  sudden  peak  indicates  that  the  panel  layout  may  not  pick  up  the  correct  behavior  in 
these  area  of  greater  gradient. 
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Figure  3.31  Spanwise  Top  Pressure  Distribution  at  1/3  and  2/3  Chord  Locations. 
Rectangular  Wing,  AR=20,  AOA=10  Degrees,  NACA  0012  Section. 
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Figure  3.32  Spanwise  Top  Pressure  Distribution  at  1/3  and  2/3  Chord  Locations. 
Rectangular  Wing,  AR=20,  AOA=10  Degrees,  NACA  0001  Section. 
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3. 


PMARC  Input  Instructions 


Except  for  the  fields  specific  to  the  wing  geometry  (CBAR,  SREF,  SSPAN,  SCALE,  SIX, 
and  STY),  the  test  and  panel  prescription  (ALDEG,  TNPS,  and  TMPC)  nothing  is  changed  from  the 
previous  test  case  input  file,  the  delta  wing,  to  the  rectangular  wings  of  various  aspect  ratios.  Figure 
3 .34  shows  the  input  file  for  the  rectangular  wing  of  aspect  ratio  7.  Similarly,  the  input  files  of  the 
other  aspect  ratios  studied  vary  only  by  the  above  fields. 

4.  PMARC  Output  -  Data  Retrieval 

As  mentioned  in  the  previous  section,  the  data  retrieval  required  in  the  PMARC  output  file 
when  interested  in  plotting  pressure  distributions  necessitates  a  complete  understanding  of  the  panel 
numbering  method.  Figure  3 .33  illustrates  this  protocol  for  the  case  of  the  rectangular  wing.  The 
data  retrieval  for  the  lift  coefficient,  the  chordwise  and  spanwise  pressure  distribution  was  discussed 
in  the  previous  section. 


Figure  3.33  PMARC  Panel  Layout  for  a  Rectangular  Wing  and  Its  Tips. 
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F.  F5  WING 

The  F5  wing  is  the  subject  of  an  independent  study.  The  results  of  the  latter  will  be 
combined  and  analyzed  with  those  of  this  paper  in  subsequent  work.  Also,  some  unsteady 
aerodynamic  results  will  be  introduced  later  in  this  paper.  Here  we  merely  present  the  PMARC 
predictions  for  steady  flow.  Figure  3.36  shows  the  half  wing  non-dimensionalized  coordinates. 
From  these  coordinates,  all  the  geometry  was  defined,  an  aspect  ratio  of  2.98  calculated,  and  a 
thickness  ratio  of  0.0480  derived.  Figure  3.35  is  the  PMARC  version  of  the  F5  wing.  The  airfoil 
cross-section  corresponds  to  a  NACA  65004.8  whose  coordinates  were  provided  by  Ref.  [13].  The 
leading  and  trailing  edge  points  were  not  presented  in  that  reference  and  therefore  had  to  be 
introduced  using  a  curve  fit  of  the  available  points.  Figure  3.37  shows  the  given  and  the  fitted  points 
of  the  wing  section. 

The  dimensionless  chord  coordinates  exceed  unity  on  Figure  3.37  to  show  the  relative 
magnitude  of  the  appended  edges.  However  all  data  points  were  divided  by  the  new  chord  length 
in  the  PMARC  section  definition  in  order  to  conserve  a  chord  length  of  unity. 


Figure  3.35  PMARC  Representation  of  The  F5  Wing. 
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Figure  3.37  Leading  and  Trailing  Edge  Fitting.  F5  Wing  Cross-Section. 


1.  Theory,  Background  and  Reference  Data 

a.  Pressure  Distribution  -  Experimental  Results 

The  data  used  for  comparison  with  the  PMARC  results  are  from  experiments 
presented  in  Ref.  [12].  Figure  3.36  displays  the  wing  layout  for  the  experiment  with  the  pressure 
orifices,  in-situ  transducers,  and  accelerometers  (Kulite  and  Endevco).  One  of  the  results  presented 
in  this  study  is  the  pressure  distribution  at  different  spanwise  locations  at  zero  degree  geometric 
angle  of  attack.  The  experimental  results  considered  were  taken  at  a  Mach  speed  of  0.6  and 
accoimted  for  compressibility  effects.  As  will  be  shown,  the  Prandlt-Glauert  compressibility 
correction  factor  will  be  applied  to  the  PMARC  incompressible  results  in  order  to  account  for 
compressibility  effects. 

b.  Velocity  Distribution  -  Area  Rule  Analytical  Results 

In  the  early  days  of  supersonic  flights,  aircraft  designed  to  penetrate  the  sound 
barrier  ran  into  difficulty  as  the  thrust  of  their  engines  could  not  overcome  the  large  peak  drag  near 
Mach  1 .  It  is  after  the  advent  of  the  area  rule  that  dramatic  breakthrough  in  high-speed  flight  started 
to  flourish.  Richard  T.  Whitcomb,  the  father  of  the  area  rule,  reasoned  that  the  variation  of  cross- 
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sectional  area  for  an  airplane  should  be  smooth,  with  no  discontinuities.  That  is  how  the  famous 
Convair  F-102  finally  broke  the  sound  barrier;  after  its  fuselage  was  modified  to  allow  a  reduction 
in  its  cross-sectional  area  in  the  region  where  the  wings  were  located,  in  such  a  way  as  to  smoothen 
the  total  cross-sectional  area  distribution. 

The  analytical  derivations  behind  the  area  rule  is  presented  later  in  Chapter  HI, 
section  G  where  we  verified  the  area  rule  by  comparison  with  Ref  [15].  We  will  see  there  that  the 
flow  can  be  represented  by  a  cross-sectional  flow  (or  cross-flow)  to  which  is  added  a  spatial 
influence  component.  The  velocity  potential  of  the  cross-sectional  flow  and  the  spatial  influence  are 
respectively  expressed  as  follows: 


<P  = 


lit 


Ill[  5(X) 


H—f  ] 
s(xy 


2nJ-s(x)^  ^-z  ^  ^ 


(3.2) 


QJx)[/n2x+ln2(x^-x)]  + 


1 

4Tzh  |x-^| 


(3.3) 


Chapter  lEL,  section  G  presents  various  forms  of  the  velocity  components  (velocity  fields,  or  partial 
derivatives)  of  the  above  equations. 

In  summary,  the  cross-sectional  flow  is  produced  in  each  y-z  plane  (Figure  3.38) 
by  a  distribution  of  two-dimensional  sources  along  a  spanwise  axis.  Each  cross-flow  is  completely 
independent  of  conditions  at  any  other  y-z  plane. 
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Figure  3.38  Low-Aspect-Ratio  Wing  and  its  Cross-Flow  Potential. 

(From  Ref.  [14]) 


Figure  3.39  Substitution  of  Source  Panels.  (From  Ref.  [14]) 
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The  spatial  influence  expression  accounts  for  those  parts  of  the  wing  situated  farther 
upstream  (and  downstream  in  subsonic  flow)  from  the  y-z  plane  under  consideration  (Figure  3.39). 
The  essential  aspect  of  this  theory  is  that  the  spatial  term  is  a  function  only  of  the  total  cross- 
sectional  area  distribution  with  x.  The  shape  of  the  wing  (or  wing-body  combination)  is  irrelevant, 
only  the  cross-section  area  is  important.  This  will  become  a  vital  characteristic  when  complicated 
shapes  are  to  be  modeled  by  computer  software.  An  equivalent  area  body  of  revolution  could  be 
modelled  instead  which  is  the  easiest  geometry  to  build.  Of  course,  the  cross-flow  still  has  to  be 
solved. 


In  the  rest  of  this  chapter,  we  will  partially  apply  the  concept  of  the  two-component 
velocity  potential  on  the  F5  wing  introduced  earlier.  We  will  discover  that  PMARC  does  not 
recreate  valid  potential  distributions.  For  completeness,  the  results  of  the  area  rule  using 
Equations  (3.2)  and  (3.3)  for  the  triangular  wing  will  be  presented  in  the  next  section.  It  will  be 
shown  that  once  more,  PMARC  can  not  reproduce  an  adequate  velocity  potential  distribution. 

2.  Results  and  Discussion 

a.  Pressure  Distribution 

Figure  3.40  presents  the  chordwise  pressure  distribution  on  the  F5  wing  when  at  an 
angle  of  attack  of  zero  degree.  The  spanwise  location  that  produced  this  distribution  is  very  close 
to  the  root.  It  corresponds  to  the  control  points  of  the  first  row  of  panels  whose  edges  define  the 
root.  Although  not  shown,  these  results  agree  with  the  experimental  results. 

The  next  set  of  figures  (Figures  3.41  to  3.44)  better  represents  the  agreement 
between  PMARC  and  the  experimental  results.  For  four  spanwise  locations  of  the  wing,  the  pressure 
distribution  as  calculated  by  PMARC  is  plotted  against  the  experimental  results.  Also,  since  the 
latter  account  for  compressibility  and  the  former  does  not,  the  PMARC  outputs  were  corrected  with 
the  Prandlt-Glauert  factor  and  the  results  superposed  on  the  respective  figures. 
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Figure  3.40  Pressure  Distribution  as  a  Function  of  The  Root  Chord  Location.  F5  Wing, 
NACA  65A004.8  Section,  AR=2.98,  AOA=0  Degree. 


The  PMARC  output  for  the  F5  wing  did  not  allow  for  spanwise  definition  such  as 
to  have  perfect  correspondence  with  the  spanwise  locations  selected  for  the  experiment.  However, 
the  disparity  is  sufficiently  small  to  be  considered  negligible  in  all  but  the  first  case  where  some  care 
should  be  exercised  as  there  is  1 1%  difference  between  the  PMARC  and  the  experimental  spanwise 
location. 


With  this  observation  in  mind,  the  reader  can  still  appreciate  that  all  the  results 
display  correspondent  behavior,  and  for  the  last  three  figures  especially,  the  Prandlt-Glauert 
correction  seems  to  bring  the  PMARC  results  progressively  closer  to  the  experimental  results  as  we 
approach  the  wing  tip. 
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Chordwise  Location  >  Percent  of  Local  Chord 

Figure  3.41  Comparison  of  The  Chordwise  Pressure  Coefficient  at  a  Spanwise  Area.  F5 
Wing,  NACA  65A004.8,  AR=2.98,  AOA=0  Degree. 
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Figure  3.42  Comparison  of  The  Chordwise  Pressure  Coefficient  at  a  Spanwise  Area.  F5 
Wing,  NACA  65A004.8,  AR=2.98,  AOA=0  Degree. 
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Figure  3.43  Comparison  of  The  Chordwise  Pressure  Coefficient  at  a  Spanwise 
Area.  F5  Wing,  NACA  65A004.8,  AR=2.98,  AOA=0  Degree. 


Chordwise  Location  -  Percent  of  Local  Chord 

Figure  3.44  Comparison  of  The  Chordwise  Pressure  Coefficient  at  a  Spanwise  Area.  F5 
Wing,  NACA  65A004.8,  AR=2.98,  AOA=0  Degree. 
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Another  result  presented  is  the  spanwise  pressure  distribution.  Although  we  had  no 
comparative  data  for  these  results,  an  intuitive  analysis  seems  to  validate  them.  Indeed,  the  general 
behavior  combines  those  of  the  rectangular  and  delta  wings.  Spanwise  cuts  of  the  F5  wing  in  the 
region  between  the  forward  apex  and  the  tip  leading  edge  and  similarly  between  the  tip  trailing  edge 
and  the  aft  apex  correspond  to  cuts  on  a  delta  wing  for  which  we  see  a  bucket  shape  Cp  distribution 
(Figure  3.45).  Cuts  in  the  region  between  the  wing  tips  correspond  to  cuts  on  a  rectangular  wing  for 
which  we  observe  an  elliptical  Cp  distribution  (Figure  3.46).  The  PMARC  spanwise  pressure 
distribution  on  the  F5  wing  at  zero  degree  angle  of  attack  combines  these  behaviors  and  is  shown 
on  Figure  3.47. 


Figure  3.45  Typical  Spanwise  Pressure 
Distribution  for  a  Delta  Wing. 


Figure  3.46  Typical  Spanwise  Pressure  Distribution 
for  a  Rectangular  Wing. 
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Figure  3.47  Approximate  Pressure  Distribution  as  a  Function  of  Spanwise  Location.  F5 
Wing,  NACA  65A004.8,  AR=2.98,  AOA=0  Degree. 


b.  Velocity  Potential  Distribution 

Comparison  with  the  analytical  results  of  Ref.  [15]  shows  that  the  PMARC  velocity 
potential  output  is  erroneous.  Figure  3.48  is  the  PMARC  distribution  of  the  x-component 
perturbation  velocity  u  for  various  spanwise  locations.  It  shows  that  the  flow  is  faster  as  the 
spanwise  location  gets  closer  to  the  tip  (that  can  be  seen  by  noting  that  the  peak  velocities  are  higlier 
as  z  increases).  Here,  z  is  the  spanwise  location.  (This  notation  does  not  correspond  to  the  PMARC 
one  but  rather  matches  the  technical  note  from  which  the  area  rule  was  taken  (Ref.  [15])).  The 
perturbation  velocity  u  is  computed  as  shown  in  Equation  (3.4). 


X  °° 
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Figure  3.48  Perturbation  Velocity  on  The  Top  Surface  at  Various  Spanwise  Locations. 
F5  Wing,  NACA  65A004.8  Section,  AR=2.98,  AOA=0  Degree. 


This  behavior  is  obviously  not  realistic.  As  the  wing  thickness  tapers  down  toward 
the  tips,  the  velocity  perturbation  should  decrease  towards  zero.  Figure  3.49  shows  the  distribution 
of  the  same  parameters  but  for  the  lower  side  of  the  wing.  Once  again,  the  velocity  peaks  are 
wrongfully  increasing  toward  the  tip. 

Although  we  do  not  have  any  comparative  data  on  the  F5  wing  spatial  influence 
portion  of  the  total  flow,  it  was  analytically  solved  using  one  form  of  Equation  (3.3)  shown  later  as 
Equation  (3.23),  and  the  results  are  plotted  in  Figure  3.50  for  future  consideration.  The  solved 
analytical  steps  can  be  found  in  Appendix  A. 
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Figure  3.49  Perturbation  Velocity  on  The  Bottom  Surface  at  Various  Spanwise 
Locations.  F5  Wing,  NACA  65A004.8  Section,  AR=2.98,  AOA=0  Degree. 


Spatial  Influence 


Figure  3.50  Spatial  Influence  u-Component  of  The  Velocity  Potential. 
F5  Wing,  NACA  65A004.8  Section,  AR=2.98,  AOA=0  Degree. 
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The  cross-sectional  portion  of  the  total  flow  (Equation  (3.2))  could  not  be  solved 
because  one  of  its  terms,  the  two-dimensional  cross-sectional  segment  area  P(x,z),  could  not  be 
constructed.  Only  unidimensional  expressions  could  be  derived  with  the  numerical  data  available. 

3.  PMARC  Input  Instructions 

At  this  point,  the  reader  should  have  acquired  enough  knowledge  about  PMARC  input  to 
adapt  it  to  various  configurations.  The  F5  wing  does  not  show  any  differences  from  what  we  have 
seen  before,  except  for  the  usual  basic  test  condition  and  geometry  fields  (  Figure  3.51). 

One  aspect  of  consideration  which  has  already  been  mentioned  is  the  importance  of  the 
order  in  which  the  airfoil  section  coordinates  are  defined.  Figure  3.52  shows  a  pressure  distribution 
calculated  by  PMARC  for  the  same  F5  wing,  but  with  a  section  definition  starting  at  the  leading 
edge  instead  of  the  trailing  edge  (Figure  3.53)  as  required  by  the  PMARC  protocol  for  a  wing.  As 
it  can  be  seen,  although  some  of  the  results  are  relatively  close  from  the  valid  distribution,  major 
irregularities  exist,  as  the  Kutta  condition  is  not  met  at  the  trailing  edge.  (Note  that  on  Figure  3.52 
the  chord  coordinates  are  not  non-dimensionalized.) 

4.  PMARC  Output  -  Data  Retrieval 

Again,  all  techniques  for  retrieval  of  the  chordwise  and  spanwise  pressure  distribution  data 
were  introduced  in  a  previous  section.  Section  E-4  of  this  chapter  mentionned  the  importance  of  a 
complete  understanding  of  the  panel  numbering  method  in  order  to  refer  to  the  correct  data. 
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Figure  3.53  Bad  PMARC  Input  File  With  Reverse  Section  Definition.  F5  Wing, 
65A004.8  Section,  AR=2.98.  (Last  Few  Lines  Not  Shown). 
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TRIANGULAR  WING 


The  following  wing  is  also  called  delta  wing.  We  changed  the  terminology  to  triangular 
to  ensure  distinction  from  the  delta  wings  studied  earlier  and  also  to  match  the  literature  which  was 
used  as  reference.  The  triangular  wing  modeled  here  has  a  thickness  distribution  represented  by  the 
Equation  (3.5). 

y  .  KxA  ■  26(x-M)(1.;c)  (3,5) 

a 


which  is  a  parabolic  arc  in  all  X,Y  planes  and  which  has  linearly  decreasing  h  in  all  Y,Z  planes. 
(Note  that  the  roll,  pitch,  and  yaw  axis  correspond  to  the  x,  y,  and  z  axis  respectively  in  PMARC  and 
to  the  X,  z,  and  y  axis  respectively  in  our  reference  literature  (Ref.  [15]).  Figure  3.54  illustrates  the 
wing.  Its  aspect  ratio  is  2/3  and  for  the  purposes  of  analytical  and  computer  analysis,  its  root  chord 
will  be  set  at  one.  Using  these  data,  it  can  be  shown  with  basic  relationships  that  the  semi-span  at 
the  trailing  edge  is  1/6,  and  that  it  varies  in  the  x-direction  like  x/6,  where  x  is  the  non- 
dimensionalized  root  chord  location. 


Figure  3.54  Triangular  Wing,  Parabolic  Arc 
Section,  AR=2/3. 
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Unlike  the  previous  sections,  the  reference  theory  will  be  discussed  in  greater  length  before 
being  compared  to  the  PMARC  results.  The  theory,  introduced  in  the  previous  section,  will  first  be 
presented  in  its  general  form,  then  we  will  analytically  rework  an  example  presented  in  the 
Reference  in  order  to  ascertain  the  validity  and  understanding  of  the  theory.  Once  our  results 
corroborate  those  of  the  author,  we  will  finally  proceed  with  the  PMARC  computation  and 
comparison.  As  it  will  be  seen,  major  disparities  exist  in  the  predicted  velocity  distributions. 

1.  Theory,  Bacl^round  and  Reference  Data 

In  the  early  fifties,  Friedrich  Keune  (Ref.  [15])  presented  a  theory  on  low  aspect  ratio  wings 
with  small  thickness  at  zero  lift  in  subsonic  and  supersonic  flow.  He  showed  that  a  first 
approximation  of  the  flow  around  low  aspect  ratio  wings  can  be  obtained  by  dividing  the  flow  into 
two  components:  the  cross-sectional  flow  {phi)  and  the  spatial  influence  {omega). 


=  {(9(x,y^yQ{x,y^)]U^  (3.6) 


a.  Approximate  Theory 

The  starting  point  of  this  approximate  theory  is  the  equation  representing  the  exact 
velocity  potential  of  a  flow.  With  a  thickness  distribution  over  the  wing  plane  represented  by: 

y  =  ±h{x,z)  (3.7) 

the  source  distribution  in  the  wing  area  (Figure  3.55)  is  depicted  by: 

.  IhJX.OU’^  (3.8) 
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Therefore  the  exact  velocity  potential  is  described  with  the  following  double  integral: 
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(3.9) 


where  F(x,z)  is  the  wing  area. 


Figure  3.55  Two-Dimensional  Flow  Around  a  Low- Aspect-Ratio  Wing  with 
Thickness  but  Zero  Lift.  (From  Ref.  [15]). 


When  considering  low  aspect  ratio  wings,  the  ordinate  zeta  of  the  location  of  any 
source  element  may  be  considered  small  in  comparison  with  the  root  chord  of  the  wing,  therefore 
allowing  a  Taylor  series  development  of  the  exact  equation  for  small  zeta.  Only  the  first  order  terms 
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are  kept,  the  higher  order  terms  being  negligible.  It  is  then  possible  to  develop  the  expression  into 
a  two-dimensional  (Equation  (3.10))  and  a  rotational  symmetric  (Equation  (3.11))  component  of  the 
potential  flow  .  The  derivation  is  detailed  in  Chapter  H  of  Ref.  [15]. 


—  =  —r^^\h  (x,Q  -  O(a^lno) 

[/«•  nJsi(x) 
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(3.11) 


The  two-dimensional  potential  flow  (Equation  (3.10))  is  caused  by  a  two- 
dimensional  source  distribution  2*Ui,rf*(d/dx)[h(x,zeta)]  on  a  line  of  constant  chord  location  and  at 
a  thickness  location  of  zero,  between  two  outer  spanwise  limits  (Figure  3.56).  These  sources  inflate 
the  flat  plate  to  the  thickness  y(z)=+-h(x=const,z),  or  to  the  cross-sectional  area  defining  the  section. 

The  axisymmetric  flow  (Equation  (3.11))  is  induced  by  a  source  distribution  along 
the  x-axis  within  the  wing  (Figure  3.56).  This  flow  is  identical  with  that  around  a  slender  body  of 
revolution. 


Now  the  sum  of  the  two  terms  shown  above  represents  the  first  order  term  of  the 
approximation  of  the  velocity  potential  in  the  incompressible  flow: 
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Figure  3.56  Cross-Sectional  Flow  and  Spatial  Influence.  (From  Ref.  [15]). 


The  potential  function  for  incompressible  flow  shown  above  is  given  for  any 
distance  from  the  wing.  On  the  surface  of  the  wing  (y=0),  the  ordinate  z  has  the  same  order  of 
magnitude  as  the  length  zeta^  so  that  z  must  consequently  be  neglected  with  zeta  in  this 
approximation.  This  leads  to  a  transformation  of  the  potential  function  as  detailed  in  Chapter  IV  of 
Ref.  [15].  The  velocity  field  of  the  cross-sectional  flow  at  y=0  is  then  given  by: 
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The  u  component  in  the  above  equation  (phQ  can  also  be  expressed  the  following  non-dimensional 
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The  spatial  influence  is  represened  by: 


n  1  ^  ^  Nr/  /  or  M  1  r^o  „ 

Q  =  —Q^(x)[ln2xM(x  -x)]  +  —  /  "  - — — - (3.19) 
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The  u  component  in  the  above  equation  (pmega^)  can  also  be  written; 
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One  marked  particularity  of  this  theory  is  that  the  spatial  influence  (axisymmetric 
flow)  is  only  a  function  of  the  cross-sectional  area.  Therefore  the  wing  can  be  represented  by  a 
simpler  geometiy  with  equivalent  sectional  area  (like  its  equivalent  body  of  revolution)  in  order  to 
simplify  the  calculations 

b.  Analytical  Application  Example  Using  The  Triangular  Wing 

Using  a  computer  analytical  solver  (Ref.  [16]),  the  u  component  of  both  terms  of 
the  flow  were  separately  computed  for  the  triangular  wing  described  at  the  beginning  of  this  section. 
For  this  wing,  the  segment  of  the  cross-sectional  area  for  positive  ordinates  z  is: 

P(x^)  .  46o(1-x)(j:-M)J£L  (3.24) 

2a  a 


and  the  cross-sectional  area  is  given  by: 

Q(x)  =  46ax\l-x)  _  (3.25) 


(1)  Spatial  Influence.  The  u  component  of  the  spatial  influence  was 
computed  using  two  equations.  The  first  one  is  the  general  equation  (Equation  (3.23))  and  the  second 
one  is  the  specific  equation  for  the  triangular  wing  (Equation  (3.26))  as  given  by  Ref.  [15]. 


^  =  -i~l-  4  -^(l-2x)  -  —(l-3x)ln4x(l-x)  (3.26) 

6  a  n  l-x  n  n 
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Both  equations  would  be  expected  to  produce  the  same  results  as  Equation  (3.26)  is  directly  derived 
from  Equation  (3.23).  As  Figure  3.57  shows  however,  small  disparities  are  observed,  especially  past 
the  mid-chord  point.  These  differences  are  explained  by  the  accuracy  of  the  integration  methods 
used.  But  more  importantly,  the  figure  also  indicates  that  we  were  able  to  exactly  reproduce  the 
results  of  Ref.  [15]  using  the  general  Equation  (3.23).  Appendix  B  presents  the  Equations  (3.21)  and 
(3.26)  at  different  stages  of  their  solution. 


Spatial  Influence 


Figure  3.57  Comparison  of  Spatial  Influence  u-Component  of  The  Velocity  Potential  as 
a  Function  of  Chord  Location.  Triangular  Wing,  Parabolic  Arc  Spanwise  Section, 
AR=2/3,  AOA=0  Degree. 


(2)  Cross-sectional  flow.  Once  again,  two  equations  were  solved  for 
comparison  purposes.  The  specific  equation  (Equation  (3.27))  as  given  in  Ref.  [15]  and  the  general 
equation  (Equation  (3.17)).  Appendix  C  presents  the  equations  at  different  stages  of  their  solution. 
The  results  were  identical  £ind  in  good  agreement  with  those  obtained  by  Ref.  [15].  Figure  3.58 
presents  these  results. 
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Cross-sectional  Flow 


Non-dimensionalized  Spanwise  Location  (z/sigma) 


Figure  3.58  Comparison  of  The  Chordwise  Cross-Sectional  Flow  u-Component 
Velocity  Potential  as  a  Function  of  Non-Dimensionalized  Spanwise  Location.  Triangular 
Wing,  Parabolic  Arc  Section,  AR=2/3,  AOA=0  Degree. 
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(3)  Total  Flow.  When  the  spatial  influence  and  the  cross-sectional  flow 
results  are  added  up,  we  obtain  the  curves  depicted  in  Figure  3.59,  which  correspond  very  well  to 
those  computed  by  Ref.  [15]. 


Total  Flow 


Figure  3.59  Comparison  of  The  Combined  u-Components  of  The  Chordwise  Basic 
Flow  as  a  Function  of  Non-Dimensionalized  Spanwise  Location.  Triangular  Wing, 
Parabolic  Arc  Spanwise  Section,  AR=2/3,  AOA=0  Degree. 


2.  Results  and  Discussion 

As  for  the  F5  wing,  PMARC  seems  not  capable  of  recreating  a  realistic  x-component 
velocity  distribution  on  the  triangular  wing.  Figure  3.60  once  again  shows  velocity  peaks  that 
increase  as  the  spanwise  location  gets  towards  the  tip.  This  is  a  behavior  that  would  be  more 
representative  of  a  wing  of  constant  thickness.  In  our  case,  the  tip  is  modelled  to  be  of  virtually  zero 
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thickness,  therefore  the  velocity  perturbation  should  approach  zero  as  we  get  closer  to  the  tip.  On 
Figure  3.60,  zbar  is  the  non-dimensional  spanwise  location  divided  by  the  dimensionless  semi-span. 
As  mentioned  before,  z  corresponds  to  they  axis  in  the  PMARC  notation  and  was  left  as  is  to  match 
our  reference  notation  (Ref.  [15]). 


Figure  3.60  Spanwise  Perturbation  Velocity  as  a  Function  of  Chord  Location. 
Triangular  Wing,  Parabolic  Arc  Spanwise  Section,  AR=2/3,  AOA=0  Degree. 


In  order  for  the  reader  to  better  appreciate  what  behavior  is  expected,  we  extracted  a  figure 
from  Ref.  [15]  which  shows  the  perturbation  velocity  distribution  on  the  triangular  wing. 
Interpretation  of  this  figure  corroborates  the  intuitive  expectation  of  progressively  smaller 
perturbation  velocity  from  root  to  tip.  See  Figure  3.61  below. 
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Figure  3.61  Isotaches  in  The  Wing  Plane.  Triangular  Wing,  Parabolic  Arc  Spanwise 
Section,  AR=2/3,  AOA=0  Degree,  z/x  Proportions  Not  to  Scale.  (From  Ref.  [15]) 


To  ensure  that  the  erroneous  results  were  not  due  to  our  geometry  which  in  turn  could  have 
caused  matrix  singularities,  a  series  of  variation  of  our  triangular  wing  were  modeled.  The  tips  were 
truncated,  the  leading  and  trailing  edges  were  rounded,  and  the  panel  density  and  distribution  were 
varied  Larger  aspect  ratio  delta  wings  were  also  modelled  and  all  without  any  positive  effect.  The 
same  erroneous  behavior  was  always  present. 

Interestingly  however,  the  problem  is  not  apparent  in  the  case  of  rectangular  geometries. 
As  shown  in  Appendix  D,  a  thickness  tapered  rectangular  wing  was  modelled  and  the  velocity 
perturbation  decreases  towards  the  tips,  as  it  should. 
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3. 


PMARC  Input  Instructions 


The  input  file  for  the  triangular  wing  is  no  different  than  any  of  the  delta  wings  previously 
studied.  Only  the  airfoil  section,  the  tip  location,  the  geometric  fields,  and  the  scaling  vary.  The 
points  defining  the  airfoil  cross-section  were  separately  obtained  substituting  chordwise  locations 
in  the  triangular  wing  expression  of  thickness  distribution  (Equation  (3.5)).  In  the  following 
paragraphs,  one  of  the  above  differences  will  be  discussed,  the  tip  scaling,  and  an  important 
observation  will  be  brought  to  light  regarding  the  spanwise  cross-section  definition  in  PMARC. 

a.  The  Tip  Scaling 

As  it  can  be  observed  in  the  triangular  wing  input  file  (Figure  3.62),  the  tip  scaling 
was  made  smaller  than  what  was  done  for  the  delta  wings.  It  was  mentioned  then  that  we  wanted  as 
small  a  tip  as  we  could  get  without  creating  too  many  panels  with  virtually  zero  area.  In  the  case 
of  the  delta  wings,  SCALE=0.05  was  an  optimum  value  which  would  generate  a  relatively  low 
number  of  zero  area  panels  .  With  the  triangular  wing,  this  scaling  would  change  the  aspect  ratio 
of  the  wing  as  it  corresponds  to  a  larger  proportion  of  the  root  chord.  We  then  had  to  go  down  as  low 
as  SCALE=0.01  to  conserve  the  triangular  shape,  although  that  increased  the  number  of  zero  area 
panels.  However,  these  panels  account  for  a  small  proportion  of  the  total  number  of  panels,  still 
allowing  for  a  good  flow  solution. 

b.  Spanwise  Cross-Section  Definition 

The  spanwise  cross-section  of  the  wing  as  defined  in  the  input  file  of  Figure  3.62 
has  a  linearly  decreasing  thickness.  The  cross-section  would  then  have  a  diamond  shape.  This  is  due 
to  the  fact  that  we  prescribed  the  coordinates  for  the  origin  of  the  root  chordwise  cross-section  and 
also  for  the  tip  cross-section  and  from  there,  PMARC  seems  to  linearly  interpolate  all  the  origin 
coordinates  of  all  chordwise  cross-sections  in  between.  The  delta  wings  studied  earlier  also  had  this 
diamond  shape  spanwise  cross-section. 

If  the  spanwise  cross-section  has  a  non-linearly  decreasing  thickness  from  root  to 
tip  creating  an  oval  shape  or  double  parabolic  arc  shape,  then  the  PMARC  construction 
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would  be  more  intricate.  Although  not  tried,  it  is  suspected  that  several  different  airfoil  section 
definitions  would  be  required  at  various  spanwise  locations,  each  showing  a  different  thickness  to 
chord  ratio.  Enough  of  these  intermediate  section  definitions  would  be  required  to  create  a  smooth 
spanwise  cross-section. 

4.  PMARC  Output  -  Data  Retrieval 

a.  Direct  Method 

In  order  to  reproduce  the  format  of  the  results  presented  in  our  reference  literature 
(Figure  3.63),  an  intricate  data  retrieval  is  required.  As  can  be  seen  from  that  figure,  we  need  arrays 
of  velocities  taken  at  constant  jc  values,  for  various  increasing  z  values  (y  in  PMARC).  To  do  that, 
we  need  to  find  the  beginning  of  the  steady  state  results  in  DATA6.  Then,  for  each  block  of  results 
(which  corresponds  to  columns,  or  strips  of  panels)  the  velocity  Fx  at  x=0. 1  is  recorded.  The  same 
process  is  repeated  for  x=0.2  and  all  other  chordwise  locations  desired. 

Unless  the  panel  density  is  very  high,  one  won't  find  in  the  output  the  chordwise 
locations  that  exactly  correspond  to  locations  plotted  in  Figure  3.63,  (x=0.1,  0.2,  0.3,  0.4  ...). 
Therefore,  an  interpolation  will  be  required.  In  that  regard,  it  was  observed  that  quadratic 
interpolations  produce  slightly  different  results  than  linear  interpolations.  The  difference  being 
small,  the  latter  method  would  be  acceptable  for  the  type  of  comparative  study  in  question. 

b.  Indirect  Method 

Another  data  retrieval  technique  is  a  lot  simpler  but  some  interpretation  of  the 
resulting  plots  is  needed  to  compare  with  the  reference  plots.  All  that  is  required  here  is  to  plot  the 
X  column  against  the  Vx  column  for  selected  blocks  of  data  (or  panel  strips).  This  produces  plots  as 
presented  in  Figure  3.60.  The  difference  is  that  instead  of  dealing  with  perturbation  velocities  as  a 
function  of  spanwise  stations  for  fixed  chordwise  locations,  we  deal  with  perturbation  velocities  as 
a  function  of  chordwise  stations  for  fixed  spanwise  locations. 
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Figure  3.63  Flow  Results  From  F.  Keune  in  The  Plane  y=0  (PMARC  z).  Triangular 
Wing,  Parabolic  Arc  Spanwise  Section,  AR=2/3,  AOA=0  Degree.  (Fjom  Ref.  [15]). 


c.  Perturbation  Velocity 

The  data  retrieved  in  DATA6  are  x-component  velocities.  Plotting  these  values 
instead  of  those  shown  in  the  reference  plots  (perturbation  velocities)  would  only  make  a  difference 
in  the  abscissa  and  ordinate;  the  behavior  of  the  curves  would  be  the  same.  However,  we  were 
interested  in  the  perturbation  velocities  which  are  simply  the  velocities  found  in  the  DATA6  minus 
the  ffeestream  velocity,  all  divided  by  the  product  of  the  freestream  velocity  and  the  dimensionless 
maximum  thickness. 


u 


U6 

oo 


(3.4) 
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H.  CONCLUSION 


In  this  chapter,  we  demonstrated  that  lift  and  pressure  coefficients  generated  by  PMARC 
accurately  compare  with  the  wide  range  of  experimental  and  theoretical  results  consulted,  or  with 
higher  order  panel  code  results.  However,  disparities  were  found  in  the  area  of  velocity  distribution 
as  the  PMARC  results  showed  increasing  velocity  perturbation  towards  the  tips  of  the  triangular 
wing  whose  thickness  actually  decreases  towards  the  tips.  On  the  other  hand,  it  was  also  shown  that 
the  perturbation  velocity  behavior  was  adequate  in  the  case  of  a  thickness  tapered  rectangular  wing. 
Another  area  that  requires  more  attention  is  the  wing  tip  modelling.  It  would  have  to  be  ascertain 
that  panels  with  virtually  no  area  indeed  do  not  corrupt  any  of  the  output. 

The  cases  studied  in  this  chapter  include  swept  back  wings,  delta  wings,  rectangular  wings 
and  the  F5  wing.  Instructions  were  given  on  how  to  build  each  model  and  make  the  input  file.  Also, 
details  on  data  retrieval  from  the  PMARC  output  file  were  provided  for  each  case  studied. 
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IV.  SOFTWARE  VALIDATION  -  BODIES  OF  REVOLUTION 

A.  INTRODUCTION 

The  previous  chapter  studied  the  PMARC  capability  to  recreate  existing  aerodynamic 
parameters  of  standard  shape  wings  in  steady  state.  This  chapter  has  the  same  objective  for  bodies 
of  revolution.  The  body  of  revolution  is  the  longitudinal  integration  of  circular  cross-sections. 
Computational  advantages  are  associated  with  this  technique  as  was  introduced  earlier.The  body  of 
revolution  can  be  one  of  independent  shape  or  one  which  is  an  equivalent  representation  of  the 
spanwise  cross-sectional  area  of  a  given  wing.  The  two  types  will  be  treated  here.  More  specifically, 
the  equivalent  area  body  of  revolution  for  the  F5  wing  and  the  triangular  wing  studied  previously 
will  be  modelled.  A  spindle  will  also  be  analyzed.  The  spindle  is  a  body  of  revolution  symmetric 
about  all  axes  originating  at  the  mid-chord  point.  In  all  cases,  the  emphasis  will  be  on  the  pressure 
distributions.  It  is  because  there  exist  important  relationships  between  the  aerodynamic  coefficients 
of  a  wing  and  those  of  its  equivalent  body  of  revolution  that  we  will  build  and  analyze  the  latter. 
Note  that  the  demonstration  of  the  above  mentioned  relationships  and  their  utilities  will  be  the 
subject  of  follow-up  work.  For  now,  we  will  only  present  and  discuss  the.  chordwise  pressure 
distribution  at  various  circumferential  locations  for  angles  of  attack  of  zero  and  five  degrees,  which 
will  be  compared  against  the  outcome  of  the  slender  body  theoiy. 

Based  on  the  format  of  the  previous  chapter,  each  body  of  revolution  will  have  a  dedicated 
section.  The  respective  bodies  will  first  be  described,  followed  by  an  overview  of  the  applicable 
theory  serving  as  reference.  Then,  the  results  as  obtained  by  PMARC  will  be  presented  and 
discussed.  Also,  the  specific  PMARC  input  steps  for  modeling  the  subject  body  of  revolution  will 
be  explained.  Finally,  the  output  data  retrieval  techniques  are  presented. 
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B.  F5  WING  -  EQUIVALENT  AREA  BODY  OF  REVOLUTION 


The  F5  wing  was  modelled  earlier  in  this  paper  and  several  pressure  distributions  were 
obtained  for  some  specific  conditions  of  flight.  The  correspondent  equivalent  area  body  of  revolution 
will  have  the  same  main  chord  length  of  0.6439  unit.  Figure  4.1  shows  the  F5  wing  superposed  to 
its  equivalent  body  of  revolution.  The  spanwise  cross-sectional  areas  of  the  wing  needed  to  build  the 
body  of  revolution  were  obtained  using  an  indepth  computer  program  designed  by  the  author  and 
presented  later  in  this  section. 


Figure  4.1  F5  Wing  With  Corresponding  Equivalent  Area 
Body  of  Revolution. 


1.  Theory,  Background  and  Reference  Data 


In  order  to  validate  the  PMARC  results  of  pressure  distribution  for  the  F5  equivalent  body 
of  revolution,  we  will  rely  on  past  experience  for  the  general  expected  behavior,  and  on  the 
approximate  slender  body  theoiy  for  the  magnitude  of  the  pressures.  Our  chordwise  pressure 
distribution  analysis  will  be  performed  at  several  circumferential  locations  to  ensure  consistency  in 
the  expected  trends. 
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a.  Expected  General  Behavior 

As  a  general  behavior,  the  total  lift  on  a  body  of  revolution  should  be  zero  as  the 
upward  lift  on  the  front  portion  of  the  body  is  cancelled  by  the  downward  lift  on  the  rear  portion.  In 
terms  of  pressure  distribution,  it  will  be  observed  that  for  a  change  of  angle  of  attack  from  zero  to 
an  arbitrary  value,  the  pressure  differential  (from  zero  to  a  new  value  of  angle  of  attack)  at  a  specific 
aft  chordwise  location  will  be  equal  in  magnitude  but  in  the  reverse  direction  than  that  of  a  related 
rearward  location.  This  means  that  somewhere  along  the  chord  there  is  a  point  where  this  differential 
is  bound  to  be  zero.  This  point  is  not  at  the  mid-chord,  except  for  a  symmetrical  body  of  revolution 
like  the  spindle,  as  we  will  see  later. 

The  behavior  described  above  would  be  modified  for  the  chordwise  distributions 
that  lie  along  the  0  and  270  degree  circumferential  positions.  Not  only  is  the  magnitude  of  the 
pressure  differential  for  a  zero  to  five  degrees  shift  in  angle  of  attack  expected  to  be  the  same  at  two 
related  forward  and  rearward  positions,  but  also  the  direction  of  the  pressure. 

b.  Slender  Body  Theory 

According  to  the  slender  body  theoiy,  the  pressure  on  a  body  at  angle  of  attack  is 
approximated  by  the  pressure  at  zero  angle  of  attack  plus  the  pressure  differential  between  its  value 
at  zero  and  at  the  new  angle  of  attack.  This  differential  can  be  calculated  using  the  following  slender 
body  theory  derivation  in  which  R  is  the  radius  distribution  and  Q  the  cross-sectional  area 
distribution. 


=  ^Pa.O  "  ^Pda  (41) 


=  2ai?^(x)cos6  (4.2) 
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where: 


m  ■  i(.^) 
% 


(4.3) 


It:  7u 


(4.4) 


Qix)  =  71.3x‘®-352.00x^743.17x*-875.60x’.630.85x^ 


-286.66X  ^81 .70x  ^-14. 13x^1 .42x  2-0.04x 


(4.5) 


Using  the  expected  behavior  description  given  above  and  the  numerical  results  of 
the  slender  body  theory  applied  to  the  F5  wing  equivalent  body  of  revolution,  we  will  now  present 
the  PMARC  results. 

2.  Results  and  Discussion 

Figures  4.3  to  4.8  presented  below  display  the  chordwise  pressure  distribution  at  various 
circumferential  positions  for  two  angles  of  attack,  which  will  allow  comparison  with  intuitive 
expectations  .  Figures  4.9  to  4.12  show  the  pressure  perturbation  when  going  from  zero  angle  of 
attack  to  5  degrees,  as  calculated  by  PMARC  and  the  slender  body  theory.  It  should  be  noted  that 
the  circumferential  positions  of  the  chordwise  body  line  of  analysis  were  chosen  to  correspond  with 
the  control  point  coordinates  of  the  chordwise  strips  of  panels  which  were  the  closest  to  the  points 
separating  the  body  cross-section  in  12  pie  sections  (Figure4.2).  Finally,  although  the  first  and  last 
point  on  all  the  figures  were  plotted,  they  should  be  disregarded.  The  reasons  are  not  well 
understood  but  PMARC  has  difficulty  with  the  pointy  nose  and  tail  when  the  model  is  built  by 
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repetitive  scaling  like  for  the  F5  and  triangular  wing  equivalent  bodies  of  revolution,  by  opposition 
to  individual  cross-section  definition  like  for  the  spindle.  These  two  methods  will  be  reviewed  later. 


Figure  4.2  Circumferential  Angular  Locations  Studied  on  The 
F5  Wing  Equivalent  Body  of  Revolution  Cross-Sections. 
(Dotted  Lines) 


a.  Chordwise  Pressure  Distributions 

The  first  figure  (Figure  4.3)  shows  the  expected  results  that  the  presswe  distribution 
will  be  the  same  for  any  circumferential  position  when  the  body  of  revolution  is  at  zero  angle  of 
attack.  Figure  4.4  presents  the  same  parameters  but  with  the  body  at  five  degrees  angle  of  attack. 
We  now  observe  the  progression  in  pressure  from  top  to  bottom  at  a  particular  aft  section  and 
regression  at  a  particular  rear  section.  This  plot  also  shows  the  relative  chordwise  position  of  the 
minimum  pressure  for  each  chordwise  distribution. 
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Figure  4.3  Chordwise  Pressure  Distribution  for  Various  Circumferential  Angles.  F5 
Wing  Equivalent  Body  of  Revolution.  AOA=0  Degree. 


Figure  4.4  Chordwise  Pressure  Distribution  for  Various  Circumferential 
Angles.  F5  Wing  Equivalent  Body  of  Revolution,  AOA=5  Degrees. 
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Figures  4.5  to  4.8  compare  the  chordwise  pressure  distribution  at  various 
circumferential  positions  for  zero  and  five  degrees  angle  of  attack.  Each  figure  displays  the 
variation  for  two  corresponding  circumferential  positions  (longitudinal  mirror  image)  of  the  body. 
In  all  cases,  the  leeward  pressure  perturbation  is  equal  and  opposite  to  the  windward  pressure 
perturbation,  as  expected.  At  theta  in  the  vicinity  of  90  degrees,  a  pressure  perturbation  exists  but 
remains  approximatively  constant  along  the  chord. 


Figure  4.5  Comparison  of  Chordwise  Pressure  Distribution  Along  Selected 
Circumferential  Angular  Positions  for  Two  Angles  of  Attack.  F5  Wing  Equivalent  Body 
of  Revolution. 
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Figure  4.6  Comparison  of  Chordwise  Pressure  Distribution  Along  Selected 
Circumferential  Angular  Positions  for  Two  Angles  of  Attack.  F5  Wing 
Equivalent  Body  of  Revolution. 


Figure  4.7  Comparison  of  Chordwise  Pressure  Distribution  Along 
Selected  Circumferential  Angular  Positions  for  Two  Angles  of  Attack.  F5 
Wing  Equivalent  Body  of  Revolution. 
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Figure  4.8  Comparison  of  Chordwise  Pressure  Distribution  Along  a  Selected 
Circumferential  Angular  Position  for  Two  Angles  of  Attack.  F5  Wing  Equivalent  Body 
of  Revolution. 


b.  Pressure  Perturbations 

Using  the  PMARC  results  on  Figures  4.5  to  4.8  above,  the  difference  in  pressme 
between  the  body  at  zero  and  five  degrees  angle  of  attack  was  calculated  for  various  circumferential 
positions  and  plotted  against  the  pressure  perturbation  as  computed  with  the  slender  body  theory  at 
the  same  respective  circumferential  positions  and  the  same  change  of  angle  of  attack  (Figures  4.9 
to  4.12). 


For  reasons  given  in  the  previous  subsection,  the  extremities  of  these  plots  should 
be  discarded.  However,  this  restriction  still  leaves  relatively  large  disparities  between  the  two 
perturbation  calculations  for  all  the  cases  considered.  These  disparities  are  larger  towards  the  nose 
and  the  tail  and  progressively  reduce  to  zero  towards  the  shifting  point. 
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The  shifting  point  is  defined  as  the  chordwise  location  where  there  is  no  pressure 
perturbation  following  the  change  of  angle  of  attack.  Ahead  of  the  shifting  point  the  pressure 
perturbations  will  be  in  the  opposite  direction  than  those  after.  It  can  be  observed  that  the  shifting 
points  as  calculated  with  the  slender  body  theory  do  not  always  correspond  to  those  obtained  by 
PMARC.  The  differences  for  the  theta  angles  4.4, 31.5,  58.5,  121.5,  148.5,  and  175.6  degrees  are 
1.5%,  0.0%,  12.7%,  3.1%,  3.3%,  and  6.9%  respectively. 

For  the  theta  position  of  86  degrees,  we  observe  that  the  perturbations  are  close  to 
zero.  A  reading  at  exactly  90  degrees  would  in  fact  reveal  zero  perturbation  when  computed  with 
the  slender  body  theory. 


X/C 

Figure  4.9  Comparison  of  Chordwise  Pressure  Perturbation  at  Selected 
Circumferential  Angular  Positions  for  a  0  to  5  Degree  Change  in  AO  A.  F5 
Wing  Equivalent  Body  of  Revolution. 
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Cp  Variation  Cp  Variation 


x/c 

Figure  4.10  Comparison  of  Chordwise  Pressure  Perturbation  at  Selected 
Circumferential  Angular  Positions  for  a  0  to  5  Degree  Change  in  AO  A.  FS.Wing 
Equivalent  Body  of  Revolution. 
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Figure  4.11  Comparison  of  Chordwise  Pressure  Perturbation  at  Selected 
Circumferential  Angular  Positions  for  a  0  to  5  Degree  Change  in  AO  A.  F5  Wing 
Equivalent  Body  of  Revolution. 


Figure  4.12  Comparison  of  Chordwise  Pressure  Perturbation  at  Selected 
Circumferential  Angular  Positions  for  a  0  to  5  Degree  Change  in  AOA.  F5  Wing 
Equivalent  Body  of  Revolution. 
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3. 


PMARC  Input  Instructions 


The  basic  steps  to  construct  a  body  of  revolution  are  very  similar  from  those  required  for 
a  wing.  The  main  difference  resides  in  the  section  definition  of  the  circles  distributed  along  the 
longitudinal  axis.  There  exists  more  than  one  way  to  construct  a  body  of  revolution.  One  will  be 
presented  here  and  a  second  when  introducing  the  spindle  in  a  later  section.  A  common  requirement 
for  each  method  is  a  set  of  spanwise  cross-sectional  areas  sufficiently  dense  to  properly  define  the 
body  of  revolution.  Section  6  will  present  a  computer  program  designed  by  the  author  that  generates 
41  cross-sectional  areas  and  the  chordwise  locations  of  their  centroids  from  the  geometry  of  the  F5 
wing.  Each  of  the  two  body-of-revolution  construction  methods  presented  have  advantages  and 
disadvantages  which  will  be  discussed  at  the  end  of  each  section.  Figure  4.13  presents  the  input  file 
for  the  first  method  and  the  next  paragraphs  discusses  it. 

a.  Basic  Input  Section 

The  only  fields  which  may  differ  from  the  standard  wings  studied  so  far  are  CBAR, 
SREF,  and  SSPAN.  CBAR  is  taken  to  be  the  length  of  the  longitudinal  axis  of  the  body  or 
revolution.  SREF  is  the  reference  area  upon  which  all  the  coefficients  are  calculated,  and  therefore 
its  selection  could  be  different  depending  on  the  analysis  performed  or  the  results  desired.  In  our 
case,  it  was  the  area  of  the  largest  cross-section.  Finally,  for  the  pressure  and  lift  coefficient  analysis 
we  did  up  to  now,  the  SSPAN  field  did  not  seem  to  be  a  useful  factor.  For  completeness,  we 
selected  it  to  be  half  the  diameter  of  the  largest  cross-section. 

b.  Assembly  and  Component  Coordinate  System  Information 

Nothing  is  different  from  the  standard  wings  studied  so  far. 

c.  Patch  and  Section  Coordinate  System  Information 

This  section  is  significantly  different  from  what  we  have  seen  so  far.  This  section 
also  offers  some  flexibility  which  we  will  demonstrate  by  building  the  F5  wing  body  of  revolution 
one  way,  and  later  the  spindle  using  another  way. 
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Oc{31  199S14:S8 


F5eqArea_DATA5 


Page  3 


fcSECTl 

STX*  0.6006, 

STY*  0.0000 

ALT*  0.0, 

THETA*  0 . 1 

INMODE*  0 , 

TNODS*  0, 

tSECTl 

STX*  0.6028, 

STY*  0.0000 

ALT*  0.0, 

THETA*  0 . 

INMODE*  0 , 

TNODS*  0, 

&SECT1 

STX*  0.6057, 

STY*  0.0000 

ALF*  0.0, 

THETA*  0 . 

INMODE*  0, 

TNODS*  0, 

JrSECTl 

STX*  0.6086, 

STY*  0.0000 

ALF*  0.0, 

THETA*  0 . 

INMODE*  0, 

TNODS-  0, 

fcSECTl 

STX*  0.6122, 

STY*  0.0000 

ALF-  0.0, 

THETA*  0 . 

INMODE*  0, 

TNODS*  0, 

tSECTl 

STX*  0.6160, 

STY*  0.0000 

ALF-  0.0, 

THETA*  0 . 

INMODE*  0, 

TNODS*  0, 

&SECT1 

STX*  0.6201, 

STY*  0.0000 

ALF*  0.0, 

THETA*  0 . 

INMODE*  0, 

TNODS*  0, 

fcSECTl 

STX*  0.6247, 

STY*  0.0000 

ALF*  0.0, 

THETA*  0 . 

INMODE*  0, 

TNODS*  0, 

tSECTl 

STX*  0.6291, 

STY*  0.0000 

ALF-  0.0, 

THETA*  0 . 

INMODE*  0 , 

TNODS*  0, 

tSECTl 

STX*  0.6341, 

STY*  0.0000 

ALF*  0.0, 

THETA*  0 . 

INMODE-  0 , 

TNODS*  0, 

SSECTl 

STX*  0.6389, 

STY*  0.0000 

ALF*  0.0, 

THETA*  0 . 

INMODE*  0 , 

TNODS-  0, 

tSECTl 

STX-  0.6439, 

STY*  0.0000 

ALF-  0.0, 

THETA*  0 . 

INMODE-  0 , 

TNODS*  5, 

&BFNODE 

TNODE-3 , 

TNPC-0,  T 

&WAKE1 

IDWAK*0 , 

IFLXW-l, 

WING  WAKE 

&WAXE2 

KWPACH-2 , 

KWSIDE*4, 

KWPAN2-0, 

NODEW-0 , 

&WAKE2 

KWPACH*!, 

KWSIDE*2  , 

KWPAN2-0, 

NODEW=5 , 

SCALE=  0.0329, 


SCALE*  0.0299, 


SCALE*  0.0276, 


SCALE*  0.0252, 


SCALE*  0.0226, 


SCALE*  0.0195, 


SCALE*  0.0195, 


SCALE*  0.0138, 


SCALE*  0.0113, 


SCALE*  0.0080, 


SCALE*  0.0002, 


SCALE*  0.0000, 


KWLINE=0, 
INITIAL=0, 
KWLINE=0 , 
INITIAL-0, 


fcONSTRM  NONSL  =0,  KPSL  = 


RN  =37210000, 

vise 

=  0.00015'J 

^  NSLBL 

=  1/ 

&END 

NVOLR* 

0, 

NVOLC* 

0, 

&END 

XO* 

-0.1000, 

YO* 

1.5000, 

zo* 

-0.1000, 

INTVSR*  1, 

&END 

XI* 

-0.1000, 

Yl= 

1,5000, 

21* 

0.1000, 

NPTl*  0, 

SEND 

X2  = 

-0.1000, 

Y2  = 

1.5000, 

Z2* 

0.1000, 

NPT2*  20, 

&END 

X3* 

1.1000, 

Y3* 

1.5000, 

Z3* 

-0.1000, 

NPT3*  25, 

&END 

XRO* 

2.0000, 

YRO* 

2.0000, 

ZRO* 

0.0000, 

INTVSC*  1, 

&END 

XRl* 

4.0000, 

YRl* 

2.0000, 

2R1* 

0.0000, 

XR2* 

2.0000, 

YR2* 

2.0000, 

ZR2* 

1.0000, 

S.END 

Rl= 

0.1000, 

R2  = 

1.0000, 

PHIl* 

0.0, 

PHI2*360.0, 

6  END 

NRAD- 

5, 

NPHI* 

12, 

NLEN* 

3, 

&END 

NSTLIN' 

*0, 

&END 

SXO* 

-2.0000, 

SYO* 

1.0000, 

szo* 

-0.5000, 

SU- 

0.0000, 

SD= 

10.0000, 

DS* 

0.2500, 

INTSL*  1, 

&END 

SXO- 

-2.0000, 

SYO* 

1.0000, 

SZO* 

-0.4000, 

su* 

0.0000, 

SD* 

10.0000, 

DS* 

0.2500, 

INTSL*  1, 

&END 

SXO* 

-2.0000, 

SYO* 

1.0000, 

SZO* 

-0.3000, 

su* 

0.0000, 

SD* 

10.0000, 

DS* 

0.2500, 

INTSL*  1, 

&END 

SXO* 

-2.0000, 

SYO* 

1.0000, 

SZO* 

-0.2000, 

su* 

0.0000, 

SD* 

10.0000, 

DS* 

0.2500, 

INTSL*  1, 

&END 

SXO* 

-2.0000, 

SYO* 

1.0000, 

SZO* 

-0.1000, 

su* 

0.0000, 

SD* 

10.0000, 

DS* 

0.2500, 

INTSL*  1, 

SEND 

Figure  4.13  (Cont’d)  PMARC  Input  File  for  The  F5  Wing  Equivalent  Area  Body  of 
Revolution. 
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The  patch  fields  are  not  varying  from  before.  It  is  the  section  fields  that  do.  In  the 
first  method,  a  circle  of  radius  one  is  defined  by  19  points.  As  shown  on  Figure  4.13,  that  circle  is 
then  scaled  to  correspond  to  the  cross-section  radius  at  a  particular  chordwise  location  of  the  body 
of  revolution.  Thus,  we  scale  the  circle  to  zero  at  the  nose  tip  and  progressively  augment  and  then 
reduce  the  scaling  back  to  zero  along  the  longitudinal  axis.  This  process  is  done  by  repeating  the 
SECTl  block  as  many  times  as  we  have  cross-section  definitions. 

The  number  of  panels  on  the  circumference  is  prescribed  in  the  first  break  point 
node  (BPNODE)  following  the  initial  circle  definition.  The  value  attributed  to  INPC  with  this 
method  determines  the  number  of  panels  on  one  half  of  the  circumference,  as  the  model  is  built  by 
joining  two  patches  (Fig  4.14).  The  number  of  panels  longitudinally  from  nose  to  tail  is  prescribed 
in  the  last  SECTl  block  by  the  TNPS  value.  The  panel  distribution  is  selected  as  before  with  the 
TINTC  and  TINTS  fields  next  to  the  respective  panel  fields. 

One  last  note  is  that  the  TNODS  value  is  set  at  0  in  all  the  SECTl  blocks  to  indicate 
that  they  constitute  the  first  or  an  intermediate  section  of  the  patch.  SECTl  of  the  last  block  is  set 
to  5  to  indicate  the  last  section  definition. 

d.  Wake  and  Streamlines 

It  is  important  for  a  body  of  revolution  to  turn  the  wake  option  off.  That  is  done  by 
setting  IDWAK  to  0  in  the  WAKEl  block.  Indeed,  the  body  of  revolution  theoretically  does  not 
generate  a  wake. 
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Figure  4.14  PMARC  Panels  and  Patches  Definition  for  The  Body  of  Revolution. 


4.  Scaled  Circles  Model  Construction  Method  Characteristics 

The  main  advantage  of  this  method  is  its  ease  of  programing.  The  input  file  is  no  longer  than 
for  a  simple  wing.  A  circle  is  defined  once  and  it  is  scaled  at  will.  PMARC  will  create  two  patches 
to  model  the  body  of  revolution  which  are  mirror  images  of  each  other.  One  drawback  when 
compared  to  the  next  method  is  the  apparent  difficulty  with  the  pointed  nose  and  tail.  Even  with  a 
low  number  of  panels  which  allow  a  finite  panel  area  at  the  nose  and  tail,  the  same  erratic  behavior 
is  present  in  the  output  as  the  pressures  at  the  nose  and  tail  break  down.  The  rest  of  the  results 
however  are  good. 
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5. 


PMARC  Output  -  Data  Retrieval 


For  the  purpose  of  our  study,  only  the  pressure  coefficient  distributions  at  various 
circumferential  positions  were  required.  Since  patch  one  and  two  generate  the  same  results  for  the 
control  points  in  mirror  image  positions,  only  the  Cp  from  patch  one  were  extracted.  In  order  to 
obtain  a  Cp  distribution  all  along  the  chord  for  the  theta  position  selected,  one  reading  per  column 
of  data  is  required.  One  or  more  columns  of  data  can  be  skipped  but  care  should  be  exercised  as  to 
conserve  the  proper  smoothness  of  the  distribution.  In  the  DATA6  output  file,  one  column  of  data 
details  the  coordinates  and  aerodynamic  parameters  of  one  circumferential  strip  of  panels  running 
from  bottom  to  top  for  patch  one,  and  from  top  to  bottom  for  patch  two,  as  shown  on  Figure  4.14. 

The  circumferential  position  theta  is  translated  to  a  panel  number  on  the  output  file  by 
looking  through  one  column  of  data  for  the  arctangent  of  the  y  over  -z  ratio  which  will  give  the 
angle  needed.  Once  the  position  of  the  panel  relative  to  the  first  panel  of  that  column  is  determined, 
the  corresponding  panel  in  each  column  of  data  will  be  approximatively  at  the  same  angle  theta.  In 
our  case,  we  wanted  to  study  the  pressure  distribution  at  seven  different  circumferential  positions: 
0,  30, 60, 90, 120,  150,  and  180  degrees.  As  mentioned  earlier,  we  did  not  have  to  study  the  mirror 
image  position  of  these  angles  (0, 330, 300, 270, 240, 210  and  180  degrees)  as  the  results  would  be 
identical.  Because  of  the  automatic  panel  distribution  option  used  in  PMARC,  there  was  virtually 
no  way  to  make  the  control  points  correspond  with  the  theta  desired.  Consequently,  we  had  to  select 
values  close  to  the  ones  above:  4.4, 31.5,  58.5,  86,  121.5,  148.5,  and  175.6  degrees.  But  that  did  not 
matter  for  our  analytical  comparison  as  we  used  the  same  values  in  the  slender  body  theory 
application.  It  should  be  noted  however,  that  by  varying  the  number  of  panels  around  the  body,  one 
would  eventually  find  the  divisions  that  best  suit  the  requirements. 

6.  Computing  The  Cross-Sectional  Area  of  The  Body  of  Revolution 

The  following  subsections  will  review  the  method  employed  to  build  the  F5  wing  equivalent 
body  of  revolution  using  the  F5  wing  PMARC  output,  which  is  our  only  source  of  data.  We  will  first 
go  through  the  steps  indicating  what  data  are  required  from  DATA6,  and  how  this  selection  can 
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become  tremendously  complicated  for  some  types  of  wing,  like  the  F5  wing,  compared  to  simple 
rectangular  wings.  We  will  then  discuss  the  program  that  was  written  to  assist  in  the  construction 
of  the  body  of  revolution,  along  with  some  of  its  output. 

a.  Using  The  PMARC  Output  to  Build  The  Body  Of  Revolution  -  Code 

The  first  objective  while  building  an  equivalent  body  of  revolution  is  to  obtain  the 
section  definition  of  the  many  spanwise  cross-sections  needed.  Once  the  cross-sectional  spanwise 
surface  coordinates  are  identified  for  each  'cut',  all  that  is  left  is  to  calculate  their  respective  area. 
In  the  case  of  the  rectangular  wing,  constructing  spanwise  cross-sections  using  the  PMARC  output 
is  very  easy  and  can  be  done  in  a  matter  of  minutes.  That  is  so  because  all  the  panel  rows  are 
parallel  to  the  leading  edge  and  the  output  provides  the  coordinates  for  each  panel  control  point 
(Figure  4.15).  It  then  becomes  very  simple  to  use  the  coordinates  of  the  top  and  bottom  row  of 
panels  to  define  the  spanwise  cross-section  of  this  wing.  In  cases  where  the  rows  of  panels  are  not 
parallel  to  the  leading  edge,  like  for  the  F5  wing,  it  becomes  more  complicated  as  the  panel  rows  do 
not  match  with  the  spanwise  cross-section  anymore  (Figure  4.15). 


Figure  4.15  PMARC  Panel  Distributions  for  The  Rectangular  and  The  F5  Wings. 


A  series  of  steps  is  then  required  to  create  the  data  points  for  spanwise  cross-sections  fi'om  the  panel 
control  points  provided  by  the  output  file  (DATA6)  of  PMARC.  These  steps  are  summarized  below 
and  each  corresponds  to  a  section  of  the  MATLAB  code  written  to  digitally  perform  them.  The 
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code  is  presented  in  Appendix  E.  It  is  to  be  noted  that  steps  (1)  to  (5)  were  built  to  simplify  the  data 
retrieval  required  for  this  process.  These  first  few  steps  reconstruct  the  chordwise  airfoil  cross- 
section  along  the  span;  an  operation  already  performed  by  PMARC.  The  advantage  of  our  method 
is  that  we  only  need  to  input  the  general  non-dimensional  airfoil  cross-sectional  definition  once 
instead  of  having  to  input  the  several  cross-sectional  definitions  required  along  the  span  in  order 
to  proceed  with  the  determination  of  the  spanwise  cross-sections. 

(1)  Chordwise  Section  Definition.  The  first  step  is  to  define  the  airfoil 
cross-section  in  three  dimensions.  The  x  and  z  coordinates  (in  the  PMARC  axis  definition)  are  in 
this  case  taken  from  the  input  file  (DATA5),  themselves  taken  from  the  literature.  They  coordinates 
(spanwise  locations  of  the  various  two-dimensional  airfoil  sections)  is  taken  in  the  PMARC  output 
file,  DATA6.  They  correspond  to  the  coordinates  of  the  chordwise  panel  control  points  calculated 
by  PMARC  from  the  prescribed  airfoil  section  and  wing  geomehy.  Strictly  on  the  basis  of  the  above 
coordinates,  all  the  following  operations  will  be  performed. 

(2)  Chordwise  Panel  Width.  Because  the  spanwise  panels  are  not  always 
the  same  width,  we  needed  a  routine  that  calculates  the  width  of  each  chordwise  strip  of  panels.  This 
information  is  later  used  along  with  the  sweep  angle  to  determine  the  chordwise  origin  of  the  airfoil 
cross-sections  covering  the  control  points  of  each  strip  of  panels. 

(3)  Origin  Of  The  Chordwise  Panel  Control  Point  Cross-Section. 
Trigonometric  operations  on  the  sweep  angle  and  the  width  and  cumulative  width  of  the  columns 
allow  to  determine  the  origin  in  global  coordinates  of  each  of  the  chordwise  cross-sections  that  cover 
the  control  points  of  each  strip  of  panels. 

(4)  Non-Translated  Coordinates  Of  The  Control  Points  Cross-Sections. 
This  step  constructs  fifteen  coordinate  sets  of  airfoil  chordwise  cross-section.  Each  cross-section 
is  the  representation  of  the  corresponding  panel  strip  control  points.  At  this  stage,  all  the  sections 
constructed  have  their  leading  edge  x  value  at  zero.  However,  all  y  and  z  values  are  representative 
of  the  actual  wing  proportional  geometry. 
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(5)  Translating  The  Cross-Sections.  Using  the  results  of  step  c,  each  cross- 
section  is  now  translated  to  its  actual  relative  position,  completing  the  chordwise  cross-sectional 
definition  of  the  wing.  As  mentioned  earlier,  these  coordinates  are  available  in  DATA6,  but  the  data 
management  required  to  isolate  the  points  needed  is  by  far  more  tedious  than  using  a  single  non- 
dimensional  airfoil  cross-section  definition  along  with  this  program. 

(6)  Fitting  the  Chordwise  Sections  Along  The  Span.  In  order  to  create 
spanwise  cross-sections,  we  needed  to  have  a  means  to  identify  the  coordinates  of  any  surface 
locations  on  a  chordwise  cross-section  instead  of  only  the  coordinates  of  the  panel  control  points. 
To  achieve  this  goal,  we  fitted  the  top  and  bottom  surface  of  each  chordwise  cross-sections  with  high 
order  polynomials.  This  resulted  in  thirty-two  polynomial  expressions  for  one  half  wing. 

(7)  Creating  the  Spanwise  Cross-Sections.  For  different  chordwise 
locations  (forty-three  in  this  case)  and  using  the  polynomials  obtained  above,  we  now  determine 
different  z  coordinates  for  selected  y  locations  (fifteen).  This  will  create  a  series  of  data  points 
defining  forty-three  spanwise  cross-sections.  Each  of  these  spanwise  cross-sections  then  was  fitted 
with  a  polynomial  in  order  to  smooth  out  some  small  irregularities  introduced  by  the  previous  curve 
fitting. 


(8)  Plotting  the  Results.  This  step  helps  visualize  the  spanwise  cross- 
sections  from  leading  edge  apex  to  trailing  edge  apex.  Some  of  these  plots  are  shown  below.  Figure 
4.16  shows  eight  of  the  forty-one  spanwise  cross-sections  for  one  half  wing.  Figure  4. 17  presents 
seven  cross-sections  that  end  at  the  wing  tip  parallel  to  the  chord.  The  flat  look  of  that  tip  reflects 
reality  as  this  geometry  allows  for  the  mounting  of  a  tip  tank  or  missile.  Figure  4.18  illustrates  a 
selection  of  cross-sections  located  in  the  back  portion  of  the  wing  past  the  tips.  The  curvature  of 
some  of  these  cross-sections  is  due  to  the  very  high  order  polynomial  used  in  the  program  which 
were  required  to  get  convergence  at  the  extremity  between  the  top  and  bottom  curves.  These 
tortuous  curves  could  in  fact  be  replaced  by  straight  lines. 
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Unit  Vertical  Location  (z  direction) 


0.02 


FWD  Apex  to  Tip  L.E 


Figure  4.16  Selected  F5  Wing  Spanwise  Cross-Sections  From  Apex  to  Tip  Leading 
Edge.  (Half  Wing). 
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0.01  s 


Tip  L.E.  to  Tip  T.E. 


Figure  4.17  Selected  F5  Wing  Spanwise  Cross-Sections  From  Tip  Leading  Edge  to  Tip 
Trailing  Edge.  (Half  Wing). 


Tip  T.E.  to  Reer  Apex 


Figure  4.18  Selected  F5  Wing  Spanwise  Cross-Sections  From  Tip  Trailing  Edge  to 
Rear  Apex.  (Half  Wing). 


Ill 


(9)  Calculating  the  Area  of  Each  Spanwise  Cross-Section.  The  final  step 
before  constructing  the  equivalent  body  of  revolution  is  to  obtain  the  area  of  each  spanwise  cross- 
section.  Two  methods  were  used  in  order  to  ascertain  the  results.  The  first  one  calculated  the  area 
between  the  curves  representing  the  spanwise  cross-sections  (top  and  bottom  surface  respectively) 
and  the  abscissa  by  integrating  the  polynomials  describing  each  curve.  A  translation  of  these  curves 
was  necessary  before  the  integration  was  performed  in  order  to  eliminate  any  intersection  with  the 
abscissa  which  would  have  rendered  erroneous  areas.  The  second  method  used  a  trapezoid  approach 
to  directly  get  the  area  between  the  two  curves  defining  each  spanwise  cross-sections.  Both  methods 
gave  the  same  results,  as  expected.  Appendix  F  summarizes  these  results. 

(10)  Constructing  the  Fxjuivalent  Body  of  Revolution.  Now  that  all 
spanwise  cross-sectional  areas  are  known,  the  radius  of  the  related  equivalent  area  circles  can  be 
determined  and  located  at  the  corresponding  chordwise  locations,  thus  producing  the  wing 
equivalent  body  of  revolution  (Figure  4.19). 


Figure  4.19  Wing  Spanwise  Cross-Section  and  Its  Equivalent  Area  Circle. 
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(11)  Cross-Sectional  Area  Distribution  It  is  possible  now  to  establish  an 
expression  describing  the  distribution  of  the  spanwise  cross-sectional  area.  All  we  need  to  do  is  to 
find  the  best  fit  for  the  arrays  of  data  points  defining  the  area  as  a  function  of  chordwise  locations. 
In  the  case  of  the  F5  wing,  this  expression  is  the  tenth  order  polynomial  shown  in  Equation  (4.5). 

b.  Note  on  the  Equivalent  Body  of  Revolution  Computer  Code 

In  the  conceptual  phase  of  the  program  described  above,  several  data  management 
complications  arose  that  needed  additional  sub-routines  that  kept  adding  to  its  length  and 
complexity.  However,  as  this  program  was  initially  intended  for  a  one-time  use  tool,  no  effort  was 
expended  in  trying  to  keep  it  computationally  efficient  as  the  subroutines  were  getting  added.  The 
end  result  is  a  program  that  requires  a  relatively  long  running  time,  and  experiences  difficulty  in 
plotting  certain  three-dimensional  representations. 

This  program  is  now  thought  to  be  an  excellent  tool  for  the  application  at  hand. 
Providing  some  spanwise  locations  are  prescribed  as  well  as  the  sweep  angle  and  the  chordwise 
cross-sections  scaling ,  this  program  will  produce  forty-three  spanwise  cross-sections.  If  it  is  to  be 
used  extensively  in  the  future,  it  could  be  modified  to  operate  more  efficiently,  therefore  becoming 
quicker.  The  modification  required  is  to  combine  all  arrays  of  data  produced  by  the  code  into  a  large 
matrix.  All  the  variables  would  then  hold  double  indices  instead  of  one.  This  program  could  also 
work  for  any  other  types  of  wing  for  which  the  general  non-dimensional  chordwise  cross-sectional 
coordinates  are  known.  However,  it  would  first  have  to  be  carefully  reviewed  in  order  to  customize 
it  to  the  wing  under  study. 

c.  Camber  of  The  Body  of  Revolution 

One  aspect  still  needing  some  investigation  is  the  question  of  the  camber.  If  the 
wing  camber  needs  to  be  reflected  in  the  body  of  revolution,  one  possibility  is  to  place  the  centroid 
of  each  spanwise  cross-section  equivalent  area  circles  along  the  mean  camber  line  of  the  wing. 
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c. 


TRIANGULAR  WING  -  EQUIVALENT  AREA  BODY  OF  REVOLUTION 


The  body  of  revolution  making  the  subject  of  this  section  is  made  out  of  equivalent  spanwise 
cross-sections  of  the  triangular  wing  studied  in  Chapter  Ed,  section  G.  The  longitudinal  axis  of  the 
body  of  revolution  will  then  correspond  to  the  root  chord  of  the  wing  which  was  made  unity  through 
non-dimensionalization.  The  spanwise  cross-sectional  area  of  the  wing  was  a  lot  easier  than  for  the 
F5  wing  as  an  analytical  expression  for  the  cross-sectional  area  was  given  in  Ref.  [15]  and 
reproduced  in  Equation  (4.6)  below. 


Q(x)  =  46ax^(l-x) 


(4.6) 


Sigma  is  the  dimensionless  semi-span  and  is  equal  to  1/6.  Delta  is  the  dimensionless  maximum 
thickness  and  will  be  set  at  0.098.  Figure  4.20  shows  the  PMARC  representation  of  the  body  of 
revolution  with  its  parent  wing. 


Figure  4.20  PMARC  Top  and  Side  Representation  of  The  Triangular  Wing  Equivalent 
Area  Body  of  Revolution  and  its  Parent  Wing.  Maximum  Thickness.  Ratio  =  0.09. 

(Two  Figures  not  at  Same  Scale). 
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1. 


Theory,  Background  and  Reference  Data 


In  order  to  validate  the  PMARC  results  of  pressure  distribution  for  the  triangular  wing 
equivalent  body  of  revolution,  we  will  once  again  rely  on  past  experience  for  the  general  expected 
behavior,  and  on  the  approximate  slender  body  theory  for  the  magnitude  of  the  pressures.  Like  for 
the  F5  equivalent  body  of  revolution,  our  chordwise  pressure  distribution  analysis  will  be  performed 
at  several  circumferential  locations  to  ensure  consistency  in  the  expected  trends. 

a.  Expected  General  Behavior 

As  for  the  F5  equivalent  body  of  revolution,  the  total  lift  on  the  triangular  wing  body 
of  revolution  should  be  zero  as  the  upward  lift  on  the  fi-ont  portion  of  the  body  is  cancelled  by  the 
downward  lift  on  the  rear  portion.  In  terms  of  pressure  distribution,  it  will  once  again  be  observed 
that  for  a  change  of  angle  of  attack  from  zero  to  5  degrees,  the  pressure  differential  at  a  specific  aft 
chordwise  location  will  be  equal  in  magnitude  but  in  the  reverse  direction  than  that  of  a  related 
rearward  location.  This  means  that  somewhere  along  the  chord  there  is  a  point  where  this  differential 
is  zero.  This  point  is  not  at  the  mid-chord,  except  for  a  symmetrical  body  of  revolution  like  the 
spindle,  as  we  will  see  later. 

As  for  the  F5  equivalent  body  of  revolution,  the  behavior  described  above  would 
be  modified  for  the  chordwise  distributions  that  lie  along  the  90  and  270  degree  circumferential 
positions.  Not  only  is  the  magnitude  of  the  pressure  differential  for  a  zero  to  five  degrees  shift  in 
angle  of  attack  expected  to  be  the  same  at  two  related  forward  and  rearward  positions,  but  also  the 
direction  of  the  pressure. 

b.  Slender  Body  Theory 

As  already  mentioned  in  the  previous  section,  the  slender  body  theory  dictates  that 
the  pressure  on  a  body  at  angle  of  attack  is  approximated  by  the  pressure  at  zero  angle  of  attack  plus 
the  pressure  differential  between  its  value  at  zero  and  at  the  new  angle  of  attack.  This  differential 
can  be  calculated  using  the  following  slender  body  theory  derivation  in  which  R  is  the  radius 
distribution  and  Q  the  cross-sectional  area  distribution. 
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(4.1) 


O’  -  *  Cp^ 


=  2a/?'(x)cos0  (4.2) 

R(x)  =  V(^)  (4.3) 

71 

«'W  -  (4.4) 

ZTl  71 

0(x)  =  46ax2(l-x)  (4.6) 


Using  the  expected  behavior  description  given  above  and  the  numerical  results  of 
the  slender  body  theory  applied  to  the  triangular  wing  equivalent  body  of  revolution,  we  will  now 
present  the  PMARC  results. 

2.  Results  and  Discussion 

Figures  4.22  to  4.27  presented  below  display  the  chordwise  pressure  distribution  at  various 
circumferential  positions  for  two  angles  of  attack,  which  will  allow  comparison  with  intuitive 
expectations  .  Figures  4.28  to  4.3 1  show  the  pressure  perturbation  when  going  from  zero  angle  of 
attack  to  a  higher  value,  as  calculated  by  PMARC  and  the  slender  body  theory.  It  should  be  noted 
that  the  circumferential  positions  of  the  chordwise  body  line  of  analysis  were  chosen  to  correspond 
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with  the  control  point  coordinates  of  the  chordwise  strips  of  panels  which  were  the  closest  to  the 
points  separating  the  body  cross-section  in  12  pie  sections  (Figure4.21).  Finally,  although  the  first 
point  on  all  the  figures  was  plotted,  it  should  be  disregarded.  As  mentioned  earlier,  the  reasons  are 
not  well  understood  why  PMARC  has  difficulty  with  the  pointed  nose  and  tail  when  the  model  is 
built  by  repetitive  scaling  like  for  the  F5  and  the  triangular  wing  equivalent  bodies  of  revolution, 
as  opposed  to  individual  cross-section  definition  like  for  the  spindle. 


Figure  4.21  Circumferential  Angular  Locations  Studied  on  The 
Triangular  Wing  Equivalent  Body  of  Revolution  Cross-Sections  (Dotted 
Lines). 


a.  Chordwise  Pressure  Distributions 

The  first  figure  (Figure  4.22)  shows  the  expected  results  that  the  pressure 
distribution  will  be  the  same  for  any  circumferential  position  when  the  body  of  revolution  is  at  zero 
angle  of  attack.  Figure  4.23  presents  the  same  parameters  but  with  the  body  at  five  degrees  angle 
of  attack.  We  now  observe  the  progression  in  pressure  from  top  to  bottom  at  a  particular  aft  section 
and  regression  at  a  particular  rear  section.  This  plot  also  shows  the  relative  chordwise  position  of 
the  minimum  pressure  for  each  chordwise  distribution. 
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Figure  4.22  Chordwise  Pressure  Distribution  for  Various  Circumferential  Angles. 
Triangular  Wing  Equivalent  Body  of  Revolution,  Delta=0.098,  AOA=0  Degree. 
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Figure  4.23  Chordwise  Pressure  Distribution  for  Various  Circumferential  Angles. 
Triangular  Wing  Equivalent  Body  of  Revolution,  Delta=0.098,  AOA=5  Degrees. 
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Figures  4.24  to  4.27  compare  the  chordwise  pressure  distribution  at  various 
circumferential  positions  for  zero  and  five  degrees  angle  of  attack.  Each  figure  displays  the 
variation  for  two  corresponding  circumferential  positions  (longitudinal  mirror  image)  of  the  body. 
In  all  cases,  the  leeward  pressure  perturbation  is  equal  and  opposite  to  the  windward  pressure 
perturbation,  as  expected.  At  theta  in  the  vicinity  of  90  degrees,  a  pressure  perturbation  exists  but 
remains  approximately  constant  along  the  chord. 


Figure  4.24  Comparison  of  Chordwise  Pressure  Distribution  Along  Selected 
Circumferential  Angular  Positions  for  Two  Angles  of  Attack.  Triangular  Wing 
Equivalent  Body  of  Revolution,  Delta=0.098. 
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Figure  4.25  Comparison  of  Chordwise  Pressure  Distribution  Along  Selected 
Circumferential  Angular  Positions  for  Two  Angles  of  Attack.  Triangular  Wing 
Equivalent  Body  of  Revolution,  Delta=0.098. 
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Figure  4.26  Comparison  of  Chordwise  Pressure  Distribution  Along 
Selected  Circumferential  Angular  Positions  for  Two  Angles  of  Attack. 
Triangular  Wing  Equivalent  Body  of  Revolution,  Delta=0.098. 
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Figure  4.27  Comparison  of  Chordwise  Pressure  Distribution  Along  a  Selected 
Circumferential  Angular  Position  for  Two  Angles  of  Attack.  Triangular  Wing  Equivalent 
Body  of  Revolution,  Delta=0.098. 


b.  Pressure  Perturbations 

Using  the  PMARC  results  on  Figures  4.24  to  4.27  above,  the  difference  in  pressure 
between  the  body  at  zero  and  five  degrees  angle  of  attack  was  calculated  for  various  circumferential 
positions  and  plotted  against  the  pressure  perturbation  as  computed  with  the  slender  body  theory  at 
the  same  respective  circumferential  positions  and  the  same  change  of  angle  of  attack  (Figures  4.28 
to  4.31). 


For  reasons  given  in  the  previous  subsection,  the  nose  results  of  these  plots  should 
be  discarded.  However,  this  restriction  still  leaves  relatively  large  disparities  between  the  two 
perturbation  calculations  for  all  the  cases  considered.  These  disparities  are  larger  towards  the  nose 
and  the  tail  and  progressively  reduce  to  zero  towards  the  shifting  point. 
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As  described  before,  the  shifting  point  is  defined  as  the  chordwise  location  where 
there  is  no  pressure  perturbation  following  the  change  of  angle  of  attack.  Ahead  of  the  shifting  point 
the  pressure  perturbations  will  be  in  the  opposite  direction  than  those  after.  It  can  be  observed  that 
the  shifting  points  as  calculated  with  the  slender  body  theory  do  not  always  correspond  to  those 
obtained  by  PMARC.  The  differences  for  the  theta  angles  3.4,  32.6, 62.6, 1 17.4, 147.4,  and  176.6 
degrees  are  8.5%,  0.0%,  225.0%,  17.0%,  0.0%,  and  13.8%  respectively.  Unlike  the  F5  wing 
equivalent  body  of  revolution  where  all  the  differences  were  within  acceptable  range,  we  now 
observe  complete  disagreement  for  the  theta  angle  of  62.6  degrees.  As  before,  the  30  degree  range 
shows  perfect  agreement.  For  the  theta  position  of  87  degrees,  we  observe  that  the  perturbations  are 
close  to  zero.  A  reading  at  exactly  90  degrees  would  in  fact  reveal  zero  perturbation  when  computed 
with  the  slender  body  theory. 


Figure  4.28  Comparison  of  Chordwise  Pressure  Perturbation  at  Selected 
Circumferential  Angular  Positions  for  a  0  to  5  Degree  Change  in  AOA.  Triangular  Wing 
Equivalent  Body  of  Revolution,  Delta=0.098. 
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Figure  4.29  Comparison  of  Chordwise  Pressure  Perturbation  at  Selected 
Circumferential  Angular  Positions  for  a  0  to  5  Degree  Change  in  AO  A. 
Triangular  Wing  Equivalent  Body  of  Revolution,  Delta=0.098. 


Figure  4.30  Comparison  of  Chordwise  Pressure  Perturbation  at  Selected 
Circumferential  Angular  Positions  for  a  0  to  5  Degree  Change  in  AO  A.  Triangular 
Wing  Equivalent  Body  of  Revolution,  Delta=0.098. 
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Figure  4.31  Comparison  of  Chordwise  Pressure  Perturbation  at  a  Selected 
Circumferential  Angular  Position  for  a  0  to  5  Degree  Change  in  AOA.  Triangular  Wing 
Equivalent  Body  of  Revolution,  Delta=0.098 

3.  PMARC  Input  Instructions 

The  steps  used  to  construct  the  PMARC  model  of  the  triangular  wing  equivalent  body  of 
revolution  are  the  same  as  for  the  F5  wing  equivalent  body  of  revolution.  The  reader  should  refer 
to  the  previous  section  for  the  details.  The  only  difference  is  in  the  acquisition  of  the  spanwise  cross- 
sectional  areas.  Unlike  for  the  F5  wing,  we  did  not  need  the  assistance  of  a  computer  program  to 
get  them  as  an  analytical  expression  was  available  (Equation  (4.6)).  Appendix  G  shows  the  results 
of  the  analytical  expression  at  various  chord  location  for  use  in  the  slender  body  theoiy.  Figure  4.32 
presents  the  input  file  for  the  triangular  wing  equivalent  body  of  revolution. 

4.  Scaled  Circles  Model  Construction  Method  Characteristics 

The  same  advantages  and  disadvantages  of  the  method  used  for  the  construction  of  the  F5 
wing  equivalent  body  of  revolution  discussed  in  the  previous  section  stand  for  the  triangular  wing 
equivalent  body  of  revolution. 
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The  data  retrieval  techniques  used  here  to  obtain  the  pressure  coefficient  distribution  at 
different  circumferential  theta  angles  are  the  same  as  for  the  F5  wing  body  of  revolution  studied  in 
the  previous  section. 
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D.  THE  SPINDLE 


The  body  of  revolution  making  the  subject  of  this  section  is  symmetrical  about  the 
longitudinal  axis  and  its  mid-chord.  The  length  of  the  longitudinal  axis  was  made  unity  through  non- 
dimensionalization.  The  thickness  distribution  represented  by  the  radius  distribution  is  shown  in 
Equation  (4.7)  below. 


R(x)  =  2x(x-x^) 


(4.7) 


where  tau  is  the  maximum  thickness.  Figure  4.33  shows  the  PMARC  representation  of  the  spindle. 


Figure  4.33  PMARC  Representation  of  The  Spindle. 


1.  Theory,  Background  and  Reference  Data 

In  order  to  validate  the  PMARC  results  of  pressure  distribution  for  the  spindle,  we  will  once 
again  rely  on  past  experience  for  the  general  expected  behavior,  and  on  the  approximate  slender 
body  theory  for  the  magnitude  of  the  pressures.  Like  for  the  F5  and  the  triangular  wing  equivalent 
body  of  revolutions,  our  chordwise  pressure  distribution  analysis  will  be  performed  at  several 
circumferential  locations  to  ensure  consistency  in  the  expected  trends. 

a.  Expected  General  Behavior 

As  for  the  F5  and  the  triangular  wing  equivalent  body  of  revolutions,  the  total  lift 
on  the  spindle  should  be  zero  as  the  upward  lift  on  the  front  portion  of  the  body  is  cancelled  by  the 
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downward  lift  on  the  rear  portion.  In  term  of  pressure  distribution,  it  will  once  again  be  observed  that 
for  a  change  of  angle  of  attack  from  zero  to  5  degrees,  the  pressure  differential  at  a  specific  aft 
chordwise  location  will  be  equal  in  magnitude  but  in  the  reverse  direction  than  that  of  the 
corresponding  rearward  location  (mirror  image  fi-om  mid-chord).  This  means  that  this  differential 
will  be  zero  exactly  at  the  mid-chord  point,  unlike  for  the  F5  and  triangular  wing  equivalent  body 
of  revolutions  where  it  was  further  back. 

The  behavior  described  above  would  be  modified  for  the  chordwise  distributions 
that  lie  along  the  90  and  270  degree  circumferential  positions.  Not  only  is  the  magnitude  of  the 
pressure  differential  for  a  zero  to  five  degrees  shift  in  angle  of  attack  expected  to  be  the  same  at  two 
related  forward  and  rearward  positions,  but  also  the  direction  of  the  pressure. 

b.  Slender  Body  Theory 

As  already  mentioned  in  the  previous  sections,  the  slender  body  theory  dictates  that 
the  pressure  on  a  body  at  angle  of  attack  is  approximated  by  the  pressure  at  zero  angle  of  attack  plus 
the  pressure  differential  between  its  value  at  zero  and  at  the  new  angle  of  attack.  Recall  that  this 
differential  can  be  calculated  using  the  following  slender  body  theory  derivation  in  which  R  is  the 
radius  distribution  and  tau  the  maximum  thickness. 


Cp  =  Cp^.o  ^  ^Pda  (41) 

=  2aJ?I'(jc)cos0  (4.2) 

R(x)  =  2x(x-x^)  (4.8) 
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Jl'(x)  =  2t(1-2x) 


(4.9) 


Using  the  expected  behavior  description  given  above  and  the  numerical  results  of 
the  slender  body  theory  applied  to  the  spindle,  we  will  now  present  the  PMARC  results. 

2.  Results  and  Discussion 

Figures  4.35  to  4.40  presented  below  display  the  chordwise  pressure  distribution  at  various 
circumferential  positions  for  two  angles  of  attack,  which  will  allow  comparison  with  intuitive 
expectations  .  Figures  4.41  to  4.44  show  the  pressure  perturbation  when  going  from  zero  angle  of 
attack  to  a  higher  value,  as  calculated  by  PMARC  and  the  slender  body  theory. 

It  should  be  noted  that  the  circumferential  positions  of  the  chordwise  body  line  of  analysis 
were  chosen  to  exactly  correspond  with  the  control  point  coordinates  of  the  chordwise  strips  of 
panels  which  separate  the  body  cross-section  in  12  pie  sections  (Figure4.34). 

This  exact  selection  was  not  possible  in  the  two  previous  cases.  It  was  made  possible  here 
due  to  the  method  employed  to  build  the  model  as  will  be  discussed  in  the  next  subsection.  Also  as 
a  result  of  that  particular  method,  the  end  points  (nose  and  tail)  were  not  corrupted  and  generated 
valid  results. 


Figure  4.34  Circumferential  Angular  Locations 
Studied  on  The  Spindle  Cross-Sections. 


a.  Chordwise  Pressure  Distributions 

The  first  figure  (Figure  4.35)  shows  the  expected  results  that  the  pressure 
distribution  will  be  the  same  for  any  circumferential  position  when  the  body  of  revolution  is  at  zero 
angle  of  attack.  Figure  4.36  presents  the  same  parameters  but  with  the  body  at  five  degrees  angle 
of  attack.  We  now  observe  the  progression  in  pressure  from  top  to  bottom  at  a  particular  section  aft 
of  the  shifting  point  and  regression  at  a  particular  section  rearward  of  the  shifting  point.  This  plot 
also  shows  the  relative  chordwise  position  of  the  minimum  pressure  for  each  chordwise  distribution. 

Figures  4.37  to  4.40  compare  the  chordwise  pressure  distribution  at  various 
circumferential  positions  for  zero  and  five  degrees  angle  of  attack.  Each  figure  displays  the 
variation  for  two  corresponding  circumferential  positions  (longitudinal  mirror  image)  of  the  body. 
In  all  cases,  the  leeward  pressure  perturbation  is  equal  and  opposite  to  the  windward  pressure 
perturbation,  as  expected.  At  90  degrees,  a  pressure  perturbation  exist  but  remains  approximately 
constant  along  the  chord. 
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Figure  4.35  Chordwise  Pressure  Distribution  for  Various  Circumferential  Angles. 
Spindle,  Delta=0.160,  AOA=0  Degree. 
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Figure  4.36  Chordwise  Pressure  Distribution  for  Various  Circumferential  Angles. 
Spindle,  Delta=0.160,  AOA=5  Degrees. 


Figure  4.37  Comparison  of  Chordwise  Pressure  Distribution  Along  Selected 
Circumferential  Angular  Positions  for  Two  Angles  of  Attack.  Spindle,  Delta=0.160. 
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Figure  4.38  Comparison  of  Chordwise  Pressure  Distribution  Along  Selected 
Circumferential  Angular  Positions  for  Two  Angles  of  Attack.  Spindle,  Delta=0.160. 


Figure  4.39  Comparison  of  Chordwise  Pressure  Distribution  Along  Selected 
Circumferential  Angular  Positions  for  Two  Angles  of  Attack.  Spindle,  Delta=0.160. 


Figure  4.40  Comparison  of  Chordwise  Pressure  Distribution  Along  a  Selected 
Circumferential  Angular  Position  for  Two  Angles  of  Attack.  Spindle,  Delta=0.160. 


b.  Pressure  Perturbations 

Using  the  PMARC  results  on  Figures  4.37  to  4.40  above,  the  difference  in  pressure 
between  the  body  at  zero  and  five  degrees  angle  of  attack  was  calculated  for  various  circumferential 
positions  and  plotted  against  the  pressure  perturbation  as  computed  with  the  slender  body  theory  at 
the  same  respective  circumferential  positions  and  the  same  change  of  angle  of  attack  (Figures  4.41 
to  4.44). 


Like  for  the  F5  and  triangular  wing  equivalent  body  of  revolution,  there  exist  for 
the  spindle  some  disparities  between  the  two  perturbation  calculations  for  almost  all  the  cases 
considered.  These  disparities  are  larger  towards  the  nose  and  the  tail  and  progressively  reduce  to 
zero  towards  the  shifting  point. 
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Once  again,  it  can  be  observed  that  the  shifting  points  as  calculated  with  the  slender 
body  theory  do  not  always  correspond  to  those  obtained  by  PMARC.  The  differences  for  the  theta 
angles  0,  30,  60,  120,  150,  and  180  degrees  are  6.0%,  3.0%,  28.0%,  19.0%,  4.0%,  and  8.0% 
respectively.  Here,  we  observe  the  largest  deviation  at  a  theta  of  60  degrees.  It  appears  like  the 
deviation  gets  larger  as  theta  approaches  90  degrees  from  either  side. 

For  the  theta  position  of  90  degrees,  we  observe  that  the  perturbations  are  zero 
when  computed  with  the  slender  body  theory,  but  the  PMARC  results  follow  an  elliptic  distribution. 


x/c 


Figure  4.41  Comparison  of  Chordwise  Pressure  Perturbation  at  Selected 
Circumferential  Angular  Positions  for  a  0  to  5  Degree  Change  in  AOA.  Spindle, 
Maximum  Thickness  Ratio  =  0.160. 
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Figure  4.42  Comparison  of  Chordwise  Pressure  Perturbation  at  Selected 
Circumferential  Angular  Positions  for  a  0  to  5  Degree  Change  in  AOA.  Spindle, 
Maximum  Thickness  Ratio  =  0.160. 


Figure  4.43  Comparison  of  Chordwise  Pressure  Perturbation  at  Selected 
Circumferential  Angular  Positions  for  a  0  to  5  Degree  Change  in  AOA.  Spindle, 
Maximum  Thickness  Ratio  =  0.160. 
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Figure  4.44  Comparison  of  Chordwise  Pressure  Perturbation  at  a  Selected 
Circumferential  Angixlar  Position  for  a  0  to  5  Degree  Change  in  AOA.  Spindle, 
Maximum  Thickness  Ratio  =  0. 1 60. 


3.  PMARC  Input  Instructions 

The  steps  used  to  construct  the  PMARC  model  of  the  spindle  significantly  vary  from  those 
used  for  the  F5  and  triangular  wing  equivalent  body  of  revolution.  We  will  present  here  the  second 
method  of  modelling  a  body  of  revolution  in  PMARC.  One  similarity  however  is  the  acquisition  of 
the  cross-sectional  areas.  Although  we  did  not  need  the  assistance  of  the  computer  program  as  with 
the  F5  wing,  an  analytical  expression  was  available  (Equation  (4.7)).  Figure  4.45  presents  an 
extract  of  the  input  file  for  the  spindle.  The  following  subsections  will  review  this  second  approach. 
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a.  Basic  Input  Section 

No  differences  exist  in  the  format  and  content  of  the  basic  input  section  of  this 
method  and  the  first  one  used  for  the  F5  and  triangular  wing  equivalent  body  of  revolution. 

b.  Assembly  and  Components  System  Information 

Once  again,  the  two  methods  present  no  differences  in  this  area  of  the  input  file. 

c.  Patch  and  Section  Coordinate  System  Information 

This  section  is  the  one  that  differs  for  the  two  methods.  Instead  of  defining  the 
coordinates  of  one  circle  and  scaling  the  latter  to  the  dimensions  corresponding  to  the  area  at  a 
prescribed  chord  location,  this  method  redefines  new  circle  coordinates  at  each  of  these  chord 
locations.  Like  for  the  first  method,  the  SECTl  block  is  repeated  each  time  a  circle  is  defined. 
However,  the  reader  will  recall  that  in  that  method  the  field  INMODE  was  set  to  zero  in  all 
succeeding  blocks  to  indicate  that  the  circle  section  defined  in  the  first  block  was  to  be  copied.  It 
was  then  scaled  by  changing  the  SCALE  field. 

In  the  second  method,  INMODE  is  set  at  four  in  all  blocks  to  indicate  that  the  x,  y, 
z  coordinates  of  each  circle  will  be  inputted.  In  this  case  the  field  SCALE  is  always  set  to  zero  since 
the  circles  are  defined  in  global  coordinates. 

In  the  approach  presented  above,  all  the  panel  fields  (TNPS  and  TNPC)  were  set  to 
zero.  This  tells  PMARC  to  use  the  input  section  definitions  to  define  the  panel  edges,  longitudinally 
and  along  the  circumference.  This  was  possible  here  because  all  the  circle  points  were  generated 
so  they  longitudinally  corresponded  to  each  other.  Otherwise,  TNPS  and  TNPC  can  be  prescribed 
and  PMARC  disregards  the  points  coordinates. 

d.  Wake  and  Streamlines 

As  for  the  previous  method,  the  wake  shall  be  turned  off  for  the  spindle.  As  for  the 
streamlines,  they  are  not  required. 
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4.  Global  Coordinate  Circles  Model  Construction  Method  Characteristics 

For  the  studies  performed  in  this  paper,  the  only  advantage  that  this  method  seems  to  present 
is  that  it  allows  PMARC  to  compute  more  easily  the  flow  at  the  pointed  ends.  This  is  said  with  some 
reservation  because  until  we  understand  why  the  behavior  of  the  first  method  is  encountered  in  the 
vicinity  of  the  pointed  ends,  there  is  no  certitude  about  the  actual  superiority  of  the  second  method. 
Maybe  it  is  possible  to  correct  the  input  file  in  the  first  method  to  allow  for  a  smoother  computation 
at  the  ends.  That  would  make  the  first  method  better  all  around  for  the  type  of  studies  performed 
because  the  input  file  for  the  spindle  method  is  about  thirty  times  longer.  Despite  this  apparent 
disadvantage,  this  approach  could  however  be  necessary  in  some  other  configurations. 

5.  PMARC  Output  -  Data  Retrieval 

For  the  purpose  of  our  study,  only  the  pressure  coefficient  distributions  at  various 
circumferential  positions  were  required.  In  order  to  obtain  a  Cp  distribution  all  along  the  chord  for 
the  theta  positions  selected,  one  reading  per  column  of  data  is  required.  One  or  more  columns  of 
data  can  be  skipped  but  care  should  be  exercised  as  to  conserve  the  proper  smoothness  of  the 
distribution.  In  the  DATA6  output  file,  one  column  of  data  details  the  coordinates  and  aerodynamic 
parameters  of  one  circumferential  strip  of  panels  running  all  around  the  body  from  bottom  to  bottom 
as  shown  on  Figure  4.46. 


Figure  4.46  PMARC  Panel  Distribution  on  The  Spindle  Single  Patch  Model. 
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The  circumferential  position  theta  which  need  a  pressure  distribution  analysis  is  translated 
to  a  panel  number  on  the  output  file  by  looking  through  one  column  of  data  for  the  arctangent  of  the 
y  over  -z  ratio  which  will  give  the  angle  needed.  Once  the  position  of  the  panel  relative  to  the  first 
panel  of  that  column  is  determined,  the  corresponding  panel  in  each  column  of  data  will  be 
approximately  at  the  same  angle  theta.  In  this  case,  we  were  able  to  study  the  pressure  distribution 
at  circumferential  positions  dividing  the  half  body  in  six  exact  parts:  0,  30,  60,  90, 120,  150,  and 
180  degrees.  As  mentioned  earlier,  we  did  not  have  to  study  the  mirror  image  position  of  these 
angles  (0,  330,  300,  270, 240,  210  and  180  degrees)  as  the  results  would  be  identical.  Recall  that 
it  is  the  cross-section  point  definitions  that  were  judiciously  made  so  that  we  ended  up  with  control 
points  exactly  at  these  theta  angles. 

E.  CONCLUSION 

From  our  study  of  three  different  bodies  of  revolution,  it  became  apparent  that  PMARC 
generates  pressure  distributions  that  match  the  expected  trends.  When  compared  to  the  slender  body 
theory,  our  results  indicate  good  agreement  near  the  chordwise  location  where  the  pressure 
differential  between  the  zero  and  five  degrees  angle  of  attack  switches  sign.  As  we  proceed  towards 
either  end  of  the  body,  the  difference  between  PMARC  results  and  the  slender  body  theory  in  the 
computation  of  these  pressure  differentials  (or  perturbations)  progressively  gets  larger  to  an 
unacceptable  level  at  the  ends. 

Interpretation  of  the  plots  presented  also  reveals  that  the  shifting  point  is  different  depending 
on  the  method  of  calculation  (PMARC  or  slender  body  theory).  The  differences  however  are  in  most 
cases  in  the  vicinity  of  5%.  In  all  cases,  PMARC  is  expected  to  be  more  precise  than  the  slender 
body  theory. 
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V.  UNSTEADY  CASE  -  OSCILLATION 


A.  INTRODUCTION 


The  intent  for  this  chapter  is  to  validate  the  PMARC  results  for  wings  and  bodies  of 
revolution  in  oscillation.  Time  constraint  prevented  a  complete  analysis  of  the  unsteady  case.  It  was 
possible  to  obtain  results  for  the  F5  and  the  triangular  wings  and  their  respective  equivalent  bodies 
of  revolution  as  well  as  for  the  spindle  for  ten-degree  amplitude  oscillations  about  their  mid-chord 
at  various  reduced  frequencies.  The  results  will  be  presented  with  a  warning  to  the  reader  about  their 
validity.  It  is  left  for  future  work  to  analytically  or  experimentally  verify  these  results. 

B.  F5  WING  AND  EQUIVALENT  BODY  OF  REVOLUTION 


The  F5  wing  was  shown  and  described  in  section  F  of  Chapter  HI,  and  its  equivalent  body 
of  revolution  in  section  B  of  Chapter  IV.  Recall  that  for  the  purpose  of  the  PMARC  modelling,  the 
main  chord  was  set  at  0.6439  unit  and  the  corresponding  maximum  thickness  was  0.0309  unit.  The 
dimensionless  maximum  thickness  ratio  is  then  0.0480.  The  wing  and  its- equivalent  body  of 
revolution  are  reproduced  below  in  Figure  5.1. 


Figure  5.1  F5  Wing  and  its  Corresponding  Equivalent  Area  Body  of  Revolution. 
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This  section  will  first  lead  the  reader  through  the  various  calculations  of  frequency,  time¬ 
stepping  and  other  parameters  related  to  each  oscillatory  case.  The  results  will  then  be  presented 
with  the  lift  and  moment  coefficients  plotted  for  two  cycles  of  oscillation.  Following  the  results,  the 
specific  input  instructions  to  run  an  oscillatory  case  in  PMARC  will  be  described.  Finally,  a 
technique  will  be  shown  to  retrieve  the  numerous  data  required  from  DATA6. 

1.  Oscillatory  Frequency  and  Time-Stepping  Calculations 

Three  reduced  frequencies  (A:)  were  selected  for  the  oscillations:  0.1,  0.5,  and  1.0.  The 
oscillations  had  an  amplitude  of  plus  and  minus  five  degrees  about  the  mid-chord.  Two  full 
oscillatoiy  cycles  were  required.  Here  is  a  sample  calculation  for  k=0.1  and  the  results  for  the  other 
cases.  All  parameters  shown  below  apply  for  both  the  wing  and  the  equivalent  body  of  revolution. 

In  the  following  equations,  c  was  taken  to  be  the  root  chord,  and  the  freestream  velocity  was 
arbitrarily  selected  at  1 1  units  per  second.  Omega  is  the  oscillatory  frequency  in  rad/sec  and  its 
PMARC  notation  is  WRY.  T  is  the  period,  iVthe  number  of  time-steps  and  dt  the  size  of  the  time- 
step  in  seconds.  The  PMARC  notations  of  the  last  two  parameters  are  NTSTPS  and  DTSTEP 
respectively. 


yfc=0. \rad—=  ^--^-^(xi^WRYAjmradls 

U  11 

00 


(5.1) 


271  ,3  57g0sec/m;c/g 

(D  1.7083 


(5.2) 


t=N*dt=NTSTFS*DTSTEP^3.67i0=NTSTPS*0.1^NTSTPS=36.7Ssteps/cycle  (5.3) 


NTSTPS=73.56  time-steps  for  two  cycles 
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As  noticed,  the  number  of  time-steps  (NTSTPS)  times  the  size  of  each  time-step  (DTSTEP) 
must  equal  the  period  (T).  The  combination  selected  is  the  one  which  will  render  a  smooth 
distribution.  In  PMARC,  this  becomes  a  trial  and  error  exercise.  The  following  table  provides  the 
values  of  the  above  parameters  for  all  values  of  reduced  frequency  studied. 


k 

Amplitude 

Osc.  About... 

Omega  (WRY) 

NTSTPS 

DTSTEP 

0.1 

+-  5  Degrees 

Mid-Chord 

1.7083  rad/s 

74 

0.1  sec 

0.5 

+-  5  Degrees 

Mid-Chord 

8.5417  rad/s 

30 

0.05  sec 

1.0 

+-  5  Degrees 

Mid-Chord 

17.0834  rad/s 

37' 

0.02  sec 

Table  5.1  Test  Parameters  for  The  Oscillation  of  The  F5  Wing  Equivalent  Body  of  Revolution. 


2.  Results 

Figures  5.3  to  5.5  below  show  the  combined  CL  and  CM  results  of  the  wing  and  the 
equivalent  body  of  revolution  in  oscillation  at  the  three  reduced  frequencies  listed  in  Table  5.1.  The 
analysis  of  the  relative  maximums  and  minimums  of  the  CL  and  CM  curves  as  a  function  of  the 
reduced  frequency,  as  well  as  the  general  behavior  will  be  left  to  future  worL  It  will  only  be  said 
that  the  lift  and  moment  coefficients  for  the  body  of  revolution  seem  to  behave  as  expected.  The 
lift  vectors  aft  of  the  pivot  point  are  balanced  by  those  behind  it  (Figure  5.2  ).  This  causes  the  net 
value  of  the  lift  to  be  close  to  zero  (zero  in  steady  state)  and  the  moments  to  increasingly  become 
larger  with  angle  of  attack  as  larger  forces  each  side  of  the  pivot  point  and  in  opposite  directions 
torque  the  body.  Some  lift  is  expected  due  to  oscillation  effects.  All  the  oscillation  figures  are 
plotted  against  reduced  time  (tZ»<2r)  which  is  obtained  as  follows: 


<0  =  NTSTPS({)  *  DTSTEP 


Figure  5.2  Typical  Balanced  Pressure  Distribution  Along  The 
Longitudinal  Camber  of  a  Body  of  Revolution. 


Figure  5.3  Comparison  of  The  Lift  and  Moment  Coefficients  in  Unsteady  State  as  a 
Function  of  Reduced  Time.  F5  Wing  and  Equivalent  Body  of  Revolution,  Oscillation 
Amplitude=10  Degrees  About  c/2,  Ujnf^ll,  Omega=1.7083  rad/s,  k=0.1. 
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Figure  5.4  Comparison  of  The  Lift  and  Moment  Coefficients  in  Unsteady  State  as  a 
Function  of  Reduced  Time.  F5  Wing  and  Equivalent  Body  of  Revolution,  Oscillation 
Amplitude=10  Degrees  About  c/2,  Uj„f=l  1,  Omega=8.5417  rad/s,  k=0.5. 
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Figure  5.5  Comparison  of  The  Lift  and  Moment  Coefficients  in  Unsteady  State  as  a 
Function  of  Reduced  Time.  F5  Wing  and  Equivalent  Body  of  Revolution,  Oscillation 
Amplitude=10  Degrees  About  c/2,  Ui„^l  1,  Omega=17.0834  rad/s,  k=1.0. 
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3. 


PMARC  Input  Instructions 


In  the  unsteady  case  like  in  the  steady  case,  the  most  difficult  part  about  the  conception  of 
the  input  file  is  the  geometry  modelling.  Throughout  this  paper  we  reviewed  most  of  the  techniques 
for  simple  wings  and  bodies  of  revolution.  With  this  background,  we  got  virtually  all  that  it  needed 
to  run  an  unsteady  case  in  PMARC.  Only  a  few  fields  have  to  be  changed  to  reflect  the  oscillation 
parameters.  Theses  fields  are  concentrated  in  the  first  two  major  sections  of  the  input  file:  the  Basic 
Input  Section,  and  the  Assembly  and  Component  Coordinate  Section.  Figures  5.6  and  5.7  show 
these  sections  extracted  from  the  input  file  of  the  F5  wing  and  its  equivalent  body  of  revolution 
respectively.  The  fields  modified  to  run  the  oscillations  are  boxed.  The  following  paragraphs  will 
review  these  fields.  The  rest  of  the  input  file  won't  be  discussed  as  there  is  no  difference  with  the 
input  file  of  the  steady  cases  we  ran  previously. 

a.  Basic  Input  Section 

In  this  section,  four  fields  are  directly  related  to  the  oscillatory  case.  They  are: 
NTSTPS,  DTSTEP,  THEMAX,  and  WRY.  Except  for  THEMAX,  definitions  and  calculation 
methods  for  these  fields  were  given  in  subsection  1  above.  THEMAX  is  simply  the  amplitude  of 
the  oscillatory  rotation  about  the  PMARC  global  Y  axis,  in  degree.  Two  other  fields  are  not  inherent 
to  the  oscillation  but  may  need  adjustment  to  allow  a  successful  run.  They  are:  MAXIT  and 
SOLRES.  The  first  one  is  the  limit  on  number  of  solver  iterations.  The  second  is  the  convergence 
criteria  for  the  matrix  solver.  The  recommended  setting  are  150  to  200  for  MAXIT  and  0.0005  for 
SOLRES. 


b.  Assembly  and  Component  Coordinate  System  Information 
In  this  section  of  the  input  file,  only  one  field  requires  to  be  changed  from  the  steady 
case:  ASEMX.  That  indicates  in  PMARC  the  axis  of  rotation  for  the  oscillation.  As  can  be  noticed 
on  Figures  5.6  and  5.7,  this  value  needs  to  be  negative  in  order  to  have  the  rotation  axis  on  the  body. 
That  is  because  the  body  is  defined  in  assembly  coordinates  and  by  putting  a  negative  value  to 
ASEMX  we  are  prescribing  the  origin  of  the  assembly  coordinates  to  begin  somewhere  behind  the 
global  coordinates  origin,  the  latter  being  taken  for  the  axis  of  rotation. 


148 


c.  Customizing  The  Source  Code  If  Failure  to  Run 
Two  other  parameters  require  attention.  However,  they  are  not  found  in  the  input 
file.  They  are  in  the  parameter  file  of  the  source  code  (PARAM.DAT).  One  is  the  number  of  surface 
panels  allowed  (NSPDIM)  and  the  other  one  is  the  number  of  wake  panels  allowed  (NWPDIM). 
NSPDIM  should  at  least  be  set  to  the  number  of  panels  covering  the  surface  studied.  NWPDIM  was 
boosted  up  form  1900  to  5000  in  order  to  get  a  successful  run.  Although  it  has  not  yet  been 
confirmed  or  scientifically  explained,  it  appears  like  this  latter  parameter  should  be  set  at  a  value  at 
least  equal  to  the  product  of  the  total  number  of  spanwise  panels  and  the  number  of  time-steps. 
When  the  parameter  file  is  modified,  the  source  code  has  to  be  recompiled. 

4.  PMARC  Output  -  Data  Retrieval 

In  the  oscillatory  case,  large  files  of  output  data  are  created;  sometimes  up  to  6  Megabytes. 
Unless  the  number  of  time-steps  is  very  small,  the  data  retrieval  should  not  be  done  manually  as  it 
would  become  very  time  consuming.  The  output  file  for  the  oscillation  case  follows  the  same  format 
as  for  the  steady  case.  At  time-step  zero  and  at  the  very  last  time-step,  the  x,  y,  z  coordinates,  the  Cp, 
Vx,  Vy,  Vz,  V,  DUB,  DPHI/DT  and  MACH  number  are  provided  for  the  control  points  of  all  panels. 
For  the  intermediate  time-steps,  the  control  point  coordinates  from  the  leading  edge  are  given  with 
the  section  chord,  CL,  CD,  CY,  Cm,  Cn,  and  Cl.  Figure  5.8  shows  selected  portions  of  a  thirty  page 
output  file. 

Not  shown  on  Figure  5.8  are  the  total  coefficients  for  each  time-steps.  These  are  the  values 
we  needed  in  our  study.  More  specifically  CL  and  CM.  In  order  to  efficiently  extract  them,  small 
programs  written  in  the  UNIX  operating  system  AWK  language  were  written.  The  data  retrieval 
becomes  relatively  simple  as  the  output  file  repeats  the  total  coefficient  block  under  the  same  format 
after  each  time-step.  Thus,  these  programs  look  for  a  prescribed  recurring  feature,  and  from  it 
counts  an  indicated  number  of  lines  and  then  take  the  desired  value  in  the  prescribed  field  of  that 
line.  In  our  case,  the  program  was  searching  for  the  words  WIND  AXES  from  which  it  counted  32 
lines  and  then  stored  the  first  field  for  the  CL  value  or  the  fourth  for  the  CM  value.  This  process  was 
repeated  at  each  time-step.  These  geometry  specific  programs  can  be  found  in  Appendix  H. 
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Figure  5.8  Selected  Portion  of  a  PMARC  Unsteady  Case  Output  File  Showing  The 
Data  Produced  at  The  Initial,  First  and  Last  Time-Step. 


c. 


TRIANGULAR  WING  AND  EQUIVALENT  BODY  OF  REVOLUTION 


The  triangular  wing  was  shown  and  described  in  section  G  of  Chapter  m,  and  its  equivalent 
body  of  revolution  in  section  C  of  Chapter  IV.  For  these  cases,  the  main  chord  and  longitudinal  axis 
respectively  were  made  unity.  Like  for  the  steady  case,  two  thickness  ratios  will  be  treated  in 
oscillation:  0.098,  and  0.196.  The  wing  and  its  equivalent  body  of  revolution  are  reproduced  below 
in  Figure  5.9.  Then  the  oscillation  parameter  calculations,  the  oscillation  results,  and  the  PMARC 
specifics  are  presented. 


Figure  5.9  PMARC  Top  and  Side  Representation  of  The  Triangular  Wing 
and  Its  Corresponding  Equivalent  Area  Body  of  Revolution.  Maximimi 
Thickness  Ratio  =  0.098.  (Each  View  Not  at  Same  Scale). 


1.  Oscillatory  Frequency  and  Time-Stepping  Calculations 

The  same  three  reduced  frequencies  (k)  as  for  the  F5  were  selected  for  these  oscillation  runs: 
0.1, 0.5,  and  1 .0.  The  movement  had  the  same  plus  and  minus  five  degrees  amplitude  about  the  mid¬ 
chord,  and  two  full  oscillatory  cycles  were  operated.  Again,  here  is  a  sample  parameters  calculation 
for  k=0.1  and  the  results  for  the  other  cases.  All  parameters  shown  below  apply  for  both  the  wing 
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and  the  equivalent  body  of  revolution,  c  is  the  root  chord  (or  longitudinal  axis),  omega  is  the 
oscillatory  frequency  in  rad/sec  (WRY  in  PMARC),  T  is  the  period,  N  the  number  of  time-steps  and 
dt  the  size  of  the  time-steps  in  seconds.  The  PMARC  notations  for  the  last  two  parameters  are 
NTSTPS  and  DTSTEP  respectively.  The  freestream  velocity  was  arbitrarily  selected  at  1 1  unit  per 
second. 


=.  «=fFi?r=l.l  rad/s 

11 


(5.4) 


=YY=5.71  SQc/cycle 


(5.5) 


T=N*dt=NTSTPS*DTSTEP  -  5.7UNTSTPS*0.2  -  NTSTPS=2i.55  stepsicycle 

(5.6) 


NTSTPS=57.1  time-steps  for  two  cycles 


As  mentioned  in  the  previous  section,  the  number  of  time-steps  (NTSTPS)  times  the  size 
of  each  time-step  (DTSTEP)  must  equal  the  period  (T).  The  combination  selected  is  the  one  which 
will  render  a  smooth  distribution  in  the  PMARC  output.  The  following  table  now  provides  the  value 
of  the  parameters  shown  above  for  all  values  of  reduced  frequency  studied.  Note  that  all  the 
parameters  not  only  apply  for  both  the  wing  and  the  equivalent  body  of  revolution,  but  also  for  the 
two  thickness  distributions  considered. 


k 

Amplitude 

Osc.  About... 

Omega  (WRY) 

NTSTPS 

DTSTEP 

0.1 

+-  5  Degrees 

Mid-Chord 

1.1  rad/s 

58 

0.2  sec 

0.5 

+-  5  Degrees 

Mid-Chord 

5.5  rad/s 

23 

0.1  sec 

1.0 

+-  5  Degrees 

Mid-Chord 

1 1 .0  rad/s 

23 

0.05  sec 

Table  5.2  Test  Parameters  for  The  Oscillation  of  The  Triang.  Wing  Equivalent  Body  of  Revolution. 
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2. 


Results 


Figures  5.10  to  5.12  below  show  the  combined  CL  and  CM  results  of  the  wing  and  the 
equivalent  body  of  revolution  of  thickness  ratio  0.098  in  oscillation  at  the  three  reduced  frequencies 
listed  in  Table  5.2.  Figures  5.13  to  5.15  do  the  same  thing  for  a  thickness  ratio  of  0.196.  The 
analysis  of  the  relative  maximums  and  minimums  of  the  CL  and  CM  curves  as  a  function  of  the 
reduced  frequency,  as  well  as  the  general  behavior  will  be  left  to  future  work.  It  will  only  be  said 
that  the  lift  coefficients  for  the  body  of  revolution  seem  to  behave  as  expected  only  for  the  low 
reduced  frequency,  contrary  to  the  F5  where  their  behavior  seemed  legitimate  for  all  reduced 
frequencies.  As  illustrated  in  section  B-2  of  this  chapter,  the  lift  vectors  in  front  of  the  pivot  point 
are  theoretically  balanced  by  those  behind  it  (Figure  5.2).  This  should  cause  the  net  value  of  the  lift 
to  be  close  to  zero  at  any  angle  of  attack.  Some  lift  is  expected  due  to  the  oscillation  effects.  Also, 
as  observed  in  om  results,  the  moments  should  increasingly  become  larger  with  angle  of  attack.  As 
detailed  in  the  previous  section,  all  the  figures  are  plotted  against  the  reduced  time  which  is  the  ratio 
of  the  product  of  freestream  velocity  and  time  over  the  main  chord. 

3.  PMARC  Input  Instructions 

Except  for  the  values  of  the  parameters,  there  is  no  difference  between  the  input  files  for  the 
oscillations  of  the  triangular  wing  and  its  equivalent  body  of  revolution  and  those  of  the  F5  wing 
and  its  equivalent  body  of  revolution  studied  in  the  previous  section.  One  input  file  for  the  triangular 
wing  and  one  for  its  equivalent  body  of  revolution  for  the  reduced  frequency  of  0.1  are  presented 
below.  The  reader  should  consult  section  B-3  of  this  chapter  for  details. 

4.  PMARC  Output-  Data  Retrieval 

As  no  difference  exists  with  the  data  retrieval  techniques  used  in  the  previous  section,  the 
reader  should  refer  to  the  latter  for  details  on  the  program  required  to  extract  the  data  off  DATA6. 
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Figure  5.10  Comparison  of  The  Lift  and  Moment  Coefficients  in  Unsteady  State  as  a 
Frmction  of  Reduced  Time.  Triangular  Wing  of  AR=2/3  and  Equivalent  Body  of 
Revolution  of  Delta=0.098,  Oscillation  Amplitude=10  Degrees  About  c/2,  Ui„f=l  1, 
Omega=l.l  rad/s,  k=0.1. 


z 


B  0 
-1 


0  5  10  15  20  25  30 

Reduced  Time 

Figure  5.11  Comparison  of  The  Lift  and  Moment  Coefficients  in  Unsteady  State  as  a 
Fimction  of  Reduced  Time.  Triangular  Wing  of  AR=2/3  and  Equivalent  Body  of 
Revolution  of  Delta=0.098,  Oscillation  Amplitude=10  Degrees  About  c/2,  Uinf=l  1, 
Omega=5.5  rad/s,  k=0.5. 


156 


2 


i  0 


-2 

0  2  4  6  8  10  12  14 

Reduced  Time 

Figure  5.12  Comparison  of  The  Lift  and  Moment  Coefficients  in  Unsteady  State  as  a 
Fimction  of  Reduced  Time.  Triangular  Wing  of  AR=2/3  and  Equivalent  Body  of 
Revolution  of  Delta=0.098,  Oscillation  Amplitude=10  Degrees  About  c/2,  Ui„f=l  1, 
Omega=l  1.0  rad/s,  k=1.0. 
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Figure  5.13  Comparison  of  The  Lift  and  Moment  Coefficients  in  Unsteady  State  as  a 
Function  of  Reduced  Time.  Triangular  Wing  of  AR=2/3  and  Equivalent  Body  of 
Revolution  of  Delta=0.196,  Oscillation  Amplitude=10  Degrees  About  c/2,  Ujnf^l  1, 
Omega=l.l  rad/s,  k=0.1. 
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Figure  5.14  Comparison  of  The  Lift  and  Moment  Coefficients  in  Unsteady  State  as  a 
Fimction  of  Reduced  Time.  Triangular  Wing  of  AR=2/3  and  Equivalent  Body  of 
Revolution  of  Delta=0. 1 96,  Oscillation  Amplitude=l  0  Degrees  About  c/2,  Uinf=l  1 , 
Omega=5.5  rad/s,  k=0.5. 
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Figure  5.15  Comparison  of  The  Lift  and  Moment  Coefficients  in  Unsteady  State  as  a 
Fimction  of  Reduced  Time.  Triangular  Wing  of  AR=2/3  and  Equivalent  Body  of 
Revolution  of  Delta=0.196,  Oscillation  Amplitude=10  Degrees  About  c/2,  Ui„f^l  1, 
Omega=l  1.0  rad/s,  k=1.0. 
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Figure  5.17  PMARC  Input  File  for  The  Triangular  Wing  Equivalent  Body  of 
Revolution  in  Oscillation  About  The  Mid-Chord  for  k=0.1. 
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4SECT1 

STX*  0.7750, 

STY-  0.0000, 

STZ-  0.0000, 

SCALE-  0.0530, 

ALF*  0.0, 

THETA*  0 . 0 

t 

XNMODE-  0 , 

TNODS*  0, 

TNPS*  0,  TINTS* 

0, 

&END 

tSECTl 

STX=  0.8000, 

STY*  0.0000, 

STZ*  0.0000, 

SCALE*  0.0516, 

ALF=  0.0, 

THETA*  0 . 0 

INMODE=  0, 

TNODS*  0, 

TNPS*  0,  TINTS* 

0, 

fcEND 

fiSECTl 

STX=  0.8250, 

STY*  0.0000, 

STZ*  0.0000, 

SCALE*  0.0498, 

ALF=  0.0, 

THETA*  0 . 0 

INMODE=  0, 

TNODS*  0, 

TNPS*  0,  TINTS* 

0  t 

&END 

iSECTl 

STX=  0.8500, 

STY*  0.0000, 

STZ*  0.0000, 

SCALE*  0.0475, 

ALF=  0.0, 

THETA*  0 . 0 

INMODE=  0, 

TNODS*  0, 

TNPS*  0,  TINTS* 

0/ 

&EMD 

&SECT1 

STX=  0.8750, 

STY*  0.0000, 

STZ*  0.0000, 

SCALE*  0.0446, 

ALF=  0.0, 

THETA*  0.0 

INMODE==  0, 

TNODS*  0, 

TNPS*  0,  TINTS* 

0, 

4END 

&SECT1 

STX=  0.9000, 

STY*  0.0000, 

STZ*  0.0000, 

SCALE*  0.0410, 

AIJ'=  0.0, 

THETA*  0.0, 

INMODE=  0 , 

TNODS*  0, 

TNPS*  0,  TINTS* 

0, 

&END 

&SECT1 

STX*  0.9250, 

STY*  0.0000, 

STZ*  0.0000, 

SCALE*  0.0365, 

AU=  0.0, 

THETA*  0.0, 

INMODE=  0 , 

TNODS*  0, 

TNPS*  0,  TINTS* 

0/ 

&END 

&SECT1 

STX=  0.9500, 

STY*  0.0000, 

STZ*  0.0000, 

SCALE*  0.0306, 

ALF-  0.0, 

THETA*  0.0, 

INMODE=  0, 

TNODS*  0, 

TNPS*  0 ,  TINTS* 

0/ 

t.END 

&SECT1 

STX=  0,9625, 

STY-  0.0000, 

STZ*  0.0000, 

SCALE*  0.0269, 

ALF=  0.0, 

THETA*  0.0, 

INMODE=  0 , 

TNODS-  0, 

TNPS*  0 ,  TINTS* 

0, 

&END 

tSECTI 

STX=  0.9750, 

STY-  0,0000, 

STZ*  0.0000, 

SCALE*  0.0222, 

ALr=  0.0, 

THETA*  0.0, 

INMODE=  0 , 

TNODS-  0, 

TNPS*  0,  TINTS- 

0, 

4END 

tSECTl 

STXs=  0.9813, 

STY-  0.0000, 

STZ*  0.0000, 

SCALE-  0.0194, 

ALF«  0.0, 

THETA-  0.0, 

INMODE=  0 , 

TNODS*  0, 

TNPS*  0 ,  TINTS* 

0, 

SEND 

fiSECTl 

STX*  0.9876, 

STY-  0.0000, 

STZ-  0.0000, 

SCALE*  0.0159, 

ALT*  0.0, 

THETA*  0.0, 

INMODE=  0, 

TNODS*  0, 

TNPS*  0,  TINTS* 

0, 

SEND 

tSECTl 

STX«  0.9938, 

STY-  0.0000, 

STZ*  0.0000, 

SCALE-  0.0113, 

ALF=  0.0, 

THETA-  0.0, 

INMODE=  0 , 

TNODS*  0, 

TNPS*  0,  TINTS* 

0, 

SEND 

&SECT1 

STX=  0.9969, 

STY-  0.0000, 

STZ*  0.0000, 

SCALE-  0.0080, 

ALF=  0.0, 

THETA*  0.0, 

INMODE*  0 , 

TNODS*  0, 

TNPS-  0,  TINTS* 

Or 

SEND 

fcSECTl 

STX=  1.0000, 

STY*  0.0000, 

STZ*  0.0000, 

SCALE*  0.0000, 

ALF*  0,0, 

THETA-  0.0, 

INMODE*  0 , 

TNODS*  5, 

TNPS*  22,  TINTS= 

=  2, 

SEND 

&BPNODE 

TNODE=3, 

TNPC-0,  TINTC-0,  SEND 

&WAKE1  IDWAK-0,  irLXW=l ,  ITRFTZ-  0,  INTRW=  0,  ^  fcEND 

WING  WAKE 


4WAKE2  KWPACH=2,  KWSIDE^4,  KWLINE=0,  KWPAN1=0, 

KWPAN2=0,  NODEW=0,  INITIAL=0 ,  SEND 


SWAKE2 

KWPACH-1, 

KWSIDE*2, 

KWLINE*0, 

KWPANl-0, 

KWPAN2-0, 

NODEW-5 , 

INITIAL-0, 

SEND 

SONSTRM 

NONSL 

=0,  KPSL 

- 

SEND 

SBLPARAM 

RN  *35755000, 

Vise 

=  0. cools'? 

,  NSLBL  =  1, 

SEND 

SVSl 

NVOLR* 

0, 

NVOLC* 

0, 

SEND 

SVS2 

X0= 

-0.1000, 

Y0  = 

1.5000, 

Z0= 

-0.1000, 

INTVSR*  1, 

SEND 

SVS3 

XI* 

-0.1000, 

Yl= 

1.5000, 

Zl* 

0.1000, 

NPTl*  0, 

SEND 

SVS4 

X2  = 

-0.1000, 

Y2* 

1.5000, 

Z2=  . 

0.1000, 

NPT2*  20, 

SEND 

SVS5 

X3* 

1.1000, 

Y3* 

1.5000, 

Z3* 

-0.1000, 

NPT3*  25, 

SEND 

SVS6 

XRO* 

2.0000, 

YRO* 

2.0000, 

ZRO* 

0.0000, 

INTVSC*  1, 

SEND 

SVS7 

XRl* 

4.0000, 

YRl* 

2.0000, 

ZRl* 

0.0000, 

XR2* 

2.0000, 

YR2* 

2.0000, 

ZR2* 

1.0000, 

SEND 

SVS8 

Rl* 

0.1000, 

R2  = 

1.0000, 

PHIl* 

0.0, 

PHI2-360.0, 

SEND 

SVS9 

NRAD* 

5, 

NPHI* 

12, 

NLEN* 

3, 

SEND 

SSLINl 

NSTLIN^ 

=0, 

SEND 

SSLIN2 

SXO* 

-2.0000, 

SYO* 

1.0000, 

SZO* 

-0.5000, 

Figure  5.17  Cont'd  PMARC  Input  File  for  The  Triangular  Wing 
Equivalent  Body  of  Revolution  in  Oscillation  About  The  Mid-Chord  for 
k=0.1. 
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D.  SPINDLE 


The  spindle  was  shown  and  described  in  section  D  of  Chapter  IV.  The  main  chord  (or 
longitudinal  axis)  was  made  unity.  The  thickness  ratio  considered  is  0.160.  The  spindle  is 
reproduced  below  in  Figure  5.18.  Then  the  oscillation  parameter  calculations,  the  oscillation 
results,  and  the  PMARC  specifics  are  presented. 


Figure  5.18  PMARC  Representation  of  The  Spindle. 


1.  Oscillatory  Frequency  and  Time-Stepping  Calculations 

This  time,  only  two  reduced  frequencies  (A:)  will  be  studied:  0.1  and  0.5.  The  oscillation  had 
the  same  plus  and  minus  five  degrees  amplitude  about  the  mid-chord,  and  two  full  oscillatory  cycles 
were  also  computed.  Again,  here  is  a  sample  parameters  calculation  for  k=0.1  and  the  results  for 
the  other  case,  c  is  the  root  chord  (or  longitudinal  axis),  omega  is  the  oscillatory  frequency  in 
rad/sec  (WRY  in  PMARC),  T  is  the  period,  N  the  number  of  time-steps  and  dt  the  size  of  the  time- 
steps  in  seconds.  The  PMARC  notations  for  the  last  two  parameters  are  NTSTPS  and  DTSTEP 
respectively.  The  ffeestream  velocity  was  arbitrarily  selected  at  100  unit  per  second  for  the  case  of 
k=0.1  and  at  1 1  for  the  case  of  k=0.5. 
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=*^  a)=priZy=10.0  rcuHs 


(5.7) 


U  100 


T~~=0.6283  secicycle  (5.8) 

T=N*dt=NTSTPS*DTSTEP  =»  0.6283=ATOrP5*0.025  =*  NTSTPS=25A3  steps/cycle 

(5.9) 


NTSTPS=50.27  time-steps  for  two  cycles 

As  mentioned  before,  the  number  of  time-steps  (NTSTPS)  times  the  size  of  each  time-step 
(DTSTEP)  must  equal  the  period  (T).  The  combination  selected  is  the  one  which  will  render  a 
smooth  distribution  in  the  PMARC  output.  The  following  table  now  provides  the  value  of  the 
parameters  shown  above  for  the  two  values  of  reduced  frequency  studied.  Recall  that  the  freestream 
velocity  was  100  unit  per  second  for  A:=0.1  and  1 1  unit  per  second  for  ^0.5.  This  selection  was 
necessary  in  order  to  keep  the  time  step  size  from  getting  too  small. 


k 

Amplitude 

Osc.  About... 

Omega  (WRY) 

NTSTPS 

DTSTEP 

0.1 

+-  5  Degrees 

Mid-Chord 

10.0  rad/s 

51 

0.025  sec 

0.5 

+-  5  Degrees 

Mid-Chord 

5.5  rad/s 

115 

0.02  sec 

Table  5.3  Test  Parameters  for  The  Oscillation  of  The  Spindle. 


2.  Results 

Figures  5.19  and  5.20  below  show  the  combined  CL  and  CM  results  for  the  two  reduced 
frequencies  considered.  The  analysis  of  the  relative  maximums  and  minimums  of  the  CL  and  CM 
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curves  as  a  function  of  the  reduced  frequency,  as  well  as  the  general  behavior  will  be  left  to  future 
work.  It  can  be  said  however  that  the  behavior  of  the  lift  and  moment  coefficients  warrants 
interrogation  as  to  their  validity.  As  illustrated  in  the  section  B-2  of  this  chapter,  the  lift  vectors  in 
front  of  the  pivot  point  are  theoretically  balanced  by  those  behind  it  (Figure  5.2).  This  should  cause 
the  net  value  of  the  lift  to  be  close  to  zero  at  any  angle  of  attack  and  the  moment  to  vary  over  wider 
ranges.  Some  lift  is  expected  due  to  the  oscillation  effects.  Although  the  scaling  of  the  plot  varies 
widely,  we  still  observe  some  lift,  but  more  importantly,  an  equal  moment  variation.  As  detailed 
in  section  B  of  this  chapter,  all  the  oscillation  figures  are  plotted  against  the  reduced  time  which  is 
the  ratio  of  the  product  of  the  freestream  velocity  and  time  over  the  main  chord. 


Figure  5.19  Lift  and  Moment  Coefficient  in  Unsteady  State  as  a  Function  of  Reduced 
Time.  Spindle,  Delta=0.160,  Oscillation  Amplitude  =  10  Degrees  About  c/2,  Uj^lOO, 
Omega=10.0  rad/s,  k=0.1. 
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X  10“® 


Figure  5.20  Lift  and  Moment  Coefficient  in  Unsteady  State  as  a  Function  of  Reduced 
Time.  Spindle,  Delta=0.160,  Oscillation  Amplitude  =  10  Degrees  About  c/2,  Ui„f=l  1, 
Omega=5.5  rad/s,  k=0.5. 


3.  PMARC  Input  Instructions 

Except  for  the  values  of  the  parameters,  there  is  no  difference  between  the  input  files  for  the 
oscillations  of  the  spindle  and  those  of  the  previous  cases.  An  extract  from  one  input  file  for  the 
spindle  is  presented  below.  The  reader  should  consult  section  B-3  of  this  chapter  for  further  details. 

4.  PMARC  Output  -  Data  Retrieval 

As  no  difference  exists  with  the  data  retrieval  techniques  used  in  the  previous  oscillation 
sections,  the  reader  should  refer  to  the  latter  for  details  on  the  program  required  to  extract  the  data 
offDATA6. 
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Figure  5.21  Extract  From  The  PMARC  Input  File  for  a  Spindle  in  Oscillation  About 
The  Mid-Chord,  Amplitude=10  Degrees,  0mega=10.0  rad/s,  k=0.1. 
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E.  CONCLUSION 


In  this  section,  the  lift  and  moment  coefficients  variations  at  different  reduced  frequencies 
were  presented  for  the  F5  wing,  the  triangular  wing,  their  respective  equivalent  body  of  revolution, 
and  the  spindle.  The  magnitude  and  relative  behavior  of  the  results  seem  legitimate  in  some  cases 
but  require  more  study  in  others. 

In  the  event  that  the  results  are  indeed  erratic,  one  should  not  exclude  the  software  as  a 
potential  source  of  error.  However,  the  input  files  could  also  be  introducing  unrealistic  behavior. 
One  should  then  look  for  fields  not  yet  discussed,  or  in  the  source  code  and  its  parameter  sub¬ 
routines. 
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VI.  SUMMARY  AND  CONCLUSION 


It  was  the  objective  of  this  thesis  to  contribute  to  the  development  of  a  new  approach  to 
determine  the  steady  and  unsteady  airloads  on  missiles  and  on  low-aspect  ratio  aircraft 
configurations  typically  used  for  high-performance  fighter  aircraft.  This  approach  uses  a  second- 
order  approximation  to  the  steady  and  oscillatory  lifting  flow  past  low  aspect  ratio  wings  which 
exhibits  the  same  fundamental  division  as  the  slender-body/slender-wing  theory:  incompressible 
cross-flow  and  a  Mach-number  dependent  spatial  influence.  To  test  the  range  of  validity  of  this 
second-order  theoiy,  it  was  necessary  to  compare  it  with  an  exact  three-dimensional  flow  solution, 
and  a  low  speed  three-dimensional  flow  panel  code  named  PMARC  developed  by  NASA-Ames  was 
selected  for  that  goal. 

Therefore,  PMARC  was  first  validated  by  comparing  it  with  other  known  solutions.  For 
this  purpose  lifting  flow  over  swept  back,  rectangular,  and  delta  wings  was  studied  and  the  PMARC 
predictions  were  found  to  be  in  good  agreement.  A  detailed  analysis  of  the  F5  wing  was  then 
performed,  and  a  program  was  written  to  compute  its  equivalent  area  body  of  revolution.  The 
lifting  pressure  distributions  were  again  found  to  agree  well  with  the  measufed  ones.  However,  a 
discrepancy  was  detected  where  the  PMARC  predictions  of  the  non-lifting  pressure  distributions  on 
delta  wings  at  zero  angle  of  attack  were  compared  with  linearized  theory  and  with  the  Oswatitsch- 
Keune  low  aspect  ratio  wing  theoiy.  PMARC  predicts  velocity  perturbations  which  do  not  decrease 
with  decreasing  local  chordwise  thickness  ratio  of  the  wing.  Interestingly,  no  similar  discrepancy 
was  found  on  lifting  and  non-lifting  rectangular  wings  when  thickness  ratio  decreases  toward  the 
tips.  Similarly,  the  PMARC  predictions  for  lifting  and  non-lifting  bodies  of  revolution  generally 
compare  well  with  slender  body  theory.  Some  differences  in  the  magnitude  of  the  pressure 
perturbation  exist  towards  the  tips  of  the  bodies  of  revolution,  which  are  in  some  cases  larger  than 
10  percent.  The  chordwise  location  where  the  pressure  perturbation  switches  sign  also  varies  from 
the  PMARC  to  the  slender  body  theory,  but  never  by  more  than  5  percent.  Finally,  PMARC  was 
applied  to  the  computation  of  the  forces  and  moments  generated  on  oscillatory  wings  and  bodies  of 
revolution.  However,  time  constraint  prevented  any  analytical  or  experimental  comparison. 
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In  light  of  the  above  summary,  our  contribution  to  the  development  of  a  second-order 
approximation  to  the  steady  and  oscillatory  lifting  flow  past  low  aspect  ratio  wings  may  appear  dim 
as  we  were  not  able  to  obtain  correct  velocity  distributions  around  delta  wings  at  zero  angle  of 
attack  using  PMARC.  However,  through  this  difficulty  we  fulfilled  one  of  our  objectives  which 
was  to  validate  the  software.  It  proved  good  in  all  other  areas.  It  remains  to  work  towards  correcting 
PMARC  or  the  procedure  that  lead  to  the  erroneous  results.  Otherwise  it  will  be  necessary  to 
establish  the  validity  of  another  software. 
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APPENDIX  A.  F5  WING  SPATIAL  INFLUENCE  FLOW  CALCULATION 


This  appendix  shows  all  the  intermediate  calculations  leading  to  the  solution  of  the  spatial  influence 
incompressible  flow  for  the  F5  wing  studied  in  section  F  of  chapter  m.  The  form  of  the  spatial  influence 
flow  equation  solved  is  introduced  as  Equation  (3.23)  and  is  reproduced  below.  The  solution  was  obtained 
for  nine  equally  spaced  chord  locations.  The  software  utilized  to  perform  the  analytical  calculations  is  called 
MAPLE  (Ref  [16]). 


1  Q(0)  Q(x)  1 

Q  =  -  -^[QJ0)ln2x.QJx^)ln2(x^-x)] 

4-11  X  X  -X  4tc  " 


— f  ‘  e^a)i»2a-x)di 

471*' 0  471  Jjt 

(3.23) 
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>Q:=71.2901*xi''10-352.0008*xi'"9+743.1697*xi^8-875.6007*xi''7+630.8471*xi'^6-2 
86.6580*xi''5+81 .701 9*xiM-14.1307*xl''3+1.41 69*xi''2-0.0348*xi; 

Q  ;=  71.2901  -  352.0008  +  743.1697  -  875.6007  4’ 

+  630.8471  -  286.6580  +  81.7019  i,*  -  14.1307 

+  1.4169  ^^-.0348^ 

>  clQ1:=diff(Q,xi); 

dQl  :=  712.9010  -  3168.0072 +  5945.3576 

-  6129.2049  +  3785.0826  -  1433.2900 

+  326.8076  -  42.3921  ^  +  2.8338  ^  -  .0348 

>  clQ2:=diff(dQ1,xi); 

dQ2  ;=  6416.1090  ^  -  25344.0576  ^  +  41617.5032 
-  36775.2294  ^  +  18925.4130  -  5733.1600 

+  980.4228  ^  -  84.7842  4  +  2.8338 

>  dQ3:=diff(dQ2,xi); 

dQ3  ;=  51328.8720  ^  -  177408.4032  +  249705.0192 

-  183876.1470  +  75701.6520  - 17199.4800  ^ 

+  1960.8456  ^  -  84.7842 

>  ans:=subs(xi=:0,dQ1); 


>  ans:=subs(xi=1,dQ1); 


>  ans:=subs(xi=0,dQ2); 


ans  :=  -.0348 


ans  :=  .0536 


ans  :=  2.8338 
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>  ans:=subs(xi=1,dQ2); 


ans  :=  5.0506 

>  term1  :=(1/(4*Pi))*((0.0348/x)+(0.0536)/(1-x)); 

1  1 

.0348  -  +  .0536 -  ' 

,1  X  I  -  X 

term!  := - 

4  71 


>  term2:=:(1/(4*PI))*(2.8338*(ln(2*x))+5.0506*ln(2*(1  -x))); 

^  1  2.8338  ln(2;c)  + 5.0506  ln(2- 2 x) 

termz  := - — - 

4  7C 


>  term3:=dQ3*ln(2*(x-xi)); 

terms  :=  (51328.8720  -  177408.4032  +  249705.0192 

-  183876.1470  +  75701.6520  -  17199.4800 

+  1960.8456  ^-84.7842)  ln(2ji:- 2^)  - 

>  term4:=dQ3*ln(2*(xi-x));  ~~  ' 

term4  ;=  (51328.8720  ^  -  177408.4032  'C  +  249705.0192 

-  183876.1470  +  75701.6520  f  -  17199.4800 

+  1960.8456 ^ -  84.7842) ln(2 %-2x) 

>term5:=(1/(4*Pi))*lnt{term3,xi=0.0..x); 


termS  := 


f  Cx 

=  ^  ( 

4  ^ 

*/  _ 


51328.8720  -  177408.4032 


+  249705.0192  -  183876.1470  +  75701.6520 

-  17199.4800  ^  +  1960.8456  %  -  84.7842)  ln(2  jc  -  2 


>  term5:=(1/(4*Pi))*int(term3,xl=0.0..x); 


1 

terms  :=  -  (84.78420000  x  +  83970.10713  x^  -  1470.634200 x^ 
4 

+  10510.79333  x^  -  39427.94375  x'* 

-  36775.22940  ln(2.  x)  x^  -  101962.8828  x* 

+  6571 3.52078  x’  -  25344.05760  ln(2.  x)  x' 

-  84.78420000  ln(2.  x)  x  +  41617.50320  ln(2.  x)  xS 

-  17438.06768  x*  +  980.4228000  ln(2.  x)  x^ 

-  5733.160000  ln(2.  x)  x^  +  6416.109000  ln(2.  x)  x* 

+  18925.41300  ln(  2.  x)x^)/jt 


>  term6:=(1/(4*Pi))*lnt(term4,xi=x..1 ); 

\fn 

term6  :=  — 

4 


(51328.8720  -  177408.4032 


J 

V  ^ 


+  249705.0192  -  183876.1470  +  75701.6520  ^ 


-  17199.4800  'C  +  1960.8456  ^  -  84.7842)  ln(2  ^  -  2  x)  ^ 

\ 


) 

>  term6:=(1/(4*Pi))*int(term4,xi=x..1 ); 
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1 

term6  :=  -  (-6416.109000  ln(2.  -  2.  x)  x*  +  17438.06768  x* 

-  329.3686900  x  -  41617.50320  ln(2.  -  2.  x)  x« 

+  3594.581720  x^  -  980.4228000  ln(  2.  -  2.  x)  x^ 

+  62238.41 350  x^  -  1 1 5054.2845  x5  +  1 24098.8859  x*^ 

-  19868.30010x3-  18925.41300 ln( 2.- 2.x)x^ 

+  36775.22940  ln(2.  -  2.  x)  x^ 

+  5733.160000  ln(2.  -  2.  x)  x^  -  72129.62978  x’ 

+  84.78420000  lii(2.  -  2.  x)  x  -i- 1 1 .63426214 
+  2.216800000  ln(2.  -  2.  x)  +  25344.05760  ln(2.  -  2.  x)  x’ 

)ln 

>  uoverUinf:=tenn1-term2-term5+term6; 


uoverUinf:= 


1  1 

.0348  —  -I-  .0536 - 

^ _ X _ 1  -X 

4  7t 


1  2.8338  1ii(2x)-h  5.0506  ln(2- 2 x)  1 
4  7t  4 

84.78420000  X -I-  83970.10713  x^  -  1470.634200  x^ 
+  10510.79333  x^  -  39427.94375  ^ 

-  36775.22940  ln(2.  x)  x^  -  101962.8828  x*^ 

+  65713.52078  x’  -  25344.05760  ln(2.  x)  x’ 

-  84.78420000  ln(2.  x)  x  +  41617.50320  ln(2.  x)  x« 
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-  17438.06768  ;c»  +  980.4228000  ln(2.  x) 

-  5733 . 1 60000  ln(  2.  x )  +  54  j  ^  j  99000  ln(  2.  jc )  jc« 

+  18925.41300  ln( 2.  ;c)x'*)/7C  +  ^( 

-6416.109000  ln(2.  -  2.  x)  jc*  +  17438.06768  x* 

-  329.3686900  X-  41617.50320  ln(2.  -  2.  x)  x« 

+  3594.581720x2  -  980.4228000  ln(  2.  -  2.  x)x2 

+  62238.41 350  x^  -  1 15054.2845  x^  +  124098.8859  x‘ 

-  19868.30010x2-  18925.41300  ln( 2.  -  2.x)x^ 

+  36775.22940  ln(2.  -  2.  x)  x^ 

-I-  5733.160000  ln(2.  -  2.  x)  x^  -  72129.62978  x’ 

+  84.78420000  ln(  2.  -  2.  x)  x  -1- 1 1 .63426214 

-I-  2.216800000  ln(2.  -  2.  x)  +  25344.05760  ln(2.  -  2.  x) 


>  xvec:=array(1 ..  10); 
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>  for  i  from  1  to  10  do  uovrUinf[i]:s::evalf(subs(x=xvec[i],uoverUinf))  od; 

uovrUinf^  ^  ^  :=  -.0573921826 
uovrUinf^^^  :=  .0346121370 
uovrUinf^^^  :=  .0017849588 
uovrUinf  :=  -.00548285722 

L  ^  J 

uovrUinf^^^  := -.01293665332 
uovrUinf^^^  :=  -.08484371084 
uovrUinf  :=  -.1297571415 

[7  ] 

uovrUinf,  1 1 89657140 

[  8  ] 

uovrUinf :=  -.0926006053 

Error,  division  by  zero 

> 
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APPENDIX  B.  TRIANGULAR  WING  SPATIAL  INFLUENCE  FLOW 
CALCULATION 

This  appendix  shows  all  the  intermediate  calculations  leading  to  the  solution  of  the  spatial  influence 
incompressible  flow  for  the  triangular  wing  studied  in  section  G  of  chapter  ID.  Two  forms  of  the  spatial 
influence  flow  equations  were  solved  and  were  introduced  as  Equation  (3.23)  and  Equation  (3.26).  They 
are  respectively  reproduced  below.  Equation  (3.23)  is  the  general  analytical  expression,  and  Equation  (3.26) 
is  a  specific  Equation  derived  from  Equation  (3.23)  for  the  triangular  wing.  The  solutions  were  obtained  for 
nine  equally  spaced  chord  locations.  The  expression  for  the  wing  cross-sectional  area  was  provided  in  Ref. 
[15]. The  software  utilized  to  perform  the  analytical  calculations  is  called  MAPLE  (Ref  [16]). 

1  0(0)  Q(x)  1 

-  -^[QJ0)ln2x.QJx^)ln2(x^-x)] 

X  X  ^X  ^TC 

—  QJi,)ln2(%-x)dl 

471  •'0  471  •'jc 

(3.23) 

^  —(l-2x)  -  —(l-3x)ln4x(l-x)  (3.26) 

6o  K  1-X  7t  It 
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>  Q:=4*delta*sigma*xi''2*(1-xl); 


>  dQ1:=diff(Q,xi); 


>  dQ2:=diff(dQ1,xi); 


>  dQ3:=diff(dQ2,xi); 


g:=4  5a4^(l  -4) 


dQl  :=86o|(l  -^)-48aV 
dQ2  ;=86o(l-4)-16  6cy^ 


>  ans:ssubs(xb1,dQ2); 


dQ3  :=  -24  6  o 


ans  :=  -16  6  a 


>  termi  :=(1/(4*Pi))*((-4*delta*slgma)/(1-x)); 


terml  :=  - 


k(1  -x) 


>  term2:=(1/(4*Pi))*(8*delta*sigma*(ln(2*x))-1 6*delta*slgma*ln(2*(1  -x))); 

1  85oln(2jc)- 166aln(2-2x) 

terml  := - 

4  TC 


>  term3:=dQ3*ln(2*(x-xi)); 

:=  -24  6aln(2x-2^) 


>  term4:=dQ3*ln(2*(xi-x)); 


term4  :=  -24  5  a  ln(  2  ^  -  2  jc) 


>  term5:=(1/(4*Pi))*mt(term3,xi=:0.0..x); 


termS  := 


1  -24  5  a  ln(2x)  X  +  24  5  ax 


>  term6:=(1/(4*Pi))*int(term4,xl=x..1); 
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term6  :=  “  (-24  6  a  ln(2  -  2  jc)  +  24  6  a  +  24  6  a  ln(2  -  2  x)  jc 
-245ax)/7C 

>  sigma:=1/6: 

>  uoverliinf :=:term1  -term2-ternn5+term6; 

-5ln(2x)--61n(2-2x) 

1  6  13  3 

uoverUinf:=  - - - - ^ - 

67i(1-x)4  n 

1  -4  6  ln(  2  x)  X  +  4  6  X 
4  % 

1  -4  6ln(2-2x)  +  45  +  46ln(2-2x)x-45x 
4  n 
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>  xvec:=array{1  ..  10); 


xvec  :=  array(  1  ..  10,  [  ]) 
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>  for  i  from  1  to  10  do  uvel[i]:=evalf(subs(x=:xvec[i],uoverUinf))  od; 

MveZ  ;=. 27 15820284  5 

wel  :=. 1436123850  5 

\ 

uvel  :=  .05338583318  6 
Mve/,„  ;= -.02562370518  5 

[4] 

uvel  := -.1061032954  5 

L  J 

uvel  -.1997561750  3 

[6] 

uvel  :=  -.3245121960  5 
“veL,  :=  -.5225377152  5 

Lo  J 

:=  -.9694453447  5 

I"  J 


Error,  division  by  zero 


>uoverUlnf:=-(1/(1-x)/Pi)+6*(1-2*x)/Pi-2*(1-3*x)*ln(4*x)*(1-x)/Pi; 

uoverUinf  := 

1  ^\-2x  ^  (1  -  3  x)  ln(4;c)  (1  -  jc) 

- - +  6 - -  2 - 

( 1  -  x)  71  n  n 


>  xvec:=array(1 ..  10); 

xvec  :=  array(  1  ..  10,  [  ]) 

>  for  I  to  10  do  xvec[i]:=l/10  od; 


xvec^ ,  •=  — 

[1]  10 


[2] 


:■ 

[3] 


xvec^ 


xvec^ 


xvec. ..  := 


xvec._.  := 
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4 

5 


xvec^^^ : 


10 


XV^C^IO]  :=  1 

>  for  i  from  1  to  10  do  uvel[i]:=evalf(subs(x=xvec[i],uoverUinf))  od; 

uvel,-.=  1.541706944 
[1] 

uvel  :=  .7934866633 

I  ^  i 

uvel  :=  .3010904523 
uvel  :=  -.1126389811 

[4] 

uvel  ;=  -.5263019721 
uvel^^  :=  -.9993975524 

[o] 

uvel  :=  -1.608669574 
uvel^^  :=  -2.530129275 

L  o  J 

-4.572356786 

Error,  division  by  zero 


> 
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APPENDIX  C.  TRIANGULAR  WING  CROSS-SECTIONAL  FLOW 
CALCULATION 

This  appendix  shows  all  the  intermediate  calculations  leading  to  the  solution  of  the  cross-sectional 
incompressible  flow  for  the  triangular  wing  studied  in  section  G  of  chapter  IQ.  The  form  of  the  spatial 
influence  flow  equation  solved  is  introduced  as  Equation  (3.17)  and  is  reproduced  below.  The  solution  was 
obtained  for  nine  equally  spaced  chord  locations  and  the  process  repeated  for  two  spanwise  locations  (near 
the  root  and  near  the  tip).  For  conciseness,  the  intermediate  spanwise  locations  although  calculated  in  our 
study  are  not  presented  because  of  their  repetitive  character.  The  software  utilized  to  pe'rform  the  analytical 
calculations  is  called  MAPLE  (Ref.  [16]). 


o6  2it  j(x) 


P  (x^)  orv'k  I  .  P  (x,0-P  (x^)  - 
°  ln|  y  \  -  dC 

2lt  2  It  Jo  ^-2 

l+'T - 

s(x) 


(3.17) 
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>  Q:=4*delta*sigma*x'^2*(1-x); 


2  :=  4  5  o  ( 1  -  x) 


>  Qbar:=Q/(sigma*delta); 


>  dQbar1:=diff(Qbar,x); 


Qbar  :=  4  ( 1  -  x) 


dQbarl  :=8x(l  -x)-4x^ 

>  dQbar2:=:diff(dQbar1,x); 

dQbar2  :=  8  -  24  x 


>  dPbar2:=-8*zbar; 


>  dPzetabar2:=-8*zetabar; 


dPbarl  :=  -8  zbar 


dPzetabarl  :=  -8  zetabar 


>  termi  :=(dQbar2/(2*Pi))*(ln((1/6))+ln(x)+ln(sqrt(1  -(zbar/x)'^2))); 


terml  :=  - 
2 


(8-24x)  -ln(6)  +  ln(x)  +  In 


-  zbai^ 


>  term2:=(dPbar2/(2*Pi))*ln(abs((1  -(zbar/x))/(1  +(zbar/x)))); 

( I  zbar 


zbar In 


zbar 


terml  :=  -4 


>term3:=(dPzetabar2-dPbar2)/(zetabar-(zbar)); 

-8  zetabar  +  8  zbar 


termS  := 


zetabar  -  zbar 
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>term4:=:(1/Pi)*lnt(term3,zeta=:0.0..(x))," 


-8  zetahar  +  8  zhar 
zetabar  -  zhar 


term4  := 


>  term4:=(1/Pi)*int(term3,zeta=0.0..(x)); 


term4  :=  -8  — 

n 


>  sigma:=1/6: 


>  phi:=term1-term2-term4; 


(8  -  24 x)  -ln(6)  +ln(jc)  +  ln 


zbaiP- 


zbar 


zbar  In 


zbar 

1  + - 

X  J  X 

- ^  +  8- 

n  n 
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>  #  CALCULATIONS  FOR  X=1/10 

>  phianaly1:=subs(x=1/10,phi); 


phianalyl  :=  — 


zbar  In 


ill 

-ln(6)  +  ln  —  +ln( /T 
14  10 


-100  zhaP' ) 


1-10  zhar\ 

1-1-10  zbar  I  4  1 

- + - 

71  5  71 


>  zvec:=array(1  ..  10); 


zvec  :=  aiTay(  1  ..  10,  [  ]) 

>  for  i  to  10  do  zvec[i]:=6*i/600  od; 


zvec  ^  ^  := - 

[1]  100 


ZVeCj3j  :: 


100 


zveCj5]  := 
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zveCr,,  :=  — 

t6]  50 


zvec 


100 


zvec 


[8]-  25 


•=  - 

t9]  100 


^^^[10]*  10 


>  #  the  division  by  6  is  to  get  the  form  of  the  answers  which  has  no  sigma 
t  the  denominator 

>  ans:=subs(zbar=1/1 00,phianaly1 ); 


(1/6)  a 


_  -ln(6)-ln(10)  +  ln  - J99  /lOO 

14  100 ^ 


ans  := 


4  1 

+ - 

5  n 

>  final  :=:evalf(ans)/6; 


final  :=  -.5669238108 


>  ans:=subs(zbar=:1/50,phianaly1): 

>  final:=evalf(ans)/6; 
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final  :=  -.5705043033 


>  ans:=:subs(zbar=:3/1 00,phianaly1 ): 

>  final  :=evalf(ans)/6; 

final  ;=  -.5766971103 

> 

>  ans:=:subs(zbar=1/25,phianaly1): 

>  finai:=:evalf(ans)/6; 

final  :=-.5S5S9322Sl 

> 

>  ans:=subs(zbar=1/20,phianaly1 ); 


5  n  5  -  %  5  % 

>  final:=:evalf(ans)/6; 

final  :=  -.5987749621 

> 

>  ans:=:subs(zbars3/50,phianaly1 ); 

ans  := 


71  25  71  5  7C 


5 

>  final:=evalf(ans)/6; 


192 


final  :=  -.6165491328 

> 

>  ans:=:subs(zbar=7/1 00,phianaly  1 ); 


4  1 

+ - 

5  % 

>  finai:=evaIf(ans)/6; 

final  ;=  -.6415289300 

> 

>  ans:=subs(zbar=2/25,phianaly1); 

ans  := 


5  71  25  7C  5  71 


>  final:=;evalf(ans)/6; 

final  :=  -.6789331503 

> 

>  ans:=subs(zbar=9/1 00,phianaly  1 ); 
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>  final :=evalf(ans)/6; 

final  :=  -.7453324925 


> 

>  ans:=subs(zbar=1/1 0,phianaly1); 


14  -lii(6)-ln(10)  +  ln(0) 

ans  := - 

5  71 


In 


0 


V 


A 

y 


71 


4  1 

+ - 

5  K 


>  final:=evalf(ans)/6; 

Error,  (in  In)  singularity  encountered 
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>  #  CALCULATIONS  FOR  X=1 

>  phianaly10:=subs(x=1,phi); 


phianalylO  :=  -8 


•ln(6)  +  ln(  1 )  +  ln( /l  -  zhar^ 


zhar In 


1  - zhar  ^ 


1  +  zhar 


>  for  I  to  10  do  zvec[i]:=6*ri  0/600  od: 

> 

>  ans:=subs(zbar=1/1 0,phianaly1 0); 


+  8- 

n 


-ln(6)  +  In 


ans  :=  -8 


100 


99/100 


>  final:=:evalf(ans)/6; 


final  :=  1.182733902 


>ans:=:subs(zbar=1/5,phianaly10); 


ans  :=  -8 


-ln(6)  +  ln  —  .  24.  25 


5  n 


>  final:=evalf(ans)/6; 


final  1.176313741 
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> 

>  ans:=subs(zbar=3/1 0,phianaly1 0); 


>  finai:s:evalf(ans)/6; 

final  :=  1.165463606 


>  final:=evalf(ans)/6; 

final  :=  1.149937580 

> 

>  ans:=subs(zbar=1/2,phianaly10); 


K 


>  final:=evalf(ans)/6; 
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final  :=  1.129341166 


>  ans:=subs(zbar=3/5,phianaly1 0); 


ans  :=  -8 


-ln(6)  +  ln  -J16  J25 


5 

In  — 


5  % 


>  final  :=evaif(ans)/6; 


final  :=  1.103056103 


>  ans:=subs(zbar=7/10,phianaly1 0); 


-ln(6)  +  ln 


ans  :=  -8 


100 


51  /lOO 


>  final:=evalf(ans)/6; 


final  :=  1.070082027 


>  ans:=subs(zbars4/5,phianaly10); 


+  8- 
n 


ans  :=  -8 


f  1  ^ 

-ln(6)  +  ln 

_ J 


5  7C 


1 

+  8- 
n 
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>  finai:=evalf(ans)/6; 


final  :=  1.028648218 


>  ans:=subs(zbars9/1 0,phianaly1 0); 


f 


-ln(6)  +  In 


ans  :=  -8 


—Ji9jm 
100^  ^ 


In 


% 


18 

+  — 

5 


1 

10 

19 

10 

J 

n 


1 

+  8- 
K 


>  final:=:evalf(ans)/6; 


final  ■-  .9749312040 


>  ans:=8ubs{zbar^1,phianaly10); 


^  4n(6)  +  ln(0)  ^ 

ans  :=  -8 - +  4 


In 


r 

0| 

2 

J 

n 


n 


~  1 
+  8- 
K 


>  final:=evalf(ans)/6; 
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>  h:=2*delta*(x-abs(zysigma)*(1  -x); 


/!;=2  5 


(1-x) 


>  dP2bar:s:-8*zbar; 


dP2bar  :=  -8  zbar 


APPENDIX  D.  STREAMWISE  VELOCITY  DISTRIBUTION  TEST 
ON  THE  RECTANGULAR  WING 


This  appendix  is  related  to  the  discussion  in  section  G-2  of  chapter  m  about  the  velocity 
perturbation  distribution  in  the  streamwise  direction  on  the  triangular  wing.  It  was  said  for  this  case  that 
PMARC  produces  erroneous  results.  This  appendix  shows  that  for  the  rectangular  wing  the  streamwise 
velocity  distribution  seems  to  present  a  proper  behavior.  The  appendix  is  broken  down  as  follows. 

a.  Application  of  the  slender  body  theory  perturbation  expression  to  corroborate  the  results 
subsequently  obtained  with  PMARC, 

b.  Plot  showing  the  streamwise  velocity  as  a  function  of  chordwise  location  for  various 
spanwise  locations,  in  the  case  of  a  rectangular  wing  of  aspect  ratio  20  with  a  biconvex  cross-section  which 
thickness  is  tapered  from  root  to  tip. 

c.  Plot  showing  the  streamwise  velocity  as  a  function  of  chordwise  location  for  various 
spanwise  locations,  in  the  case  of  a  rectangular  wing  of  aspect  ratio  7  with  a  biconvex  cross-section  which 
thickness  is  tapered  from  root  to  tip.. 

d.  Plot  showing  the  streamwise  velocity  as  a  function  of  chordwise  location  for  various 
spanwise  locations,  in  the  case  of  a  rectangular  wing  of  aspect  ratio  4.5  with  a  biconvex  cross-section  which 
thickness  is  tapered  from  root  to  tip.. 

e.  Plot  showing  the  streamwise  velocity  as  a  function  of  chordwise  location  for  various 
spanwise  locations,  in  the  case  of  a  rectangular  wing  of  aspect  ratio  4.5  with  a  NACA  0012  cross-section 
which  thickness  is  tapered  from  root  to  tip.. 

As  can  be  observed,  the  PMARC  results  for  the  biconvex  cross-section  cases  show  similar  behavior 
than  the  one  predicted  by  the  slender  body  theory.  Furthermore,  all  the  distributions  shown  are  legitimate 
as  the  velocities  (or  perturbations)  become  smaller  toward  the  tip,  as  expected. 
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>u:=((4*tau/Pi)*(1-(0.5-x)*ln(abs((1-xyx)))*Vinf)+Vinf; 

f 

T  1  -  (.5  - x)  In -  Vinf 

u:=A  — — - +  Vinf 

n 


>  xvec:=:array(1..10); 

xvec  :=  aiTay(  1  ..  10,  [  ]) 

>  for  i  to  10  do  xvec[i]:=i/10  od; 


xvec^  — 
t*]  10 


xvec  := 

[2] 


:= 


xvec^^^ :: 


xvec^^^ 


xvec^^^:-. 
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xvec,„,  :=  — 
[9]  10 


xv<.Cj,o,:=l 


>  tau:=0.098: 

>Vinf:=11: 

>  for  i  from  1  to  10  do  uvel[i]:=evalf(subs(x=xvec[i],u))  od; 

«vgZ  :=  11.16623003 

uvel  :=  11.80172380 

uvel  ;=  12.13996012 

L  J 

MveZ  ,  :=  12.31690003 

[4] 

uvel  :=  12.37255223 

L  J 

uvel,  ;=  12.31690003 

LoJ 

uvel  :=  12.13996012 

ttveZ,;=  11.80172380 

LoJ 

Mve/  :=  11.16623003 

Error,  (in  evalf/ln)  singularity  encountered 
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Slender  Body  Theory  -  Rectangular  Wing 
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Vx  vs  X/C  at  Various  Percent  Half-Span  Locations  -  Biconvex  Set.  Thickness  Taper  Rect.  Wing,  AR=20 


o> 

d 


00 
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12.5 


(xa)  AijooiaA  asiMLueajis 
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10.5 


12.5 


_j_ 

(M 


_L 

in 


(XA)  AI!00|9A  0SjMlUB0JlS 
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10.5 


APPENDIX  E.  SPANWISE  CROSS-SECTIONAL  AREA  CODE 


This  appendix  present  the  MATLAB  computer  code  created  to  transform  the  chordwise  section 
definition  of  the  F5  wing  into  41  spanwise  cross-sectional  areas.  As  can  be  seen,  the  program 's  input 
requirements  are  the  x,  y,  and  z  coordinates  of  the  cross-sections.  The  y  coordinates  can  be  taken  from  the 
PMARC  output  file,  which  by  virtue  of  its  construction  delineate  some  of  the  geometiy  parameters.  Also 
required  in  the  input  are  the  sweep  angle  and  span.  Section  B-6-a  of  chapter  IV  details  how  the  code 
operates. 

This  code  could  be  used  on  other  wings  than  the  F5  providing  the  input  data  mentioned  above  is 
updated  as  well  as  some  other  parameters  scattered  within  the  code.  It  could  become  advantageous  to  bring 
all  variables  that  are  now  scattered  in  the  code  at  the  beginning  of  the  file  so  it  is  a  lot  easier  to  adapt  to  other 
wings.  Section  B-6-b  of  chapter  IV  explains  why  the  program  is  presently  geared  more  towards  the  F5 
wing. 
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Printed  by  carrier  from  viper 
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Printed  by  carrier  from  viper 
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APPENDIX  F.  F5  WING  SPANWISE  CROSS-SECTIONAL  AREAS 


This  appendix  provides  the  results  of  the  computer  code  described  in  section  A-6  of  chapter  IV.  It 
shows  the  area,  radius  and  diameter  of  each  of  the  41  spanwise  cross-sections  defined  along  the  chord  of  the 
F5  wing.  Also,  a  plot  is  included  that  shows  the  distribution  of  the  spanwise  cross-sectional  areas  along  the 
chord. 
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Columns  22  through  28 


APPENDIX  G.  TRIANGULAR  WING  SPANWISE  CROSS-SECTIONAL 
AREAS 


This  appendix  provides  the  radii  of  46  spanwise  cross-sectional  areas  with  the  x  coordinate  of  their 
centroid  along  the  chord  of  the  triangular  wing.  This  information  is  provided  for  the  two  cases  of  maximum 
thickness  ratio  treated:  0.098,  and  0.196.  A  plot  showing  the  distribution  of  the  spanwise  cross-sectional 
areas  along  the  chord  is  also  included  for  each  case. 
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Spanwise  Cross-sectional  Are 
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Columns  43  through  46 

0.0224  0.0160  0.0113 


XandR_lril96BR 

0.0147  0,0193  0.0238 

0.0442  0.0478  0.0512 

0.0657  0.0681  0.0702 

0.0774  0.0781  0.0784 

0.0750  0.0730  0.0704 

0.0433  0.0380  0.0314 

0 


Columns  1  through  7 

0  0.0250  0,0500 

Columns  8  through  14 

0.1750  0.2000  0,2250 

Columns  15  through  21 

0:3500  0.3750  0.4000 

Columns  22  through  28 

0,5250  0.5500  0,5750 

Columns  29  through  35 

0.7000  0.7250  0.7500 

Columns  36  through  42 

0.8750  0.9000  0.9250 

Columns  43  through  46 

0.9876  0.9938  0.9969 


0.0750  0.1000  0.1250 

0.2500  0.2750  0.3000 

0.4250  0.4500  0,4750 

0.6000  0.6250  0.6500 

0.7750  0.8000  0.8250 

0.9500  0.9625  0.9750 

1.0000 


p||ii 


0.0282 

0.0545 

0.0721 

0.0785 

0.0671 

0.0274 


0.1500 


0,3250 


0.5000 


0.6750 


0.8500 


0.9813 
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Spanwise  Cross-sectional  Area 


Triangular  Wing  AR=!2/3,  delta=0.196,  Q  vs  x 
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APPENDIX  H.  AUTOMATED  EXTRACTION  OF  CL  AND  CM 


This  appendix  presents  two  short  programs  written  in  the  UNIX  AWK  language  used  to  extract  CL 
and  CM  at  each  time-step  of  the  PMARC  output  file  'DATA6’.  These  programs  introduced  in  section  B-4 
of  chapter  V  are  specific  to  the  case  of  the  F5  wing  treated  in  this  study.  The  location  of  CL  and  CM  in  the 
output  file  can  vary  depending  on  the  configuration  modelled. 
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CM  EXTRACTION 


# ! /bin/sh 

#  AWK  Script  to  extract  the  total  cm  value  from 

#  the  PMARC  output  DATA6  file. 

#  Finds  "WIND  AXES"  and  skips  32  lines  then  picks  the  fourth  field  ($4 

#  which  is  the  total  cm  value 

# 

echo 

echo  Reads  in  DATA6  file  for  a  simple  wing  with  symmetry 

echo  and  its  wingtips  for  a  total  of  four  patches 

echo 

echo  cm  total  data  is  put  in  file  \’'outcm\* 
awk  ' 

/WIND  AXES/  { 

NL  =  32 

f or { i=0 ; i<NL; ++i ) 

getline  #  Gets  32  line  after  WIND  AXES  term 

print  $4  #  Gets  fourth  field  (cm  value) 

} 

'  DATA6  >  outcm 


CL  EXTRACTION 

***★*•*■★★★*★★*  ^ 

# i /bin/sh 

#  AWK  Script  to  extract  the  total  CL  value  from 

#  the  PMARC  output  DATA6  file. 

#  Finds  "WIND  AXES"  and  skips  32  lines  then  picks  first  field  ($1) 

#  which  is  the  total  CL  value 

# 

echo 

echo  Reads  in  DATA6  file  for  a  simple  wing  with  symmetry 

echo  and  its  wingtips  for  a  total  of  four  patches 

echo 

echo  CL  total  data  is  put  in  file  \"outCL\* 
awk  ' 

/WIND  AXES/  { 

NL  =  32 

for (i=0; i<NL;++i) 

getline  #  Gets  32  line  after  WIND  AXES  term 

print  $1  #  Gets  first  field  (CL  value) 

} 

'  DATA6  >  out CL 
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